BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,941	A→T	22.9%	intergenic (‑1024/‑144)	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,941	0	A	T	22.9%	306.5 / 91.3	147	intergenic (‑1024/‑144)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (65/48); new base T (18/16); total (83/64)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 150 total aligned reads displayed.

GATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAA  >  CP000730/1929803‑1930081

 gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                                                 >  2:863701/2‑141 (MQ=255)

 atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                                                  >  1:562697/1‑141 (MQ=255)

 atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                                                  >  1:427389/1‑141 (MQ=255)

  tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGGAGAGAGGAGGTGTA‑TGA‑G‑G‑TGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                                                 <  2:319493/141‑1 (MQ=255)

  tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                                <  2:260133/141‑2 (MQ=255)

     aataatAGAGACGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                             >  2:578323/1‑141 (MQ=255)

      ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                            <  2:1101199/137‑1 (MQ=255)

      ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                            <  2:939624/137‑1 (MQ=255)

      ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                                                             >  2:816708/1‑141 (MQ=255)

       tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGc                                                                                                                                            >  2:650653/1‑141 (MQ=255)

       taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                           >  2:827990/4‑141 (MQ=255)

        agtagAGTCGCCTATCTCGCAGGCGTCAAATTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa                                                                                                                                           >  2:473945/1‑141 (MQ=255)

          tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                                        >  2:907281/1‑139 (MQ=255)

           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                       >  1:483866/1‑138 (MQ=255)

             agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                     <  1:333207/141‑6 (MQ=255)

             agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                     >  1:659176/1‑136 (MQ=255)

             agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                     <  1:1173820/141‑6 (MQ=255)

              gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa                                                                                                                                   <  1:1601841/141‑1 (MQ=255)

              gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                                    >  1:569305/1‑135 (MQ=255)

               tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑a                                                                                                                                <  1:1464183/141‑1 (MQ=255)

                 gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAt                                                                                                                              >  2:1261408/1‑141 (MQ=255)

                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                 <  2:377700/141‑10 (MQ=255)

                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                 <  1:1499478/141‑10 (MQ=255)

                   cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAAAg                                                                                                                            >  1:828419/1‑141 (MQ=255)

                   cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                               <  2:1009074/141‑12 (MQ=255)

                   cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaataattaa                                                                                                                               >  1:389714/1‑130 (MQ=255)

                   cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                               <  2:434342/141‑12 (MQ=255)

                      tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGACGCAAT‑‑ACAA‑‑AATAGGTg                                                                                                                         >  1:627953/1‑141 (MQ=255)

                            aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑atat                                                                                                                  <  1:1058421/141‑1 (MQ=255)

                              gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑Ag                                                                                                               >  2:398026/1‑141 (MQ=255)

                                  cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gc                                                                                                          >  1:234803/1‑141 (MQ=255)

                                      ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCa                                                                                                      >  2:785120/1‑141 (MQ=255)

                                        aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATGT‑AGCC‑GCACCAAt                                                                                                    >  2:136528/1‑141 (MQ=255)

                                        aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                    <  1:1087759/141‑1 (MQ=255)

                                        aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                    <  2:1503160/141‑1 (MQ=255)

                                        aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                    >  2:1034142/1‑141 (MQ=255)

                                        aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                    <  2:184202/141‑1 (MQ=255)

                                         gACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATa                                                                                                   <  2:219729/141‑1 (MQ=255)

                                         gACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATa                                                                                                   >  2:398111/1‑141 (MQ=255)

                                           cGCAGAGAGGAGGTGTA‑TAA‑G‑T‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaa                                                                                                 <  2:148404/141‑1 (MQ=255)

                                            gCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                >  1:71365/1‑141 (MQ=255)

                                            gCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                >  1:71364/1‑141 (MQ=255)

                                            gCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                <  1:1004352/141‑1 (MQ=255)

                                            gCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                >  1:1089195/1‑141 (MQ=255)

                                            gCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTAATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                >  2:1318084/1‑141 (MQ=255)

                                                     aggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCa                                                                                       >  2:767599/1‑141 (MQ=255)

                                                     aggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCa                                                                                       >  1:1449911/1‑141 (MQ=255)

                                                              aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                               >  1:482454/1‑141 (MQ=255)

                                                              aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCATCACTACCGCa                                                                               >  2:1164721/1‑141 (MQ=255)

                                                                 g‑g‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTatttgcgg                                                                             <  1:433283/141‑8 (MQ=255)

                                                                 g‑g‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                                             <  1:1002215/141‑1 (MQ=255)

                                                                   g‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTatttgcggt                                                                            <  2:310055/141‑9 (MQ=255)

                                                                   g‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAAt                                                                            >  1:666623/1‑141 (MQ=255)

                                                                       aTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCCCTACCGCAAATAGt                                                                         >  2:1552005/1‑141 (MQ=255)

                                                                          cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aa                                                                                                                               >  2:1587018/1‑87 (MQ=255)

                                                                          cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                      <  1:24189/141‑1 (MQ=255)

                                                                            tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                    <  2:286852/141‑1 (MQ=255)

                                                                             aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGCCGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACAGCAAATAGTGAgggg                                                                   >  2:1396390/1‑141 (MQ=255)

                                                                             aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                   >  2:276834/1‑141 (MQ=255)

                                                                              ttttCGTTCGCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                  <  2:893844/141‑1 (MQ=255)

                                                                              ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTATCGCAAATAGTGAGGGGa                                                                  <  2:1382912/141‑1 (MQ=255)

                                                                              ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                  <  2:1357559/141‑1 (MQ=255)

                                                                              ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                  >  1:1374596/1‑141 (MQ=255)

                                                                                ttCGTTCACATCATAGCACCAGTAATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                                                <  1:180279/141‑1 (MQ=255)

                                                                                 tCGTTCACATCATAGCCCCATCCATCCGTGGTTTTGCCATGTCGGTGGTTTTTTATTTTCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTg                                                               <  1:1562916/141‑1 (MQ=255)

                                                                                    ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                                            >  1:76360/1‑141 (MQ=255)

                                                                                     tCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtgt                                                           <  1:1504494/141‑1 (MQ=255)

                                                                                        catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                                        <  1:307028/141‑1 (MQ=255)

                                                                                        catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                                        >  1:1042125/1‑141 (MQ=255)

                                                                                         atcatAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGCGGGGATTGgtg                                                            >  2:1641040/1‑141 (MQ=255)

                                                                                          tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAg                                                      >  2:743996/1‑141 (MQ=255)

                                                                                             tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                                   >  1:1625092/1‑141 (MQ=255)

                                                                                             tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                                   >  2:551125/1‑141 (MQ=255)

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

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

                                                                                                cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                                                >  2:1654052/1‑141 (MQ=255)

                                                                                                 aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATATTGAGGGGATTGGTGTATAAGTAAATAc                                               <  1:386347/141‑1 (MQ=255)

                                                                                                 aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                               <  2:627953/141‑1 (MQ=255)

                                                                                                  ccAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACt                                              <  1:785120/141‑1 (MQ=255)

                                                                                                   cAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt                                             >  2:1303069/1‑141 (MQ=255)

                                                                                                     gTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt                                           >  1:820262/1‑141 (MQ=255)

                                                                                                        atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                                        <  1:759162/141‑1 (MQ=255)

                                                                                                         tccgTGGCTGTGCCATTGCGTTTTTTTCTTATTTGTTTAGTAGACGCGAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATTGTGAGGGAGATGGTGTAAAAGTATGTACGTAtgtgt                                       >  2:745494/4‑137 (MQ=255)

                                                                                                            gTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTTAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCCCCAATAACAATCCCCTCACTACCGCAAATAGTTAGGGGATTGGTGTATACTTAAATTCTTTTTTTCGtt                                    >  2:845157/1‑141 (MQ=255)

                                                                                                              ggCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGt                                  >  1:1484354/1‑141 (MQ=255)

                                                                                                               gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGCGGATTGTTGTATAAGTAAATACTTATTTTCGTTGTc                                 <  1:793150/141‑1 (MQ=255)

                                                                                                                  gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                              >  2:769257/1‑141 (MQ=255)

                                                                                                                   tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGc                                                                                                              <  2:1185039/61‑1 (MQ=255)

                                                                                                                   tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGc                                                                                                              >  1:1185039/1‑61 (MQ=255)

                                                                                                                   tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTaa                             >  1:182714/1‑141 (MQ=255)

                                                                                                                    gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt                            >  2:822892/1‑141 (MQ=255)

                                                                                                                           cGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg                     <  2:1066434/141‑1 (MQ=255)

                                                                                                                             tttttttCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCt                   <  2:350340/141‑1 (MQ=255)

                                                                                                                                     ttATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTc           >  1:624558/1‑141 (MQ=255)

                                                                                                                                     ttATTGGCAAAGTAGACGCAAT‑‑ACCA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATACCCTTCCTACCGCAAATAGGGAGGGGGTCGGTGTATAAGTGAATACTTATTTTCGTTGTGTTCACTCTACTGCttcttcgta           >  2:93127/1‑133 (MQ=255)

                                                                                                                                       tttGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtt         <  2:1127000/140‑1 (MQ=255)

                                                                                                                                       aTTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtt         >  1:1488464/1‑141 (MQ=255)

                                                                                                                                           gCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTg     >  2:615529/1‑141 (MQ=255)

                                                                                                                                              aaGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAttaaaatttttataa  <  2:562697/141‑14 (MQ=255)

                                                                                                                                              aaGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaa  >  1:1345360/1‑141 (MQ=255)

GATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAA  >  CP000730/1929803‑1930081

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: `‑`