Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,068,854 C→A 42.5% E230* (GAG→TAG)  yggF ← fructose 1,6 bisphosphatase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,068,8540CA42.5% 332.7 / inf 324E230* (GAG→TAG) yggFfructose 1,6 bisphosphatase isozyme
Reads supporting (aligned to +/- strand):  ref base C (111/75);  new base A (72/66);  total (183/141)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

AAGAGGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAA  >  CP009273/3068706‑3069003
                                                                                                                                                    |                                                                                                                                                     
     gATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTg                                                                                                                                                 <  2:485853/150‑1 (MQ=255)
       tatCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCa                                                                                                                                               >  1:1070183/1‑150 (MQ=255)
        atCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCat                                                                                                                                              >  2:1150603/1‑150 (MQ=255)
          cGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCatat                                                                                                                                            <  1:73634/150‑1 (MQ=255)
           gTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATc                                                                                                                                           >  1:2714262/1‑150 (MQ=255)
            ttCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATcc                                                                                                                                          >  2:976561/1‑150 (MQ=255)
              cccccTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATcccc                                                                                                                                        <  1:3300620/150‑1 (MQ=255)
               ccccTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCg                                                                                                                                       >  2:2380072/1‑150 (MQ=255)
                cccTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGc                                                                                                                                      >  1:369339/1‑150 (MQ=255)
                   tCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTa                                                                                                                                   >  1:48430/1‑150 (MQ=255)
                    cACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTaa                                                                                                                                  >  1:2338555/1‑150 (MQ=255)
                       cAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGc                                                                                                                               <  1:1964249/150‑1 (MQ=255)
                       cAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGc                                                                                                                               >  1:2387312/1‑150 (MQ=255)
                        aGTTCATCGAGCGAGTACTCACGGTTAACGTCAACGCCCATTGCCTTGCTACGCTTGCGCTCCTGCTCTGCAATCTGACGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGATCGGCCTGCATTTCACCGCCTAATGCt                                                                                                                              <  2:354623/150‑1 (MQ=255)
                          ttCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCttt                                                                                                                            <  2:602318/150‑1 (MQ=255)
                           tCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTNTa                                                                                                                           <  1:2111524/150‑1 (MQ=255)
                            cATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTaa                                                                                                                          >  2:4059192/1‑150 (MQ=255)
                            cATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTaa                                                                                                                          >  2:3651076/1‑150 (MQ=255)
                             aTCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                                                                                         >  1:3821822/1‑150 (MQ=255)
                               cgagcgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACgg                                                                                                                       >  2:3986334/1‑150 (MQ=255)
                                 agcgagTACACCCGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGca                                                                                                                     <  2:2147121/150‑1 (MQ=255)
                                   cgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGcaca                                                                                                                   <  2:2703840/150‑1 (MQ=255)
                                    gagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                                                                                  >  1:2861479/1‑150 (MQ=255)
                                    gagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                                                                                  >  2:1875788/1‑150 (MQ=255)
                                    gagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                                                                                  <  2:3745530/150‑1 (MQ=255)
                                     agTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAgg                                                                                                                 >  1:3119090/1‑150 (MQ=255)
                                     agTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAgg                                                                                                                 <  1:1292708/150‑1 (MQ=255)
                                      gTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGc                                                                                                                >  1:3184988/1‑150 (MQ=255)
                                      gTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGc                                                                                                                >  1:1196450/1‑150 (MQ=255)
                                       tACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCa                                                                                                               <  1:536669/150‑1 (MQ=255)
                                        acacacGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAg                                                                                                              <  2:526658/150‑1 (MQ=255)
                                          acacGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                                                                                            <  2:1739547/150‑1 (MQ=255)
                                          acacGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                                                                                            >  2:759833/1‑150 (MQ=255)
                                            acGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACAGCACAGGCAGAAAt                                                                                                          <  2:3421295/150‑1 (MQ=255)
                                               gTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCAc                                                                                                       <  2:28192/150‑1 (MQ=255)
                                                 tAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGc                                                                                                     >  1:3729459/1‑150 (MQ=255)
                                                  aaCGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGcc                                                                                                    >  1:915586/1‑150 (MQ=255)
                                                    cGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCtt                                                                                                  >  2:1476915/1‑150 (MQ=255)
                                                     gTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTc                                                                                                 >  1:1485054/1‑150 (MQ=255)
                                                       cAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTg                                                                                               >  2:258637/1‑150 (MQ=255)
                                                       cAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTg                                                                                               >  2:321149/1‑150 (MQ=255)
                                                            cccATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCa                                                                                          >  1:3056315/1‑150 (MQ=255)
                                                             ccATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAc                                                                                         >  2:3274369/1‑150 (MQ=255)
                                                                 tGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgcc                                                                                     <  1:1336933/150‑1 (MQ=255)
                                                                  gCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgccg                                                                                    >  2:91745/1‑150 (MQ=255)
                                                                    cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGAt                                                                                  <  2:266940/150‑1 (MQ=255)
                                                                    cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGAt                                                                                  <  2:3826416/150‑1 (MQ=255)
                                                                      tGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCCTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATgg                                                                                >  1:3728440/1‑150 (MQ=255)
                                                                       gCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgt                                                                               >  2:2562327/1‑150 (MQ=255)
                                                                       gCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgt                                                                               <  1:3992762/150‑1 (MQ=255)
                                                                        cAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtg                                                                              >  1:1995032/1‑150 (MQ=255)
                                                                        cAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtg                                                                              >  1:4067356/1‑150 (MQ=255)
                                                                        cAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtg                                                                              <  1:3000270/150‑1 (MQ=255)
                                                                         aaCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtgt                                                                             >  1:2681211/1‑150 (MQ=255)
                                                                         aaCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtgt                                                                             >  1:1382827/1‑150 (MQ=255)
                                                                          aCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTa                                                                            >  2:691942/1‑150 (MQ=255)
                                                                              ttGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatc                                                                        <  1:823238/150‑1 (MQ=255)
                                                                                 cgcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcaca                                                                     >  2:1686786/1‑150 (MQ=255)
                                                                                  gcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcacat                                                                    <  1:3351064/150‑1 (MQ=255)
                                                                                  gcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcacat                                                                    >  2:2468513/1‑150 (MQ=255)
                                                                                    tCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCg                                                                  >  2:3141224/1‑150 (MQ=255)
                                                                                        gCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAg                                                              <  1:688435/150‑1 (MQ=255)
                                                                                         ctctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATGTCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAgg                                                             <  1:3083522/150‑1 (MQ=255)
                                                                                         ctctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAgg                                                             <  1:1035113/150‑1 (MQ=255)
                                                                                          tctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGa                                                            <  1:2742986/150‑1 (MQ=255)
                                                                                              cAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTAt                                                        >  1:2427469/1‑150 (MQ=255)
                                                                                                aTCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATcc                                                      <  2:1096252/150‑1 (MQ=255)
                                                                                                  cTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTg                                                    <  2:1424663/150‑1 (MQ=255)
                                                                                                          tttCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAgg                                            >  1:2177980/1‑150 (MQ=255)
                                                                                                          tttCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAgg                                            >  2:630211/1‑150 (MQ=255)
                                                                                                            tCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTc                                          >  1:2081627/1‑150 (MQ=255)
                                                                                                              cGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAg                                        >  1:1538136/1‑150 (MQ=255)
                                                                                                                  tAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACg                                    <  2:1638911/150‑1 (MQ=255)
                                                                                                                   aaTCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGc                                   >  2:877212/1‑150 (MQ=255)
                                                                                                                    aTCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCt                                  >  2:3854806/1‑150 (MQ=255)
                                                                                                                     tCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTg                                 <  2:536220/150‑1 (MQ=255)
                                                                                                                     tCCCCTTTAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTg                                 <  2:2618763/150‑1 (MQ=255)
                                                                                                                           ttAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGgcggcg                           <  1:426860/150‑1 (MQ=255)
                                                                                                                           ttAGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGgcggcg                           >  2:1009462/1‑150 (MQ=255)
                                                                                                                            tAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGa                          >  2:2989698/1‑150 (MQ=255)
                                                                                                                            tAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGa                          <  2:818784/150‑1 (MQ=255)
                                                                                                                             aGCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGAc                         >  1:365285/1‑150 (MQ=255)
                                                                                                                              gCCTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACa                        >  1:1786678/1‑150 (MQ=255)
                                                                                                                               ccTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACAt                       >  2:3300303/1‑150 (MQ=255)
                                                                                                                                cTGGCAAAAGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATc                      <  1:3451421/150‑1 (MQ=255)
                                                                                                                                  ggCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGc                    >  1:167767/1‑150 (MQ=255)
                                                                                                                                      aaaGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATc                <  1:2238359/150‑1 (MQ=255)
                                                                                                                                      aaaGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATc                >  2:2941331/1‑150 (MQ=255)
                                                                                                                                       aaGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATcc               <  2:1241050/150‑1 (MQ=255)
                                                                                                                                       aaGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATcc               <  2:2869205/150‑1 (MQ=255)
                                                                                                                                        aGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCg              <  2:2718746/150‑1 (MQ=255)
                                                                                                                                        aGTCAATCAGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCg              >  1:1600830/1‑150 (MQ=255)
                                                                                                                                         gTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCgg             >  2:210764/1‑150 (MQ=255)
                                                                                                                                          tCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGc            <  2:4158904/150‑1 (MQ=255)
                                                                                                                                           cAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCa           <  2:733159/150‑1 (MQ=255)
                                                                                                                                           cAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCa           >  1:753731/1‑150 (MQ=255)
                                                                                                                                              tCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTACGGATTATCCTGCCAGCCGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAggg        >  1:2927502/1‑150 (MQ=255)
                                                                                                                                                aGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCa      >  1:2576659/1‑150 (MQ=255)
                                                                                                                                                aGCTAGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCa      >  2:4004691/1‑150 (MQ=255)
                                                                                                                                                  cTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCaaa    >  2:1287327/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGAGGATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAA  >  CP009273/3068706‑3069003

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: