Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 630,742 C→A 17.3% E31* (GAA→TAA)  ybdM ← Spo0J family protein, ParB‑like nuclease domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273630,7420CA17.3% 882.0 / 111.9 280E31* (GAA→TAA) ybdMSpo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base C (115/116);  new base A (20/29);  total (135/145)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

AATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGG  >  CP009273/630593‑630891
                                                                                                                                                     |                                                                                                                                                     
aaTGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTc                                                                                                                                                       <  1:3112666/150‑1 (MQ=255)
aaTGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTc                                                                                                                                                       <  1:1655056/150‑1 (MQ=255)
     aTTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATttaattaa                                                                                                                                                  >  1:1361384/1‑150 (MQ=255)
      tttttGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAAt                                                                                                                                                 >  1:63171/1‑150 (MQ=255)
        tttGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCg                                                                                                                                               >  1:3801519/1‑150 (MQ=255)
              gcagTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGa                                                                                                                                         >  2:1562598/1‑150 (MQ=255)
                agTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATg                                                                                                                                       >  2:2782966/1‑150 (MQ=255)
                 gTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGc                                                                                                                                      <  2:3672032/150‑1 (MQ=255)
                 gTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGc                                                                                                                                      >  1:3419741/1‑150 (MQ=255)
                  ttttttttCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCg                                                                                                                                     <  2:4183458/150‑1 (MQ=255)
                   tttttttCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGa                                                                                                                                    >  1:3677478/1‑150 (MQ=255)
                   tttttttCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGa                                                                                                                                    >  1:1522989/1‑150 (MQ=255)
                        ttCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatcatc                                                                                                                               >  1:2674986/1‑150 (MQ=255)
                          cGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCtt                                                                                                                             >  1:2460200/1‑150 (MQ=255)
                          cGGGTGGTGCCACGTTGTTTGGGGTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGNTTTAATTAATCGCTTTGATGCGATCATCtt                                                                                                                             <  1:4285608/150‑1 (MQ=255)
                            ggtggtGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTcc                                                                                                                           <  2:2008914/150‑1 (MQ=255)
                              tggtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCgg                                                                                                                         >  2:3025372/1‑150 (MQ=255)
                               ggtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGc                                                                                                                        >  1:412053/1‑150 (MQ=255)
                                gtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCa                                                                                                                       <  2:4529829/150‑1 (MQ=255)
                                   ccACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGc                                                                                                                    <  1:1302630/150‑1 (MQ=255)
                                    cACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCt                                                                                                                   >  2:3934158/1‑150 (MQ=255)
                                    cACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCt                                                                                                                   >  2:1123071/1‑150 (MQ=255)
                                      cGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTgg                                                                                                                 >  2:2816806/1‑150 (MQ=255)
                                       gttgttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGc                                                                                                                >  1:439013/1‑150 (MQ=255)
                                       gttgttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGc                                                                                                                <  1:1270523/150‑1 (MQ=255)
                                        ttgttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCg                                                                                                               <  1:683178/150‑1 (MQ=255)
                                          gttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgag                                                                                                             >  1:4081971/1‑150 (MQ=255)
                                           tttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCgaga                                                                                                            <  1:1538383/150‑1 (MQ=255)
                                                 ttGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGc                                                                                                      >  1:2374154/1‑150 (MQ=255)
                                                     aaTCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                                                                                                  <  1:673128/150‑1 (MQ=255)
                                                     aaTCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                                                                                                  <  2:149884/150‑1 (MQ=255)
                                                        cATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                                                                                               <  1:2360860/150‑1 (MQ=255)
                                                        cATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                                                                                               <  1:1889217/150‑1 (MQ=255)
                                                         aTTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATc                                                                                              >  1:4540080/1‑150 (MQ=255)
                                                          ttAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                                                                                             <  2:2764128/150‑1 (MQ=255)
                                                          ttAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                                                                                             <  1:1225696/150‑1 (MQ=255)
                                                          ttAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                                                                                             >  1:3257672/1‑150 (MQ=255)
                                                            aGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg                                                                                           >  2:3084992/1‑150 (MQ=255)
                                                             ggCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGc                                                                                          >  2:263317/1‑150 (MQ=255)
                                                               cATAAGTTGGCTGATTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCATTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGt                                                                                        <  2:3371378/150‑1 (MQ=255)
                                                                aTAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGtt                                                                                       <  1:2024625/150‑1 (MQ=255)
                                                                 tAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTa                                                                                      <  2:377731/150‑1 (MQ=255)
                                                                  aaGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTaa                                                                                     <  2:861241/150‑1 (MQ=255)
                                                                  aaGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTaa                                                                                     >  1:443453/1‑150 (MQ=255)
                                                                   aGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAAt                                                                                    >  1:2540964/1‑150 (MQ=255)
                                                                   aGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAAt                                                                                    >  1:1790885/1‑150 (MQ=255)
                                                                    gTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATc                                                                                   <  2:4364998/150‑1 (MQ=255)
                                                                    gTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATc                                                                                   >  1:1430426/1‑150 (MQ=255)
                                                                       ggCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgtt                                                                                <  1:2716713/150‑1 (MQ=255)
                                                                        gCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttg                                                                               <  2:2320894/150‑1 (MQ=255)
                                                                          tGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt                                                                             <  1:3020189/150‑1 (MQ=255)
                                                                               ttAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATa                                                                        <  1:631353/150‑1 (MQ=255)
                                                                                 aaCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGt                                                                      <  2:554964/150‑1 (MQ=255)
                                                                                 aaCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGt                                                                      >  2:1829019/1‑150 (MQ=255)
                                                                                 aaCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGATAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGTCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCTTTGTTGCATAGt                                                                      <  2:517760/150‑1 (MQ=255)
                                                                                  aCCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTa                                                                     >  1:609262/1‑150 (MQ=255)
                                                                                   cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTAt                                                                    >  1:3370575/1‑150 (MQ=255)
                                                                                   cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTAt                                                                    >  2:3571008/1‑150 (MQ=255)
                                                                                    ccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATc                                                                   >  2:2705280/1‑150 (MQ=255)
                                                                                      acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATccc                                                                 >  2:2301415/1‑150 (MQ=255)
                                                                                        aGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATccccc                                                               <  2:599335/150‑1 (MQ=255)
                                                                                        aGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATccccc                                                               >  1:4212942/1‑150 (MQ=255)
                                                                                         gAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCa                                                              <  1:2588589/150‑1 (MQ=255)
                                                                                          aaCGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCAt                                                             <  2:82378/150‑1 (MQ=255)
                                                                                            cGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTc                                                           >  1:72317/1‑150 (MQ=255)
                                                                                              cAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTcat                                                         <  2:1523036/150‑3 (MQ=255)
                                                                                                aTCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGa                                                       <  1:196179/150‑1 (MQ=255)
                                                                                                  ccACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACg                                                     >  2:3284092/1‑150 (MQ=255)
                                                                                                  ccACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACg                                                     >  2:2487292/1‑150 (MQ=255)
                                                                                                   cACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGt                                                    >  1:1203451/1‑150 (MQ=255)
                                                                                                   cACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGt                                                    >  1:3816890/1‑150 (MQ=255)
                                                                                                    aCTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTGTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGtt                                                   >  1:491653/1‑150 (MQ=255)
                                                                                                     cTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGttt                                                  <  2:2369579/150‑1 (MQ=255)
                                                                                                      tGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGtttt                                                 >  1:1069180/1‑150 (MQ=255)
                                                                                                             tCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTa                                          >  1:2122856/1‑150 (MQ=255)
                                                                                                               aCGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATc                                        <  1:2285718/150‑1 (MQ=255)
                                                                                                                cGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATcc                                       <  1:903205/150‑1 (MQ=255)
                                                                                                                cGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATcc                                       <  2:3801641/150‑1 (MQ=255)
                                                                                                                    aGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgttc                                   <  1:1021350/150‑1 (MQ=255)
                                                                                                                     ggggCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgttcg                                  <  2:2674986/150‑1 (MQ=255)
                                                                                                                      gggCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgttcgt                                 <  1:875879/150‑1 (MQ=255)
                                                                                                                           aaCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATaa                            >  1:1839050/1‑150 (MQ=255)
                                                                                                                           aaCCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATaa                            >  2:2949868/1‑150 (MQ=255)
                                                                                                                                gATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTg                       >  1:1729746/1‑150 (MQ=255)
                                                                                                                                   ggATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATa                    >  1:3995274/1‑150 (MQ=255)
                                                                                                                                   ggATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATa                    >  1:4488689/1‑150 (MQ=255)
                                                                                                                                     aTCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGt                  >  2:4472811/1‑150 (MQ=255)
                                                                                                                                      tCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGtt                 >  2:1006947/1‑150 (MQ=255)
                                                                                                                                      tCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGtt                 <  2:494159/150‑1 (MQ=255)
                                                                                                                                      tCGCCATGCGGATTTAATTAATCGCTTTGATGCTATCATCTTCCGGCAAGCTGGCGAGAAATTGAGTTAAATCCTGCTTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAgat                 <  2:2088894/150‑3 (MQ=255)
                                                                                                                                        gCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCGTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGttct               <  2:4096690/150‑1 (MQ=255)
                                                                                                                                         ccATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGttctt              <  2:1328740/150‑1 (MQ=255)
                                                                                                                                          cATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGttcttc             <  2:1384551/150‑1 (MQ=255)
                                                                                                                                          cATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGttcttc             >  1:4278114/1‑150 (MQ=255)
                                                                                                                                           aTGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCg            >  1:3215542/1‑150 (MQ=255)
                                                                                                                                           aTGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCg            >  2:1831727/1‑150 (MQ=255)
                                                                                                                                              cGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAtt         <  1:3758065/150‑1 (MQ=255)
                                                                                                                                               ggATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAttt        >  1:1822253/1‑150 (MQ=255)
                                                                                                                                               ggATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAttt        >  2:3872827/1‑150 (MQ=255)
                                                                                                                                                gATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAtttt       <  1:1387134/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGG  >  CP009273/630593‑630891

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: