BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,925	T→A	27.5%	intergenic (‑1008/‑160)	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,925	0	T	A	27.5%	275.7 / 127.9	156	intergenic (‑1008/‑160)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (63/50); new base A (24/19); total (87/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 163 total aligned reads displayed.

GCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCT  >  CP000730/1929785‑1930064

gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                                                                                                                                             <  2:368184/141‑1 (MQ=255)

gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                                                                                                                                             <  2:1448315/141‑1 (MQ=255)

    gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                                                                        >  1:325601/12‑138 (MQ=255)

     tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                                                                                        <  1:611918/141‑1 (MQ=255)

        ttttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCCATTGCGTTTTTTTcgt                                                                                                                                     <  2:895126/140‑3 (MQ=255)

           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                 >  2:419108/5‑141 (MQ=255)

           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑A‑TTGCTAATCATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                 >  2:656539/5‑141 (MQ=255)

           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAACACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                 >  2:518595/5‑141 (MQ=255)

           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGCAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                 <  1:782490/137‑1 (MQ=255)

             ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                >  1:1519713/1‑141 (MQ=255)

             cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                               <  1:1619207/139‑1 (MQ=255)

             cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                               <  1:1543040/139‑1 (MQ=255)

             cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                               >  1:1257520/3‑141 (MQ=255)

              agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                              >  1:303314/2‑141 (MQ=255)

               gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                                                                                              <  1:172675/141‑1 (MQ=255)

               gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                                                                                              <  2:199270/141‑1 (MQ=255)

                 gggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACT‑TTGCTAATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                           >  1:1569239/4‑138 (MQ=255)

                   atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                          >  1:427389/1‑141 (MQ=255)

                    tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGGAGAGAGGAGGT‑GTA‑TGA‑G‑GTGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                         <  2:319493/141‑1 (MQ=255)

                    tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                        <  2:260133/141‑2 (MQ=255)

                        ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                    <  2:939624/137‑1 (MQ=255)

                        ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                    <  2:1369602/137‑1 (MQ=255)

                        ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                                     >  2:816708/1‑141 (MQ=255)

                          agtagAGTCGCCTATCTCGCAGGCGTCAAATTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa                                                                                                                   >  2:473945/1‑141 (MQ=255)

                           atagAGTCGCCTATCTCTCAGGCGTCAAGT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                 <  2:1651970/140‑2 (MQ=255)

                            tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                >  2:907281/1‑139 (MQ=255)

                             agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                                >  1:1503160/1‑141 (MQ=255)

                             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                               >  1:493091/1‑138 (MQ=255)

                              gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                               <  1:1643263/141‑1 (MQ=255)

                               agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                             <  1:1173820/141‑6 (MQ=255)

                               agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                             <  1:333207/141‑6 (MQ=255)

                                gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                                                                             >  2:875771/1‑141 (MQ=255)

                                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                            <  2:622170/141‑7 (MQ=255)

                                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                            <  1:1116210/141‑7 (MQ=255)

                                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                            >  1:569305/1‑135 (MQ=255)

                                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGATGGGG‑TAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                            >  1:798539/1‑135 (MQ=255)

                                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                          >  1:476466/1‑133 (MQ=255)

                                   gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGCGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                                                                          >  2:1102600/1‑141 (MQ=255)

                                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                         <  1:107660/141‑10 (MQ=255)

                                     cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAAAg                                                                                                        >  1:828419/1‑141 (MQ=255)

                                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                       >  1:1433287/1‑130 (MQ=255)

                                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                       <  2:1009074/141‑12 (MQ=255)

                                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaataattaa                                                                                                       >  1:389714/1‑130 (MQ=255)

                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaataattaata                                                                                                     <  2:389714/141‑14 (MQ=255)

                                        tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGACGCAATACAAAATAGGTg                                                                                                     >  1:627953/1‑141 (MQ=255)

                                              aGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                                                                               <  1:1058421/141‑1 (MQ=255)

                                                   tCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACACAATAGGTAACATATACCCg                                                                                          >  2:593632/1‑141 (MQ=255)

                                                    cAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                                                                                         >  1:234803/1‑141 (MQ=255)

                                                     aaTTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGCCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                                                                                        <  2:1303011/141‑1 (MQ=255)

                                                        ttAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                                                                     >  2:785120/1‑141 (MQ=255)

                                                         tAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCaa                                                                                    >  2:307028/1‑141 (MQ=255)

                                                          aGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                                   >  2:1034142/1‑141 (MQ=255)

                                                          aGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                                   <  2:184202/141‑1 (MQ=255)

                                                          aGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                                   <  2:1503160/141‑1 (MQ=255)

                                                           gACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATa                                                                                  >  2:398111/1‑141 (MQ=255)

                                                             cGCAGAGAGGAGGT‑GTA‑TAA‑G‑TTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaa                                                                                <  2:148404/141‑1 (MQ=255)

                                                              gCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                               >  1:71365/1‑141 (MQ=255)

                                                              gCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                               <  1:1004352/141‑1 (MQ=255)

                                                              gCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                               >  1:1089195/1‑141 (MQ=255)

                                                              gCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                               >  1:71364/1‑141 (MQ=255)

                                                              gCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTAATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                               >  2:1318084/1‑141 (MQ=255)

                                                                     ggaggT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                        >  2:1583716/1‑141 (MQ=255)

                                                                       aggT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                                                                      >  2:767599/1‑141 (MQ=255)

                                                                       aggT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                                                                      >  1:1449911/1‑141 (MQ=255)

                                                                                 aa‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa                                                              >  1:482454/1‑141 (MQ=255)

                                                                                 aa‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCATCACTACCGCa                                                              >  2:1164721/1‑141 (MQ=255)

                                                                                    g‑gTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                            <  1:1002215/141‑1 (MQ=255)

                                                                                    g‑gTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                            <  2:69780/141‑1 (MQ=255)

                                                                                    g‑gTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                            <  2:273451/141‑1 (MQ=255)

                                                                                      gTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAt                                                           >  1:666623/1‑141 (MQ=255)

                                                                                            cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                                                                           >  2:1587018/1‑87 (MQ=255)

                                                                                            cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                                                                           <  1:1587018/87‑1 (MQ=255)

                                                                                            cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                     <  1:24189/141‑1 (MQ=255)

                                                                                              tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                   <  1:633162/141‑1 (MQ=255)

                                                                                              tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                   <  2:286852/141‑1 (MQ=255)

                                                                                               aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGCCGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACAGCAAATAGTGAgggg                                                  >  2:1396390/1‑141 (MQ=255)

                                                                                                ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                 <  1:1315681/141‑1 (MQ=255)

                                                                                                ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                 <  2:1357559/141‑1 (MQ=255)

                                                                                                      ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTTGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                           >  2:843509/1‑141 (MQ=255)

                                                                                                          catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                       >  1:1042125/1‑141 (MQ=255)

                                                                                                            tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAg                                     <  1:1484012/141‑1 (MQ=255)

                                                                                                               tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                  >  2:551125/1‑141 (MQ=255)

                                                                                                                  cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                 <  2:1492794/59‑1 (MQ=255)

                                                                                                                  cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                 >  1:1492794/1‑59 (MQ=255)

                                                                                                                  cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                 >  1:1492543/1‑59 (MQ=255)

                                                                                                                  cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                               >  2:1654052/1‑141 (MQ=255)

                                                                                                                   aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATATTGAGGGGATTGGTGTATAAGTAAATAc                              <  1:386347/141‑1 (MQ=255)

                                                                                                                   aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                              <  2:627953/141‑1 (MQ=255)

                                                                                                                   aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                              <  1:781169/141‑1 (MQ=255)

                                                                                                                     cAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt                            >  2:1303069/1‑141 (MQ=255)

                                                                                                                          atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                       <  1:759162/141‑1 (MQ=255)

                                                                                                                           tccgTGGCTGTGCCATTGCGTTTTTTTCTTATTTGTTTAGTAGACGCGATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATTGTGAGGGAGATGGTGTAAAAGTATGTACGTAtgtgt                      >  2:745494/4‑137 (MQ=255)

                                                                                                                              gTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTTAGTAGACGCAATACAAAATAGGTGACATATAGCCGCCCCAATAACAATCCCCTCACTACCGCAAATAGTTAGGGGATTGGTGTATACTTAAATTCTTTTTTTCGtt                   >  2:845157/1‑141 (MQ=255)

                                                                                                                                 gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGCGGATTGTTGTATAAGTAAATACTTATTTTCGTTGTc                <  1:793150/141‑1 (MQ=255)

                                                                                                                                    gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa             >  2:769257/1‑141 (MQ=255)

                                                                                                                                     tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc                                                                                            <  2:1185039/61‑1 (MQ=255)

                                                                                                                                     tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTaa            >  1:182714/1‑141 (MQ=255)

                                                                                                                                      gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt           >  2:822892/1‑141 (MQ=255)

                                                                                                                                             cGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg    <  2:1066434/141‑1 (MQ=255)

                                                                                                                                               tttttttCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCt  <  2:350340/141‑1 (MQ=255)

GCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCT  >  CP000730/1929785‑1930064

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