Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,808 G→T 46.9% Q97K (CAG→AAG)  ypfJ ← putative neutral zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,589,8080GT46.9% 216.9 / inf 310Q97K (CAG→AAG) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (93/71);  new base T (96/49);  total (190/120)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

AGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCT  >  CP009273/2589661‑2589956
                                                                                                                                                   |                                                                                                                                                    
aGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgtt                                                                                                                                                    <  1:2973966/150‑1 (MQ=255)
 gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgttg                                                                                                                                                   <  1:3576728/150‑1 (MQ=255)
 gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgttg                                                                                                                                                   >  2:4076178/1‑150 (MQ=255)
   aGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTc                                                                                                                                                 <  1:3344426/150‑1 (MQ=255)
    ggAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcc                                                                                                                                                >  1:1281637/1‑150 (MQ=255)
    ggAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcc                                                                                                                                                <  1:4434242/150‑1 (MQ=255)
    ggAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcc                                                                                                                                                >  2:2787711/1‑150 (MQ=255)
       gagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCa                                                                                                                                             <  2:2201440/150‑1 (MQ=255)
       gagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCa                                                                                                                                             >  1:2447952/1‑150 (MQ=255)
          aTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgt                                                                                                                                          >  2:3650559/1‑150 (MQ=255)
          aTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgt                                                                                                                                          <  2:3227641/150‑1 (MQ=255)
            cGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgtgt                                                                                                                                        >  1:1821902/1‑150 (MQ=255)
            cGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgtgt                                                                                                                                        <  1:2842884/150‑1 (MQ=255)
             gATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTc                                                                                                                                       >  1:1413536/1‑150 (MQ=255)
             gATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                       >  1:3535784/1‑150 (MQ=255)
             gATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                       >  1:4117925/1‑150 (MQ=255)
              atatAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGAACGCGTCCTTCCACGGTACATCACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTATTGACCCCAGGTGTCt                                                                                                                                      >  2:641898/1‑150 (MQ=255)
               tataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCtt                                                                                                                                     >  1:64909/1‑150 (MQ=255)
                ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTc                                                                                                                                    >  2:1860456/1‑150 (MQ=255)
                ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTc                                                                                                                                    <  2:140841/150‑1 (MQ=255)
                ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTc                                                                                                                                    >  2:4021372/1‑150 (MQ=255)
                ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTc                                                                                                                                    <  1:4541661/150‑1 (MQ=255)
                ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTc                                                                                                                                    >  1:1100617/1‑150 (MQ=255)
                   aaCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGt                                                                                                                                 <  1:1411403/150‑1 (MQ=255)
                    aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTg                                                                                                                                <  2:1297362/150‑1 (MQ=255)
                    aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTg                                                                                                                                <  2:3619358/150‑1 (MQ=255)
                      cGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGt                                                                                                                              >  1:1853189/1‑150 (MQ=255)
                       gTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGtt                                                                                                                             >  1:520744/1‑150 (MQ=255)
                       gTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGtt                                                                                                                             >  2:3352997/1‑150 (MQ=255)
                       gTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGtt                                                                                                                             >  1:883334/1‑150 (MQ=255)
                        tGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTg                                                                                                                            >  1:1024381/1‑150 (MQ=255)
                         gCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGc                                                                                                                           <  2:729050/150‑1 (MQ=255)
                            aTCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAg                                                                                                                        >  2:3624154/1‑150 (MQ=255)
                              ccgccgGGCAATAGAACGGCCCCAGTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGaa                                                                                                                      <  2:1657389/150‑1 (MQ=255)
                                gccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATc                                                                                                                    <  2:2416584/150‑1 (MQ=255)
                                 ccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCa                                                                                                                   >  1:3329784/1‑150 (MQ=255)
                                  cgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCAc                                                                                                                  >  1:2311590/1‑150 (MQ=255)
                                  cgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCAc                                                                                                                  <  2:1725969/150‑1 (MQ=255)
                                      cAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAg                                                                                                              >  1:3439191/1‑150 (MQ=255)
                                       aaTAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAgg                                                                                                             <  1:1697964/150‑1 (MQ=255)
                                         tagccccgCCCCCTTCTGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGAAAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTg                                                                                                           <  2:4198462/143‑1 (MQ=255)
                                         tAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTg                                                                                                           >  2:3819027/1‑150 (MQ=255)
                                            aaCGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAt                                                                                                        <  1:816336/150‑1 (MQ=255)
                                              cGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAttt                                                                                                      >  1:2655431/1‑150 (MQ=255)
                                                  cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCg                                                                                                  >  1:3272675/1‑150 (MQ=255)
                                                       tATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTc                                                                                             >  2:1449153/1‑150 (MQ=255)
                                                          ggACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTc                                                                                          >  2:1965891/1‑150 (MQ=255)
                                                          ggACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTc                                                                                          <  1:2173076/150‑1 (MQ=255)
                                                              tGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCa                                                                                      >  2:2925416/1‑150 (MQ=255)
                                                                 cccgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAttt                                                                                   <  1:3094265/150‑1 (MQ=255)
                                                                  ccgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTg                                                                                  <  2:565912/150‑1 (MQ=255)
                                                                   cgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTgg                                                                                 >  2:2550793/1‑150 (MQ=255)
                                                                   cgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTgg                                                                                 >  2:2550717/1‑150 (MQ=255)
                                                                      cccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCt                                                                              <  2:4055117/150‑1 (MQ=255)
                                                                      cccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCt                                                                              >  1:3936113/1‑150 (MQ=255)
                                                                              cGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAcc                                                                      >  2:3102339/1‑150 (MQ=255)
                                                                                gTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                                    >  2:4121829/1‑150 (MQ=255)
                                                                                 tACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgcg                                                                   >  2:2636790/1‑150 (MQ=255)
                                                                                     cgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGa                                                               <  1:2696367/150‑1 (MQ=255)
                                                                                     cgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGa                                                               >  1:1771768/1‑150 (MQ=255)
                                                                                       tCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAGTTGGGCTAATTGACCGCGTTGAtt                                                             >  2:1792185/1‑150 (MQ=255)
                                                                                        cATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttg                                                            >  2:3912790/1‑150 (MQ=255)
                                                                                        cATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttg                                                            >  2:1219604/1‑150 (MQ=255)
                                                                                         aTTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgt                                                           >  2:3496394/1‑150 (MQ=255)
                                                                                          ttCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgtt                                                          <  2:3885089/150‑1 (MQ=255)
                                                                                          ttCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgtt                                                          >  1:2508708/1‑150 (MQ=255)
                                                                                           tCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttg                                                         >  2:723347/1‑150 (MQ=255)
                                                                                             cACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgtt                                                       >  1:505022/1‑150 (MQ=255)
                                                                                              aCGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                                      <  2:1023308/150‑1 (MQ=255)
                                                                                                 gTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGa                                                   >  2:3403305/1‑150 (MQ=255)
                                                                                                 gTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGa                                                   <  2:4508928/150‑1 (MQ=255)
                                                                                                   aCATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGaaa                                                 >  1:2742311/1‑150 (MQ=255)
                                                                                                    cATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAc                                                >  2:3168409/1‑150 (MQ=255)
                                                                                                    cATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAc                                                >  2:2355638/1‑150 (MQ=255)
                                                                                                      tGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCg                                              <  1:1974621/150‑1 (MQ=255)
                                                                                                          cAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCGTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTg                                          <  2:2889437/150‑1 (MQ=255)
                                                                                                            gCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGcc                                        >  1:695077/1‑150 (MQ=255)
                                                                                                                 ggTTGCTGATAGGTCTTACCCCTCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtc                                   >  2:3434047/1‑150 (MQ=255)
                                                                                                                 ggTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtc                                   >  2:3433858/1‑150 (MQ=255)
                                                                                                                 ggTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtc                                   >  2:3433964/1‑150 (MQ=255)
                                                                                                                   ttGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                                 >  1:3093475/1‑150 (MQ=255)
                                                                                                                    tGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcatca                               >  1:3098966/1‑150 (MQ=255)
                                                                                                                         aTAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAccc                           <  2:2729301/150‑1 (MQ=255)
                                                                                                                         aTAGGTCTTACCCATCGTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAccc                           <  2:1476180/150‑1 (MQ=255)
                                                                                                                                 tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   >  2:157646/1‑150 (MQ=255)
                                                                                                                                 tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   <  1:949475/150‑1 (MQ=255)
                                                                                                                                 tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   >  2:4485225/1‑150 (MQ=255)
                                                                                                                                     cATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAac               >  2:3722186/1‑150 (MQ=255)
                                                                                                                                     cATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAac               >  1:214383/1‑150 (MQ=255)
                                                                                                                                     cATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAac               >  1:213670/1‑150 (MQ=255)
                                                                                                                                      aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAaca              <  1:2782772/150‑1 (MQ=255)
                                                                                                                                           ctcGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAta         >  1:3765872/1‑150 (MQ=255)
                                                                                                                                           ctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAta         <  2:3272721/150‑1 (MQ=255)
                                                                                                                                           ctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAta         <  2:766080/150‑1 (MQ=255)
                                                                                                                                            tcGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGACAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACACGAGTAAATCAACACCAtag        >  2:3065427/1‑150 (MQ=255)
                                                                                                                                              gAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagta      <  2:1139545/150‑1 (MQ=255)
                                                                                                                                               aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag     >  2:2007134/1‑150 (MQ=255)
                                                                                                                                               aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag     >  2:4254289/1‑150 (MQ=255)
                                                                                                                                                 cTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGcc   <  1:1243912/150‑1 (MQ=255)
                                                                                                                                                  tgttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCt  >  2:3036963/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCT  >  CP009273/2589661‑2589956

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: