Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,808 G→T 53.8% 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,8080GT53.8% 51.9 / inf 267Q97K (CAG→AAG) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (67/56);  new base T (85/59);  total (152/115)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 284 total aligned reads displayed.

AGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCT  >  CP009273/2589661‑2589956
                                                                                                                                                   |                                                                                                                                                    
aGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGTTTGCTGATAGGTCTTACCCATCTTCTCGAACtgtt                                                                                                                                                    <  1:4182353/150‑1 (MQ=255)
 gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACttttg                                                                                                                                                   >  2:1422657/1‑150 (MQ=255)
  aaGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACttttgt                                                                                                                                                  <  2:770340/150‑1 (MQ=255)
  aaGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACttttgt                                                                                                                                                  >  2:1253452/1‑150 (MQ=255)
       gagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCa                                                                                                                                             >  2:4243879/1‑150 (MQ=255)
       gagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCa                                                                                                                                             >  2:1742220/1‑150 (MQ=255)
          aTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGgt                                                                                                                                          <  2:2576969/150‑1 (MQ=255)
           tCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGgtg                                                                                                                                         >  2:63062/1‑150 (MQ=255)
           tCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGgtg                                                                                                                                         >  2:4028650/1‑150 (MQ=255)
            cGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGgtgt                                                                                                                                        <  1:2496556/150‑1 (MQ=255)
            cGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgtgt                                                                                                                                        <  1:3209604/150‑1 (MQ=255)
             gATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                       >  1:1098196/1‑150 (MQ=255)
                  aaaCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCg                                                                                                                                  >  2:1136986/1‑150 (MQ=255)
                    aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTg                                                                                                                                <  2:4008907/150‑1 (MQ=255)
                    aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTg                                                                                                                                >  1:2211665/1‑150 (MQ=255)
                    aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTg                                                                                                                                >  2:312061/1‑150 (MQ=255)
                        tGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTg                                                                                                                            <  2:1877846/150‑1 (MQ=255)
                         gCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGc                                                                                                                           >  1:3151059/1‑150 (MQ=255)
                           cATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCa                                                                                                                         >  1:1527245/1‑150 (MQ=255)
                              ccgccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGaa                                                                                                                      >  2:3977263/1‑150 (MQ=255)
                              ccgccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGaa                                                                                                                      <  1:2059508/150‑1 (MQ=255)
                              ccgccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGaa                                                                                                                      <  1:2307805/150‑1 (MQ=255)
                                    ggCAATAGAACGGCCCCATTATGGACTGTCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCg                                                                                                                <  2:4533180/150‑1 (MQ=255)
                                    ggCAATAGAACGGCCCCATTATGGACTGTCCCGCCCCGCAGCCGGTACGCGTCAGTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCg                                                                                                                <  2:1019794/150‑1 (MQ=255)
                                    ggCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCg                                                                                                                <  2:3741834/150‑1 (MQ=255)
                                    ggCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCg                                                                                                                >  1:975845/1‑150 (MQ=255)
                                      cAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAg                                                                                                              >  1:4378134/1‑150 (MQ=255)
                                       aaTAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAgg                                                                                                             >  1:2480749/1‑150 (MQ=255)
                                        aTAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGt                                                                                                            >  1:3439276/1‑150 (MQ=255)
                                        aTAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGt                                                                                                            >  1:4074619/1‑150 (MQ=255)
                                         tAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTg                                                                                                           <  2:3924903/150‑1 (MQ=255)
                                          aGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGa                                                                                                          >  2:4586357/1‑150 (MQ=255)
                                               ggCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAtttt                                                                                                     >  2:2461924/1‑150 (MQ=255)
                                               ggCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAtttt                                                                                                     >  1:3531023/1‑150 (MQ=255)
                                               ggCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAtttt                                                                                                     <  1:1995712/150‑1 (MQ=255)
                                                  cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCg                                                                                                  >  1:1182866/1‑150 (MQ=255)
                                                  cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCg                                                                                                  <  1:1654421/150‑1 (MQ=255)
                                                  cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCg                                                                                                  <  2:4029163/150‑1 (MQ=255)
                                                      ttATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCtt                                                                                              <  1:580507/150‑1 (MQ=255)
                                                      ttATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCtt                                                                                              >  2:1363977/1‑150 (MQ=255)
                                                        aTGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCg                                                                                            >  1:3886802/1‑150 (MQ=255)
                                                        aTGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCg                                                                                            <  1:3398415/150‑1 (MQ=255)
                                                          ggACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTc                                                                                          >  2:1060033/1‑150 (MQ=255)
                                                             cTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTc                                                                                       <  1:2581872/150‑1 (MQ=255)
                                                              tGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCa                                                                                      >  2:1768228/1‑150 (MQ=255)
                                                               ggCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAt                                                                                     <  2:3193691/150‑1 (MQ=255)
                                                                 cccgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAttt                                                                                   <  2:534690/150‑1 (MQ=255)
                                                                   cgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTgg                                                                                 >  2:885926/1‑150 (MQ=255)
                                                                   cgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTgg                                                                                 >  1:1584866/1‑150 (MQ=255)
                                                                   cgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTgg                                                                                 >  1:3762995/1‑150 (MQ=255)
                                                                    gccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTggg                                                                                >  1:4308814/1‑150 (MQ=255)
                                                                    gccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTggg                                                                                >  1:2815359/1‑150 (MQ=255)
                                                                        cgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTaa                                                                            >  1:2655710/1‑150 (MQ=255)
                                                                         gcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAAt                                                                           <  2:165951/150‑1 (MQ=255)
                                                                         gcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAAt                                                                           >  2:2631450/1‑150 (MQ=255)
                                                                           aGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTg                                                                         >  1:3931208/1‑150 (MQ=255)
                                                                            gCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGa                                                                        <  2:2494582/150‑1 (MQ=255)
                                                                               ggTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcg                                                                     >  1:2253640/1‑150 (MQ=255)
                                                                                gTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                                    <  2:1570724/150‑1 (MQ=255)
                                                                                       tCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAtt                                                             <  2:330919/150‑1 (MQ=255)
                                                                                           tCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttg                                                         >  2:2115994/1‑150 (MQ=255)
                                                                                           tCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttg                                                         >  2:4205622/1‑150 (MQ=255)
                                                                                               cGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTgg                                                     <  1:1735001/150‑1 (MQ=255)
                                                                                                ggTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTggg                                                    >  1:213985/1‑150 (MQ=255)
                                                                                                  tACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGaa                                                  >  2:1945116/1‑150 (MQ=255)
                                                                                                   aCATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGaaa                                                 <  2:747314/150‑1 (MQ=255)
                                                                                                    cATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAc                                                >  1:1326711/1‑150 (MQ=255)
                                                                                                     aTGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAcc                                               <  1:3232240/150‑1 (MQ=255)
                                                                                                       gACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCgg                                             <  2:3562952/150‑1 (MQ=255)
                                                                                                         ccAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCt                                           <  2:2348032/150‑1 (MQ=255)
                                                                                                         ccAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCt                                           <  2:3283318/150‑1 (MQ=255)
                                                                                                           aGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGc                                         <  1:4486778/150‑1 (MQ=255)
                                                                                                             cTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                                       <  1:1253452/150‑1 (MQ=255)
                                                                                                             cTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                                       <  2:1221279/150‑1 (MQ=255)
                                                                                                             cTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                                       <  1:2019650/150‑1 (MQ=255)
                                                                                                                 ggTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtc                                   >  1:896898/1‑150 (MQ=255)
                                                                                                                 ggTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtc                                   >  1:1078438/1‑150 (MQ=255)
                                                                                                                  gTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtca                                  >  1:54379/1‑150 (MQ=255)
                                                                                                                   ttGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                                 >  2:2782726/1‑150 (MQ=255)
                                                                                                                   ttGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                                 <  1:721729/150‑1 (MQ=255)
                                                                                                                    tGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcatc                                >  2:3384114/1‑150 (MQ=255)
                                                                                                                     gCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcatca                               <  2:3864907/150‑1 (MQ=255)
                                                                                                                       tGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAc                             >  1:3904387/1‑150 (MQ=255)
                                                                                                                        gATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAcc                            <  1:1810728/150‑1 (MQ=255)
                                                                                                                         aTAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAccc                           <  2:4349449/150‑1 (MQ=255)
                                                                                                                          tAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCg                          <  1:831160/150‑1 (MQ=255)
                                                                                                                            ggTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCgat                        >  2:4396942/1‑148 (MQ=255)
                                                                                                                                 tACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   <  1:1363977/150‑1 (MQ=255)
                                                                                                                                 tACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   >  1:4386027/1‑150 (MQ=255)
                                                                                                                                 tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   <  2:2014067/150‑1 (MQ=255)
                                                                                                                                    ccATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCaa                <  1:4348069/150‑1 (MQ=255)
                                                                                                                                     cATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAac               <  1:1035950/150‑1 (MQ=255)
                                                                                                                                       tcttctCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAacac             >  2:4146416/1‑150 (MQ=255)
                                                                                                                                        cttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACAcc            <  1:312061/150‑1 (MQ=255)
                                                                                                                                         ttctCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCa           <  1:946871/150‑1 (MQ=255)
                                                                                                                                            tcGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag        <  1:3950207/150‑1 (MQ=255)
                                                                                                                                            tcGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag        >  1:2375328/1‑150 (MQ=255)
                                                                                                                                             cGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt       >  2:4467654/1‑150 (MQ=255)
                                                                                                                                              gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagta      >  1:373055/1‑150 (MQ=255)
                                                                                                                                                 cTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGcc   >  1:3502429/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: