Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,137,659 C→A 15.5% R305R (CGG→AGG)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,137,6590CA15.5% 1041.0 / 113.8 320R305R (CGG→AGG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (140/130);  new base A (23/27);  total (163/157)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 321 total aligned reads displayed.

CTGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAG  >  CP009273/2137512‑2137808
                                                                                                                                                   |                                                                                                                                                     
ctGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAgg                                                                                                                                                     <  1:715005/150‑1 (MQ=255)
 tgcttgctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGa                                                                                                                                                    <  1:643088/150‑1 (MQ=255)
 tgcttgctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGa                                                                                                                                                    <  1:601430/150‑1 (MQ=255)
    ttgctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAAcg                                                                                                                                                 <  1:3587715/150‑1 (MQ=255)
      gctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAAcgcg                                                                                                                                               <  1:3249545/150‑1 (MQ=255)
       ctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGa                                                                                                                                              >  1:2807907/1‑150 (MQ=255)
        tACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGAt                                                                                                                                             <  2:2340467/150‑1 (MQ=255)
         aCGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATg                                                                                                                                            <  1:2720440/150‑1 (MQ=255)
          cgccgcGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATgg                                                                                                                                           >  2:872531/1‑150 (MQ=255)
           gccgcGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGa                                                                                                                                          >  2:4413693/1‑150 (MQ=255)
            ccgcGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGaa                                                                                                                                         <  2:2029155/150‑1 (MQ=255)
            ccgcGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGaa                                                                                                                                         <  1:2398305/150‑1 (MQ=255)
                gtacgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATAtt                                                                                                                                     >  1:935820/1‑150 (MQ=255)
                 tacgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTc                                                                                                                                    <  2:2526206/150‑1 (MQ=255)
                 tacgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTc                                                                                                                                    >  1:3389313/1‑150 (MQ=255)
                   cgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCg                                                                                                                                  <  1:1255945/150‑1 (MQ=255)
                   cgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCg                                                                                                                                  <  1:2136193/150‑1 (MQ=255)
                    gtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGc                                                                                                                                 >  2:1172019/1‑150 (MQ=255)
                     tacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCt                                                                                                                                <  2:4306290/150‑1 (MQ=255)
                     tacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCt                                                                                                                                >  1:176636/1‑150 (MQ=255)
                           atgaGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGt                                                                                                                          <  1:1380609/150‑1 (MQ=255)
                            tgaGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCACCATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTa                                                                                                                         <  2:2475937/150‑1 (MQ=255)
                              aGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATg                                                                                                                       <  1:1674138/150‑1 (MQ=255)
                                ccATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGACTATTCCGCTATCGGTATGGc                                                                                                                     >  2:4622858/1‑150 (MQ=255)
                                   tatGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGtt                                                                                                                  <  1:2621671/150‑1 (MQ=255)
                                   tatGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGCATGGCGtt                                                                                                                  >  2:223985/1‑150 (MQ=255)
                                     tGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAg                                                                                                                >  1:419703/1‑150 (MQ=255)
                                             cTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCAcc                                                                                                        <  1:4000526/150‑1 (MQ=255)
                                             cTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCAcc                                                                                                        >  1:1848734/1‑150 (MQ=255)
                                             cTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCAcc                                                                                                        >  1:3864312/1‑150 (MQ=255)
                                             cTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGACTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCAcc                                                                                                        <  1:4641330/150‑1 (MQ=255)
                                                 cGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAgg                                                                                                    >  2:3757927/1‑150 (MQ=255)
                                                 cGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAgg                                                                                                    <  1:2224886/150‑1 (MQ=255)
                                                 cGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAgg                                                                                                    >  1:4612029/1‑150 (MQ=255)
                                                  gTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAggg                                                                                                   <  2:3383294/150‑1 (MQ=255)
                                                    ttttGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGAc                                                                                                 <  1:3666427/150‑1 (MQ=255)
                                                      ttGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACaa                                                                                               >  1:2289347/1‑150 (MQ=255)
                                                         cTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATgg                                                                                            >  1:1695907/1‑150 (MQ=255)
                                                           gATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCt                                                                                          >  1:3812352/1‑150 (MQ=255)
                                                             tCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGc                                                                                        >  1:4399154/1‑150 (MQ=255)
                                                             tCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGc                                                                                        <  1:1335945/150‑1 (MQ=255)
                                                             tCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGc                                                                                        <  2:2075127/150‑1 (MQ=255)
                                                             tCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTGTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGc                                                                                        >  1:1491756/1‑150 (MQ=255)
                                                              ccTGCTGCCCGCCAACATCATGACCATGGTGATTTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCa                                                                                       <  2:4087983/150‑1 (MQ=255)
                                                                 gctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAAc                                                                                    <  2:2349157/150‑1 (MQ=255)
                                                                 gctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAAc                                                                                    >  1:1157098/1‑150 (MQ=255)
                                                                  ctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTAGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAAcc                                                                                   <  1:1468243/150‑1 (MQ=255)
                                                                     cccGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACcaac                                                                                >  2:3254569/1‑150 (MQ=255)
                                                                         ccAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAg                                                                            <  2:2611321/150‑1 (MQ=255)
                                                                         ccAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAg                                                                            <  1:3989721/150‑1 (MQ=255)
                                                                          cAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAgc                                                                           >  2:675135/1‑150 (MQ=255)
                                                                              atcatGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACAGAGGGACAATGGCTGCAAACCAACAAAGCGctc                                                                       >  2:476247/1‑150 (MQ=255)
                                                                               tcatGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGctct                                                                      <  1:3622504/150‑1 (MQ=255)
                                                                               tcatGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGctct                                                                      <  2:1353918/150‑1 (MQ=255)
                                                                               tcatGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGctct                                                                      <  1:3517755/150‑1 (MQ=255)
                                                                                  tGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGcc                                                                   >  1:2571923/1‑150 (MQ=255)
                                                                                  tGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGcc                                                                   >  1:1603381/1‑150 (MQ=255)
                                                                                  tGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGcc                                                                   >  2:1920497/1‑150 (MQ=255)
                                                                                   gACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCa                                                                  <  1:3654299/150‑1 (MQ=255)
                                                                                   gACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCa                                                                  >  2:3208809/1‑150 (MQ=255)
                                                                                     ccATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGt                                                                <  1:3152260/150‑1 (MQ=255)
                                                                                       aTGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGttt                                                              >  2:4220663/1‑150 (MQ=255)
                                                                                       aTGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGttt                                                              >  1:4460861/1‑150 (MQ=255)
                                                                                        tGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtc                                                             >  2:1769627/1‑150 (MQ=255)
                                                                                        tGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtc                                                             <  1:3482577/150‑1 (MQ=255)
                                                                                         ggTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtct                                                            >  1:309069/1‑150 (MQ=255)
                                                                                         ggTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtct                                                            <  2:3928648/150‑1 (MQ=255)
                                                                                          gTGATGTATGTCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtctc                                                           <  1:456882/150‑1 (MQ=255)
                                                                                          gTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtctc                                                           <  1:112618/150‑1 (MQ=255)
                                                                                             aTGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCggg                                                        <  2:935820/150‑1 (MQ=255)
                                                                                                    ccTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTc                                                 <  2:3335486/150‑1 (MQ=255)
                                                                                                     cTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCa                                                <  1:93430/150‑1 (MQ=255)
                                                                                                          gTGCGGAACGCAAACGCATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAag                                           <  2:3680135/150‑1 (MQ=255)
                                                                                                             cGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGc                                        >  2:310411/1‑150 (MQ=255)
                                                                                                               gAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTg                                      >  2:2920500/1‑150 (MQ=255)
                                                                                                                aaCGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTgc                                     >  2:777365/1‑150 (MQ=255)
                                                                                                                 aCGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTgcg                                    >  2:2161536/1‑150 (MQ=255)
                                                                                                                  cGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTgcgc                                   >  2:1436551/1‑150 (MQ=255)
                                                                                                                  cGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTgcgc                                   >  2:360916/1‑150 (MQ=255)
                                                                                                                   gCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGcgcg                                  <  2:3900816/150‑1 (MQ=255)
                                                                                                                   gCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGcgcg                                  <  2:901819/150‑1 (MQ=255)
                                                                                                                    cAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCgg                                 >  2:965281/1‑150 (MQ=255)
                                                                                                                    cAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCgg                                 <  1:1999781/150‑1 (MQ=255)
                                                                                                                     aaaCACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGa                                >  1:2315780/1‑150 (MQ=255)
                                                                                                                     aaaCACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGa                                >  2:1423424/1‑150 (MQ=255)
                                                                                                                        cacaTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTc                             <  2:1059183/150‑1 (MQ=255)
                                                                                                                          caTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCAc                           >  1:737368/1‑150 (MQ=255)
                                                                                                                               ccGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTc                      >  2:1293014/1‑150 (MQ=255)
                                                                                                                                 gAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAg                    >  2:3393221/1‑150 (MQ=255)
                                                                                                                                      cGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGTCTCACCTTTCAGCAACt               <  2:3463789/150‑1 (MQ=255)
                                                                                                                                       gAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTg              >  2:3821045/1‑150 (MQ=255)
                                                                                                                                       gAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTg              >  1:2435706/1‑150 (MQ=255)
                                                                                                                                        aaaCCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGa             >  2:2850459/1‑150 (MQ=255)
                                                                                                                                          aCCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGAcc           <  1:275347/150‑1 (MQ=255)
                                                                                                                                            ccACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTg         >  2:2972955/1‑150 (MQ=255)
                                                                                                                                            ccACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTg         <  2:1848734/150‑1 (MQ=255)
                                                                                                                                             cACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTgg        >  2:1682420/1‑150 (MQ=255)
                                                                                                                                              aCTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGc       >  1:4257546/1‑150 (MQ=255)
                                                                                                                                               cTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGcc      >  1:1000539/1‑150 (MQ=255)
                                                                                                                                               cTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGcc      <  2:2289347/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAG  >  CP009273/2137512‑2137808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: