Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 109,050 C→T 5.0% E13K (GAA→AAA)  zapD ← FtsZ stabilizer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273109,0500CT 5.0% 1099.3 / 28.8 297E13K (GAA→AAA) zapDFtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base C (145/137);  new base T (4/11);  total (149/148)
Fisher's exact test for biased strand distribution p-value = 6.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 297 total aligned reads displayed.

TCAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAT  >  CP009273/108902‑109198
                                                                                                                                                    |                                                                                                                                                    
 cAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCAt                                                                                                                                                    <  2:2460629/150‑1 (MQ=255)
 cAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCAt                                                                                                                                                    >  2:1708383/1‑150 (MQ=255)
     gCGGACTTCGCCGCGCTCGAAAACATCCAGTAAGTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAg                                                                                                                                                <  1:498623/150‑1 (MQ=255)
      cGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGt                                                                                                                                               >  2:1027470/1‑150 (MQ=255)
       ggACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTg                                                                                                                                              >  2:3501148/1‑150 (MQ=255)
       ggACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTg                                                                                                                                              >  2:4012861/1‑150 (MQ=255)
         aCTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGa                                                                                                                                            <  1:2079233/150‑1 (MQ=255)
          cgtCGCCGCGCACGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGAt                                                                                                                                           <  2:4262508/148‑1 (MQ=255)
             cgccgcGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGtt                                                                                                                                        >  2:2197216/1‑150 (MQ=255)
             cgccgcGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGtt                                                                                                                                        >  2:4310495/1‑150 (MQ=255)
              gccgcGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTc                                                                                                                                       >  2:3145844/1‑150 (MQ=255)
                  cgcTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCaaaa                                                                                                                                   <  2:3840345/150‑1 (MQ=255)
                    cTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAg                                                                                                                                 <  2:4152220/150‑1 (MQ=255)
                       gAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGAc                                                                                                                              <  1:906055/150‑1 (MQ=255)
                       gAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGAc                                                                                                                              >  1:2252769/1‑150 (MQ=255)
                        aaaaCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGAcc                                                                                                                             >  2:4222019/1‑150 (MQ=255)
                         aaaCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCt                                                                                                                            <  1:3419884/150‑1 (MQ=255)
                         aaaCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCt                                                                                                                            >  1:43982/1‑150 (MQ=255)
                           aCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTgg                                                                                                                          >  1:3135381/1‑150 (MQ=255)
                            cATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTggg                                                                                                                         >  1:2386124/1‑150 (MQ=255)
                             aTCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGt                                                                                                                        >  1:2098553/1‑150 (MQ=255)
                             aTCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATGTAGTGGATGTTCAAAAAGGACCTGGGt                                                                                                                        >  1:2729658/1‑150 (MQ=255)
                              tCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTc                                                                                                                       >  1:282883/1‑150 (MQ=255)
                               ccAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCt                                                                                                                      >  2:2958578/1‑150 (MQ=255)
                               ccAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCt                                                                                                                      >  2:786427/1‑150 (MQ=255)
                                 aGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTATCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGc                                                                                                                    >  1:1918638/1‑150 (MQ=255)
                                   tAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCAt                                                                                                                  <  2:3608369/150‑1 (MQ=255)
                                   tAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCAt                                                                                                                  >  2:1857017/1‑150 (MQ=255)
                                    aaTTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCAtt                                                                                                                 >  1:3640205/1‑150 (MQ=255)
                                    aaTTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCAtt                                                                                                                 <  1:3423412/150‑1 (MQ=255)
                                       tCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACg                                                                                                              <  1:3482416/150‑1 (MQ=255)
                                       tCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCTTTACg                                                                                                              >  2:2528686/1‑150 (MQ=255)
                                        cACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACgg                                                                                                             <  1:3416710/150‑1 (MQ=255)
                                        cACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACgg                                                                                                             >  2:1256789/1‑150 (MQ=255)
                                         acagACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGt                                                                                                            >  1:2480928/4‑150 (MQ=255)
                                         aCTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGt                                                                                                            <  1:2997769/150‑1 (MQ=255)
                                         aCTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGt                                                                                                            >  2:2870941/1‑150 (MQ=255)
                                          cTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGtt                                                                                                           <  2:1898894/150‑1 (MQ=255)
                                           tGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGttt                                                                                                          <  1:1725483/150‑1 (MQ=255)
                                             aCATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGttttt                                                                                                        >  2:2077050/1‑150 (MQ=255)
                                              cATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTc                                                                                                       >  2:3306572/1‑150 (MQ=255)
                                              cATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTc                                                                                                       >  2:1491634/1‑150 (MQ=255)
                                               aTTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTcc                                                                                                      <  2:3150889/150‑1 (MQ=255)
                                               aTTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTcc                                                                                                      <  2:1067389/150‑1 (MQ=255)
                                                 tACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCtg                                                                                                    <  2:1020578/150‑1 (MQ=255)
                                                         aaaTGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTAAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACaa                                                                                            <  2:112917/150‑1 (MQ=255)
                                                          aaTGCCGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTGCCTGTGAGACaca                                                                                           >  2:3746006/1‑148 (MQ=255)
                                                          aaTGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACaaa                                                                                           >  1:856437/1‑150 (MQ=255)
                                                            tGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACt                                                                                         >  2:1515982/1‑150 (MQ=255)
                                                            tGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACt                                                                                         <  1:4257449/150‑1 (MQ=255)
                                                             gCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTg                                                                                        >  1:1870803/1‑150 (MQ=255)
                                                                   cgcCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGc                                                                                  >  2:3652083/1‑150 (MQ=255)
                                                                   cgcCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGc                                                                                  >  1:325001/1‑150 (MQ=255)
                                                                      cAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAg                                                                               >  2:3951779/1‑150 (MQ=255)
                                                                       aGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGc                                                                              <  1:3896137/150‑1 (MQ=255)
                                                                        gCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCt                                                                             <  1:2783205/150‑1 (MQ=255)
                                                                              tCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAAt                                                                       >  1:3848482/1‑150 (MQ=255)
                                                                              tCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAAt                                                                       <  2:2143669/150‑1 (MQ=255)
                                                                              tCAACGATGGGTAAATTAACGGTGAGTTGCTGAAACAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAAt                                                                       <  2:3026052/150‑1 (MQ=255)
                                                                               cAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATa                                                                      >  1:1092958/1‑150 (MQ=255)
                                                                               cAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATa                                                                      >  2:2473176/1‑150 (MQ=255)
                                                                               cAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATa                                                                      >  2:340325/1‑150 (MQ=255)
                                                                                aaCGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAg                                                                     >  2:860474/1‑150 (MQ=255)
                                                                                 aCGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAgt                                                                    <  1:3198598/150‑1 (MQ=255)
                                                                                  cGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAgtg                                                                   >  2:2273847/1‑150 (MQ=255)
                                                                                      gggTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTgtgcg                                                               <  1:3786809/150‑1 (MQ=255)
                                                                                        gTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTgtgcgtg                                                             >  2:593222/1‑150 (MQ=255)
                                                                                           aaTTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAg                                                          >  1:1775472/1‑150 (MQ=255)
                                                                                           aaTTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAg                                                          >  1:4246707/1‑150 (MQ=255)
                                                                                            aTTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAgg                                                         >  2:2591556/1‑150 (MQ=255)
                                                                                             ttAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGc                                                        <  1:687002/150‑1 (MQ=255)
                                                                                             ttAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGc                                                        >  2:1998491/1‑150 (MQ=255)
                                                                                                  ggTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCa                                                   <  2:4172034/150‑1 (MQ=255)
                                                                                                  ggTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCa                                                   >  2:3899809/1‑150 (MQ=255)
                                                                                                       gTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATc                                              <  1:3739488/150‑1 (MQ=255)
                                                                                                       gTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGAGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTATTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATATTGTGCGAGCAGGCGGGCAACATc                                              >  1:1173947/1‑150 (MQ=255)
                                                                                                          gCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGa                                           <  1:246252/150‑1 (MQ=255)
                                                                                                             gAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGc                                        <  2:4341498/150‑1 (MQ=255)
                                                                                                                  aaaaaCTCACTGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAg                                   <  2:3358798/150‑1 (MQ=255)
                                                                                                                  aaaaaCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAg                                   >  2:2418270/1‑150 (MQ=255)
                                                                                                                  aaaaaCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAg                                   >  1:3053857/1‑150 (MQ=255)
                                                                                                                     aaCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCAt                                >  2:843494/1‑150 (MQ=255)
                                                                                                                      aCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATc                               <  2:1351312/150‑1 (MQ=255)
                                                                                                                      aCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATc                               <  1:4012861/150‑1 (MQ=255)
                                                                                                                         cAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCgg                            >  1:1443275/1‑150 (MQ=255)
                                                                                                                          aaTGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGt                           >  2:2648022/1‑150 (MQ=255)
                                                                                                                           aTGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTg                          <  1:2086987/150‑1 (MQ=255)
                                                                                                                               gcAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGcc                      >  2:713239/1‑150 (MQ=255)
                                                                                                                                cAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCg                     >  1:1200368/1‑150 (MQ=255)
                                                                                                                                   ccATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGtta                  >  1:1672684/1‑150 (MQ=255)
                                                                                                                                    cATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGttat                 <  2:3714295/150‑1 (MQ=255)
                                                                                                                                    cATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGttat                 <  1:1027470/150‑1 (MQ=255)
                                                                                                                                       gTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTtattat              >  1:3548139/1‑150 (MQ=255)
                                                                                                                                        tACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATc             >  1:268631/1‑150 (MQ=255)
                                                                                                                                         aCGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCa            >  1:2960923/1‑150 (MQ=255)
                                                                                                                                          cGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCaa           <  2:3866793/150‑1 (MQ=255)
                                                                                                                                            cATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTCTTCTCAATg         >  2:2921241/1‑150 (MQ=255)
                                                                                                                                              ttttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGAc       <  1:3681222/150‑1 (MQ=255)
                                                                                                                                                  ttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCa   >  1:2228943/1‑150 (MQ=255)
                                                                                                                                                   tCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAt  >  1:517986/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAT  >  CP009273/108902‑109198

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: ‑