Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,137,659 C→A 38.2% 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,6590CA38.2% 390.7 / 295.2 279R305R (CGG→AGG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (83/89);  new base A (49/58);  total (132/147)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

CCCTGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAG  >  CP009273/2137510‑2137808
                                                                                                                                                     |                                                                                                                                                     
ccctGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTc                                                                                                                                                       <  2:2176560/150‑1 (MQ=255)
  ctGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCgg                                                                                                                                                     <  2:3434181/150‑1 (MQ=255)
     cttgctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAAc                                                                                                                                                  <  1:3770220/150‑1 (MQ=255)
        gctACGCCGCGTACGTACCTGTTGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAAcgcg                                                                                                                                               <  1:3227557/150‑1 (MQ=255)
        gctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAAcgcg                                                                                                                                               <  2:1059663/150‑1 (MQ=255)
        gctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAAcgcg                                                                                                                                               <  2:2945382/150‑1 (MQ=255)
         ctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGa                                                                                                                                              <  2:4489321/150‑1 (MQ=255)
         ctACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGa                                                                                                                                              >  2:1107381/1‑150 (MQ=255)
            cgccgcGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATgg                                                                                                                                           <  2:4447834/150‑1 (MQ=255)
             gccgcGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGa                                                                                                                                          >  1:3503354/1‑150 (MQ=255)
               cgcgTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAat                                                                                                                                        <  1:2993406/150‑1 (MQ=255)
                gcgTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAata                                                                                                                                       >  1:979916/1‑150 (MQ=255)
                  gtacgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATAtt                                                                                                                                     >  1:280670/1‑150 (MQ=255)
                   tacgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTc                                                                                                                                    <  1:311358/150‑1 (MQ=255)
                     cgtacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCg                                                                                                                                  <  2:1417629/150‑1 (MQ=255)
                       tacCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCt                                                                                                                                >  2:1576241/1‑150 (MQ=255)
                         ccTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTAt                                                                                                                              <  1:4508257/150‑1 (MQ=255)
                          cTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATc                                                                                                                             <  2:4269658/150‑1 (MQ=255)
                           tgatgaGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCg                                                                                                                            <  1:799541/150‑1 (MQ=255)
                           tgatgaGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCg                                                                                                                            <  1:2782242/150‑1 (MQ=255)
                            gatgaGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCgg                                                                                                                           >  2:3487638/1‑150 (MQ=255)
                              tgaGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTa                                                                                                                         >  1:4318180/1‑150 (MQ=255)
                                aGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATg                                                                                                                       <  1:2622626/150‑1 (MQ=255)
                                aGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCGTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCGATCGGTATg                                                                                                                       >  2:1767421/1‑150 (MQ=255)
                                 gCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATgg                                                                                                                      >  2:587027/1‑150 (MQ=255)
                                  ccATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGc                                                                                                                     <  1:4574061/150‑1 (MQ=255)
                                    atatGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGt                                                                                                                   >  1:918760/1‑150 (MQ=255)
                                    atatGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGAATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGt                                                                                                                   >  2:2303505/1‑150 (MQ=255)
                                       tGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAg                                                                                                                >  1:1383292/1‑150 (MQ=255)
                                       tGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAg                                                                                                                <  2:494013/150‑1 (MQ=255)
                                        gCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGt                                                                                                               >  2:1683593/1‑150 (MQ=255)
                                           gTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTggg                                                                                                            >  2:2201894/1‑150 (MQ=255)
                                             ggCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCa                                                                                                          <  1:1219185/150‑1 (MQ=255)
                                             ggCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCa                                                                                                          >  1:3462333/1‑150 (MQ=255)
                                              gCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCAc                                                                                                         <  1:4192162/150‑1 (MQ=255)
                                               cTACCGATTTTGAAGATCCTGCTTCCCGCCAACATCATGACCATGGTGAGGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACACACTTTCGGACCGCGAGGGAATATACCGCTAACGGTATGGCGTTAGTGGGCAcc                                                                                                        >  2:1303450/1‑150 (MQ=255)
                                                 aCCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGa                                                                                                      >  2:2416258/1‑150 (MQ=255)
                                                    gTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAggg                                                                                                   >  1:1960271/1‑150 (MQ=255)
                                                      ttttGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGAc                                                                                                 >  2:4621446/1‑150 (MQ=255)
                                                                ccTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCa                                                                                       <  1:3926870/150‑1 (MQ=255)
                                                                  tgctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCaaa                                                                                     >  1:3303923/1‑150 (MQ=255)
                                                                   gctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAAc                                                                                    <  2:2320589/150‑1 (MQ=255)
                                                                    ctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAAcc                                                                                   >  1:1039173/1‑150 (MQ=255)
                                                                    ctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAAcc                                                                                   <  1:2812024/150‑1 (MQ=255)
                                                                    ctgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAAcc                                                                                   >  2:3398091/1‑150 (MQ=255)
                                                                     tgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACca                                                                                  <  2:2661301/150‑1 (MQ=255)
                                                                     tgcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACca                                                                                  >  1:2686795/1‑150 (MQ=255)
                                                                          gCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACaaa                                                                             >  2:2513969/1‑150 (MQ=255)
                                                                            cAACATCATGACCATGGTGATGTAGGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAgc                                                                           <  1:2144860/150‑1 (MQ=255)
                                                                             aaCATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAgcg                                                                          >  1:2171120/1‑150 (MQ=255)
                                                                              tcatcatGACCATGGTGATGTTTGCCTTTCGTTCGGTATGCAAACACATTTCCGAAAGCGTATCCCTCTTTAGGAACGCGATGGATTATTCCGCTATCGGTATGGAGTTAGTGGGCACCGAGGGACTATGGCTGCTAAGCAACAAAgcgc                                                                         <  1:630715/149‑1 (MQ=255)
                                                                                 tcatGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGctct                                                                      <  2:4450529/150‑1 (MQ=255)
                                                                                  catGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACACTGGCTGCAAACCAACAAAGCGCTCTg                                                                     >  2:3576429/1‑150 (MQ=255)
                                                                                   atGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGc                                                                    <  1:4386858/150‑1 (MQ=255)
                                                                                    tGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGcc                                                                   >  1:3578443/1‑150 (MQ=255)
                                                                                         aTGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGttt                                                              <  2:2045796/150‑1 (MQ=255)
                                                                                         aTGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGttt                                                              >  1:2504976/1‑150 (MQ=255)
                                                                                         aTGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGttt                                                              <  2:3017527/150‑1 (MQ=255)
                                                                                          tGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCCCTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtc                                                             <  2:3121706/150‑1 (MQ=255)
                                                                                          tGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtc                                                             <  2:1142439/150‑1 (MQ=255)
                                                                                          tGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtc                                                             <  1:649659/150‑1 (MQ=255)
                                                                                            gTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtctc                                                           >  2:1519664/1‑150 (MQ=255)
                                                                                             tGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCg                                                          <  2:1140073/150‑1 (MQ=255)
                                                                                              gATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCgg                                                         <  1:1086079/150‑1 (MQ=255)
                                                                                                  tATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTAc                                                     >  1:1723019/1‑150 (MQ=255)
                                                                                                    tGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAg                                                   >  1:1378058/1‑150 (MQ=255)
                                                                                                           cGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGaa                                            >  1:1924444/1‑150 (MQ=255)
                                                                                                            gTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAag                                           <  1:3371878/150‑1 (MQ=255)
                                                                                                              gCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAagag                                         >  1:4133013/1‑150 (MQ=255)
                                                                                                               cGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGc                                        >  2:3717320/1‑150 (MQ=255)
                                                                                                               cGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGc                                        <  1:2985617/150‑1 (MQ=255)
                                                                                                                ggAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCt                                       <  2:2558042/150‑1 (MQ=255)
                                                                                                                ggAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCt                                       >  2:4086456/1‑150 (MQ=255)
                                                                                                                 gAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTg                                      <  1:267568/150‑1 (MQ=255)
                                                                                                                  aaCGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTgc                                     >  2:1389049/1‑150 (MQ=255)
                                                                                                                   aCGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTgcg                                    <  2:4517056/150‑1 (MQ=255)
                                                                                                                      cAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCgg                                 >  1:1606443/1‑150 (MQ=255)
                                                                                                                         acacATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACt                              <  1:1545068/150‑1 (MQ=255)
                                                                                                                             aTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCAcc                          >  2:1242173/1‑150 (MQ=255)
                                                                                                                             aTTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCAcc                          >  2:2720038/1‑150 (MQ=255)
                                                                                                                                tCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCttt                       >  1:874395/1‑150 (MQ=255)
                                                                                                                                 ccGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTc                      <  1:1127897/150‑1 (MQ=255)
                                                                                                                                   gAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAg                    <  1:500523/150‑1 (MQ=255)
                                                                                                                                    aaaGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGc                   <  1:517933/150‑1 (MQ=255)
                                                                                                                                       gCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAAc                <  1:4114755/150‑1 (MQ=255)
                                                                                                                                        cGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACt               <  1:1490232/150‑1 (MQ=255)
                                                                                                                                        cGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACt               <  1:176858/150‑1 (MQ=255)
                                                                                                                                        cGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACt               <  1:2207947/150‑1 (MQ=255)
                                                                                                                                        cGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACt               <  1:1672101/150‑1 (MQ=255)
                                                                                                                                        cGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACt               <  1:1700306/150‑1 (MQ=255)
                                                                                                                                         gAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTg              >  2:2484583/1‑150 (MQ=255)
                                                                                                                                         gAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTg              >  1:4547869/1‑150 (MQ=255)
                                                                                                                                          aaaCCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGa             <  2:2868742/150‑1 (MQ=255)
                                                                                                                                          aaaCCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGa             >  2:4056754/1‑150 (MQ=255)
                                                                                                                                           aaCCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGAc            >  1:571943/1‑150 (MQ=255)
                                                                                                                                            aCCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGAcc           <  1:4246750/150‑1 (MQ=255)
                                                                                                                                             cccACTTTCTGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCt          <  1:746780/150‑1 (MQ=255)
                                                                                                                                              ccACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTg         >  1:2027433/1‑150 (MQ=255)
                                                                                                                                                aCTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGc       >  1:3839944/1‑150 (MQ=255)
                                                                                                                                                aCTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGc       >  1:3755188/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCTGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAG  >  CP009273/2137510‑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: