BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,930	T→A	59.3%	intergenic (‑1013/‑155)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,930	0	T	A	59.3%	86.4 / 159.1	134	intergenic (‑1013/‑155)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (27/27); new base A (50/29); total (77/57)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 138 total aligned reads displayed.

TATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTT  >  CP000730/1929792‑1930070

tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                                                                                                                     >  1:1491/9‑141 (MQ=255)

 attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                                                                                                    >  2:170795/8‑140 (MQ=255)

 aattaccaatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                                                                                                     <  1:333251/133‑1 (MQ=255)

  tttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCATCAATTTAGACGCAGAGACACGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCACCAGTGACTGTTCCCTTGCGtttttttactt                                                                                                                                                    >  2:174548/1‑137 (MQ=255)

  ttcacagatgatATAAAAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                                   <  2:23651/135‑1 (MQ=255)

   tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                                                  >  1:736519/6‑141 (MQ=255)

     acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                                >  1:687411/4‑141 (MQ=255)

     acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                                <  2:585235/138‑1 (MQ=255)

     acagatgggataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                                <  1:71719/132‑1 (MQ=255)

       agttggtataATAATAGAGTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                              <  1:639072/135‑1 (MQ=255)

       agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                              >  2:211307/2‑141 (MQ=255)

        gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                             <  1:165398/141‑2 (MQ=255)

         atggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                                                                            <  2:753092/141‑3 (MQ=255)

          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCC‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                           >  1:297551/4‑141 (MQ=255)

          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                           >  1:324897/4‑141 (MQ=255)

          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                                           <  2:425958/138‑4 (MQ=255)

          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                           >  1:430400/4‑141 (MQ=255)

          tggtataATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGCAGAGAGGAGGAGCATT‑AC‑T‑T‑TGCTAATCATTTTAGTTCACATAATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                           >  1:646422/4‑141 (MQ=255)

          gggtataATAAGAGAGTCGCCTATCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCGAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                                           <  1:261843/138‑4 (MQ=255)

           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                          >  1:409270/3‑141 (MQ=255)

           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                          >  1:582952/3‑141 (MQ=255)

           ggtataAAAATAGAGTCGCCAATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCGCATCATAGCACCAGTCATCAGTGGCTGTGCTGTTACGTATTATACTTATTGGCTTAGt                                                                                                                                          >  1:488767/3‑141 (MQ=255)

             tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                                         >  2:697637/1‑141 (MQ=255)

             tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                        >  2:668643/1‑140 (MQ=255)

             tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                        >  2:579148/1‑140 (MQ=255)

                 ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTATACGc                                                                                                                                    <  1:525159/141‑1 (MQ=255)

                 ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                    <  1:293933/141‑1 (MQ=255)

                 agtatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                    <  1:632896/137‑1 (MQ=255)

                  taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                   <  1:416749/138‑1 (MQ=255)

                   aatagATTCGCCTATCTCTCAGGCGTCAA‑‑TGGACTAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑a                                                                                                                                <  2:770001/139‑1 (MQ=255)

                   aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑a                                                                                                                                <  2:749139/139‑1 (MQ=255)

                    atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                               >  2:565761/2‑140 (MQ=255)

                     tagAGTCGCCTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aca                                                                                                                              <  2:264399/141‑3 (MQ=255)

                     atgagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aca                                                                                                                              <  1:544764/139‑3 (MQ=255)

                      agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑a                                                                                                                            <  2:478470/141‑1 (MQ=255)

                      agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTACACGCa‑‑acaa                                                                                                                             >  2:416835/1‑138 (MQ=255)

                      agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                             <  1:244552/141‑4 (MQ=255)

                       gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                            >  1:59597/1‑137 (MQ=255)

                        agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaat                                                                                                                           >  2:72282/1‑136 (MQ=255)

                        agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaat                                                                                                                           >  2:480994/1‑136 (MQ=255)

                         gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACaa                                                                                                                         <  2:32723/141‑1 (MQ=255)

                         gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACaa                                                                                                                         <  1:72687/141‑1 (MQ=255)

                         gTCGCCTATCTCTCAGGCGTCAATT‑‑GACCAAGCGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                          <  1:315589/141‑7 (MQ=255)

                          tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataa                                                                                                                         >  2:113229/1‑134 (MQ=255)

                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTACACGCa‑‑acaaatgatt                                                                                                                       >  1:244536/1‑132 (MQ=255)

                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                       <  2:154858/141‑10 (MQ=255)

                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                       >  1:287339/1‑132 (MQ=255)

                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                       <  2:231325/141‑10 (MQ=255)

                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                       <  2:155011/141‑10 (MQ=255)

                            gCCTATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                       <  2:284388/141‑10 (MQ=255)

                             ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATa                                                                                                                   >  2:193899/1‑141 (MQ=255)

                             ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCC‑TCAC‑G‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataataa                                                                                                                      >  1:730734/1‑131 (MQ=255)

                             ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                      >  2:514869/1‑131 (MQ=255)

                             ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                      >  1:491941/1‑131 (MQ=255)

                             ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                      >  2:616714/1‑131 (MQ=255)

                             ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                      >  1:457055/1‑131 (MQ=255)

                              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                     >  2:85138/1‑130 (MQ=255)

                              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                     >  2:101946/1‑130 (MQ=255)

                              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                     >  1:320741/1‑130 (MQ=255)

                              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                     >  1:744668/1‑130 (MQ=255)

                               tATCTCTCAGGCGTCAATT‑‑GTCGAAGAGAGGAGGTGCA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaat                                                                                                                    <  2:71437/141‑13 (MQ=255)

                               tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaat                                                                                                                    <  2:613140/141‑13 (MQ=255)

                               tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaaa                                                                                                                    >  1:694320/1‑129 (MQ=255)

                                aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataagtaata                                                                                                                   >  1:171442/1‑128 (MQ=255)

                                aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                   >  2:666719/1‑128 (MQ=255)

                                 tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTg                                                                                                               >  2:205758/1‑141 (MQ=255)

                                     tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑a‑t                                                                                                         >  1:608235/1‑141 (MQ=255)

                                     tcAGGCGTCAATTTAGACGCAGAGAGGAGGAGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTAGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑CATGGGTGGC‑a‑t                                                                                                         >  2:256390/1‑141 (MQ=255)

                                     ccAGGCGGCAATTTAGACGCAGAGAGGAGGGGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAGTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑a‑t                                                                                                         <  1:342103/140‑1 (MQ=255)

                                       aGGCGTCAATTTAGACGCAGAGAGGAGGGGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGGGGC‑a‑tat                                                                                                       >  2:587540/1‑141 (MQ=255)

                                            tCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑CGCC‑g                                                                                                >  2:644272/1‑141 (MQ=255)

                                                     aCGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTGTTATTCGGCCTGCCTAAATCGGCGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCCCCAAGaa                                                                                       >  2:316106/1‑141 (MQ=255)

                                                       gctgagcgGAGGTGTA‑T‑AA‑G‑G‑GGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATaaaa                                                                                     <  2:379884/138‑1 (MQ=255)

                                                        cAGAGCGGAACTGTG‑T‑AA‑G‑C‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATaaaaa                                                                                    <  1:756816/141‑1 (MQ=255)

                                                             aggaggTGTG‑G‑AA‑G‑G‑TGATGCTTAGTTTCGTTCACATCATAGCACCCGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑CCAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATcctc                                                                               >  1:566864/1‑139 (MQ=255)

                                                               gaggTGTA‑T‑AA‑G‑G‑TGATGATAGGGTGCTGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTc                                                                             >  2:748512/1‑141 (MQ=255)

                                                                   tgtA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTa                                                                         <  2:246895/141‑1 (MQ=255)

                                                                       tt‑aA‑G‑G‑TCGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCg                                                                      <  2:539066/140‑1 (MQ=255)

                                                                           a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                                   >  1:489336/1‑141 (MQ=255)

                                                                                gttaTGCTTGTTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAAt                                                                 <  2:638656/138‑1 (MQ=255)

                                                                                    tGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTATTTCTAATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAgag                                                             >  1:658523/1‑139 (MQ=255)

                                                                                       ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAACAGTGAgg                                                          >  1:49664/1‑141 (MQ=255)

                                                                                        tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                         <  2:370522/141‑1 (MQ=255)

                                                                                              cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                                   >  1:237177/1‑141 (MQ=255)

                                                                                              cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACCCA‑‑AA‑‑ACAA‑‑AATCGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAAGGGATTgg                                                   >  1:626273/1‑141 (MQ=255)

                                                                                                ttCACATCATAGCACCAGTCATCAGTGGCGGTGCGATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGGGAC‑A‑TAC‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGgtg                                                 <  2:268565/141‑1 (MQ=255)

                                                                                                    catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                             >  2:291179/1‑141 (MQ=255)

                                                                                                        atAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAA‑‑ATAATTAATAGGCGACTT‑TATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGAGGGGATTGGTGcat                                              >  1:484282/1‑141 (MQ=255)

                                                                                                         tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCA‑‑ACAAATAATTAATAGGCGAC‑TTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcatt                                             >  2:677055/1‑140 (MQ=255)

                                                                                                            cACCGGTCATCAGTGGCTGTGCCATGGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAACTTCCCCTCACTTCCGCCAATAGTGAGGGGATTGGTGTATCAGTACATa                                     <  2:377276/141‑1 (MQ=255)

                                                                                                                aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑ACTAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGGTGTCTAAGTAAATACTTa                                 >  2:118788/1‑141 (MQ=255)

                                                                                                                  tcatcaGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTAAGTAGACGCA‑‑AT‑‑GCCA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCATTACCGCAAATAGTGAGGGGATTGTTGTATAAGTAAATACTTAtt                               <  2:423686/141‑1 (MQ=255)

                                                                                                                  tcatcaGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                               <  2:386619/141‑1 (MQ=255)

                                                                                                                  gcatcaGTGGCTGTGCCATTGCGTTTTTTGCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                               <  1:61050/140‑1 (MQ=255)

                                                                                                                   catcaGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                              >  2:302004/1‑141 (MQ=255)

                                                                                                                        gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGtt                         <  1:452381/141‑1 (MQ=255)

                                                                                                                               tgCCATTGCGTTTTTTTCTAATTTGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCT‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATGCGTAAATCCTTATTTTGGTTGTCTTaa                  >  1:409666/1‑141 (MQ=38)

                                                                                                                                       cGTTGTTTTCTTATTGGATAATTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCATAAATAAACATCCCCTCACTACCGGAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg          <  2:392736/141‑1 (MQ=255)

                                                                                                                                           tttttCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa      <  1:455407/141‑1 (MQ=255)

                                                                                                                                             tttCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt    <  1:26577/141‑1 (MQ=255)

TATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTT  >  CP000730/1929792‑1930070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`