BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,917	C→T	70.0%	intergenic (‑1000/‑168)	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,917	0	C	T	70.0%	152.6 / 96.1	129	intergenic (‑1000/‑168)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (17/22); new base T (52/38); total (69/60)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 130 total aligned reads displayed.

GTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTT  >  CP000730/1929777‑1930052

gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                                                                                                                                         >  2:358101/1‑141 (MQ=255)

  gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑GTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                                                                       >  1:553498/1‑141 (MQ=255)

   aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                                                                                      <  1:87817/141‑1 (MQ=255)

    cTGTGCTTGTATTTTTGCCGATGAAATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                                                                                                                                     <  1:626350/141‑1 (MQ=255)

    cTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                                                                                                                                     <  1:408338/141‑1 (MQ=255)

    atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCCTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGCCTGTGCTAGTGc                                                                                                                                    >  2:190326/4‑141 (MQ=255)

      atgCTTGTATATTCACAGATGGCATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                  <  2:653165/140‑1 (MQ=255)

         cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                               >  2:463506/1‑141 (MQ=255)

          ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                              <  2:356769/141‑5 (MQ=255)

          ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                              >  1:638863/1‑140 (MQ=255)

          ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                              <  1:538881/141‑5 (MQ=255)

          ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TCCTAATC‑ATTTTGGTTCACATCATAGCACCAGTCTTCAGTGGCTGTGCTGTTGCGtatta                                                                                                                              >  2:282518/1‑140 (MQ=255)

          ttGTATATTCACAGATGGTATAACAATAGAGTCTCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                              <  1:102810/141‑5 (MQ=255)

           tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGTGTA‑TAT‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                                                              <  1:424911/141‑1 (MQ=255)

           tGTATATTCACAGCTGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                                                             <  1:278405/141‑6 (MQ=255)

            gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                                                            >  1:242208/12‑138 (MQ=255)

            gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGttttata                                                                                                                            >  1:575104/12‑138 (MQ=255)

              atattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                                                                          >  1:572382/10‑136 (MQ=255)

               tattcacagatggtATAATAATGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                                                                         <  2:274578/133‑7 (MQ=255)

               tattatcagatggtATAATAATAGTGTTTGCTATCTCTCGGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                                                                         <  1:120571/133‑7 (MQ=255)

                 ttcacagatggtATAATGATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                       >  1:469921/7‑141 (MQ=255)

                 ttcacagatggtATAATAATAGAGTCGCGTATCTCTCAGGCGTCAATT‑‑GACGAAAAGAGGAGGTGCATTAC‑T‑T‑TGCTAAGC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                       >  1:429667/7‑141 (MQ=255)

                  tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                      >  2:91718/6‑141 (MQ=255)

                   cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTCTCCTTCTTg                                                                                                                     >  1:458604/5‑141 (MQ=255)

                   cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                     >  2:609170/5‑141 (MQ=255)

                   cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTACACTTATTg                                                                                                                     >  2:534203/5‑141 (MQ=255)

                     ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                    <  2:161790/141‑1 (MQ=255)

                     ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                    <  1:448480/141‑1 (MQ=255)

                     caggtggtataATAATAGAGTCGCCTATCTCTCATGCGTCAATT‑‑GACGGCGAGAGGGGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTTCGTATTATACTTATTGGc                                                                                                                   <  2:285528/134‑1 (MQ=255)

                      agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                  <  1:437605/140‑1 (MQ=255)

                       tatggtATAATAATGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCTTAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                 <  2:320702/140‑2 (MQ=255)

                       gatggtATAATAATAGAGTCGCCTATCTCTCAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                 <  2:536543/141‑2 (MQ=255)

                       gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                 >  2:272694/1‑140 (MQ=255)

                       gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                 >  2:508681/1‑140 (MQ=255)

                         tgttataATAATAGAGTCGTCTATCTCTCAGGCGTCAA‑‑TGGACGAATAGAGGAGGTCCATTAC‑T‑T‑TGCTAATC‑CTTTTCGTTCGCATCATAGCACCAGTCAGCAGTGGCTGTGCTGTGGCGTATTATACTTATTGGCTTAg                                                                                                               <  2:270883/138‑1 (MQ=255)

                         tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                               <  1:612933/138‑4 (MQ=255)

                         tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                               <  2:415240/138‑4 (MQ=255)

                         gagtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                               <  1:638521/138‑4 (MQ=255)

                          ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCttggt                                                                                                              >  1:573848/3‑137 (MQ=255)

                          gatATAATAGTAGAGCCGCCTATCTCTCAGGCGACAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTAGTTCACACCATAGAACCAGGCAACAGTGGCTGTGGCATTGCGTGTTTTTcatcttcactaat                                                                                                               >  2:625437/1‑129 (MQ=255)

                           gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                             >  1:496047/2‑141 (MQ=255)

                           gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCTGTTGCGTATTATACTGATTGGCTTAGTa                                                                                                             >  1:74474/2‑141 (MQ=255)

                           atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                              >  1:352569/1‑141 (MQ=255)

                           atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                              >  1:652215/1‑141 (MQ=255)

                             ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                           >  2:493284/1‑139 (MQ=255)

                             ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                           <  1:445057/141‑3 (MQ=255)

                              aaataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                          <  1:99481/140‑4 (MQ=255)

                               aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaaag                                                                                                         >  2:647690/1‑137 (MQ=255)

                               aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                         >  2:111394/1‑141 (MQ=255)

                                ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTGAGTAGACg                                                                                                         >  2:122135/1‑141 (MQ=255)

                                ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                         <  2:290132/141‑1 (MQ=255)

                                 taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                       >  1:566730/4‑141 (MQ=255)

                                 taatagAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                       <  2:376500/138‑1 (MQ=255)

                                 gaatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGGACATTAC‑T‑T‑TGCTAATC‑ATTTGCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTGGTAAACgga                                                                                                       >  2:96135/4‑139 (MQ=255)

                                  agtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa                                                                                                       <  1:206125/141‑1 (MQ=255)

                                  aatagGGTCGCCTATATCTCAGGCGTCAATT‑‑GACGAACAGAGGCGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAATGCaa                                                                                                      <  1:469100/139‑1 (MQ=255)

                                  aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                      <  1:613537/139‑1 (MQ=255)

                                   atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                     >  2:498300/2‑140 (MQ=255)

                                    tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                     <  1:407826/141‑1 (MQ=255)

                                    tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGCGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                    <  1:292403/141‑3 (MQ=255)

                                    tagAGTCGCCTATCCCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                     <  1:53886/141‑1 (MQ=255)

                                     agctgctccTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAAGC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                   <  1:77642/134‑4 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                   >  1:270264/1‑138 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                   >  2:571661/1‑138 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                   >  2:297534/1‑138 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                   <  1:277141/141‑4 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TTGCGAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                   >  2:629138/1‑138 (MQ=255)

                                        gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                                                                 <  2:400338/141‑1 (MQ=255)

                                        gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                <  1:14081/141‑7 (MQ=255)

                                         tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                               <  1:618995/141‑8 (MQ=255)

                                         tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                               >  1:162660/1‑134 (MQ=255)

                                         tCGCCTATCTCCCAGGCCTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                               <  1:8777/141‑8 (MQ=255)

                                          cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                              >  1:261191/1‑133 (MQ=255)

                                          cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                              >  2:379385/1‑133 (MQ=255)

                                          cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                              >  1:358919/1‑133 (MQ=255)

                                           gCCTATCTCTCAGGCGTCAATT‑‑GACGCGATGAGGGTGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaacaatt                                                                                             <  2:300611/141‑10 (MQ=255)

                                           cccTATCTCTCGGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCAAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                             <  2:59850/140‑10 (MQ=255)

                                            ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                            >  2:173791/1‑131 (MQ=255)

                                            ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                            <  1:97113/141‑11 (MQ=255)

                                            ccTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                            <  1:602673/141‑11 (MQ=255)

                                              tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatttag                                                                                          >  2:164795/1‑129 (MQ=255)

                                              tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTATTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                          >  2:13627/1‑129 (MQ=255)

                                               ttctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAT‑T‑T‑TGCTAATC‑ATTTTCGTTCACGTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                         <  1:450432/140‑14 (MQ=255)

                                               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                         >  1:168804/1‑128 (MQ=255)

                                               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                         >  2:580444/1‑128 (MQ=255)

                                               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                         >  2:577827/1‑128 (MQ=255)

                                                        gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAg                                                                                 >  2:526166/1‑141 (MQ=255)

                                                             aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                                                                            <  2:454536/141‑1 (MQ=255)

                                                               tttGGACGCAGACAGGAGGTGTA‑TCA‑G‑G‑TGACGCTT‑ATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                                                                          <  2:63656/141‑1 (MQ=255)

                                                                             ggaggTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                            >  2:26504/1‑141 (MQ=255)

                                                                                      ttaA‑G‑G‑TGATGCGT‑CTTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATCGGCTAAGTAGACGCAATACAAAATCGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg                                                    <  1:86504/140‑1 (MQ=255)

                                                                                                 aTGCTT‑ATTTTCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGt                                            >  2:190265/1‑141 (MQ=255)

                                                                                                        aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGcccctcgagagggc                                      >  2:624889/1‑128 (MQ=255)

                                                                                                            tctgtcatatcataGCACCAGTCATCATTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTg                                  <  2:663060/137‑1 (MQ=255)

                                                                                                               ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                               <  1:338440/141‑1 (MQ=255)

                                                                                                                                gTCACCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAACATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt              <  1:133167/141‑1 (MQ=255)

                                                                                                                                       gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGtt       <  2:652215/141‑1 (MQ=255)

                                                                                                                                           cTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCt   >  1:226948/1‑141 (MQ=255)

                                                                                                                                            tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCGGTCCCGGCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGGCtt  >  1:153638/1‑141 (MQ=255)

                                                                                                                                            tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCtt  >  2:615313/1‑141 (MQ=255)

GTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑ATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTT  >  CP000730/1929777‑1930052

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