Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,933 G→A 74.9% T55I (ACC→ATC)  ypfJ ← putative neutral zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,589,9330GA74.9% 417.8 / 154.1 227T55I (ACC→ATC) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base A (87/83);  total (120/107)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 227 total aligned reads displayed.

GGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGT  >  CP009273/2589785‑2590079
                                                                                                                                                    |                                                                                                                                                  
 gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGtt                                                                                                                                                  >  2:2535867/1‑150 (MQ=255)
  tCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTa                                                                                                                                                 >  1:2296899/1‑150 (MQ=255)
  tCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTa                                                                                                                                                 <  2:5134806/150‑1 (MQ=255)
     tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                                                                                                                                              <  1:2357703/150‑1 (MQ=255)
     tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTGGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                                                                                                                                              <  2:3156959/150‑1 (MQ=255)
      aCCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGTTTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATc                                                                                                                                             <  1:2915306/150‑1 (MQ=255)
      aCCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATc                                                                                                                                             >  2:4194699/1‑150 (MQ=255)
          aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAaca                                                                                                                                         >  1:3738359/1‑150 (MQ=255)
          aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAaca                                                                                                                                         >  2:2569049/1‑150 (MQ=255)
          aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAaca                                                                                                                                         >  2:661764/1‑150 (MQ=255)
           tcttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAacac                                                                                                                                        <  1:4593504/150‑1 (MQ=255)
              tctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTTGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAt                                                                                                                                     <  2:2522834/150‑1 (MQ=255)
                 cGAACTGTTGTCCCCATGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt                                                                                                                                  <  1:3829623/150‑1 (MQ=255)
                 cGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt                                                                                                                                  <  2:2005965/150‑1 (MQ=255)
                  gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCAtagta                                                                                                                                 <  1:529812/150‑1 (MQ=255)
                   aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag                                                                                                                                >  2:3666098/1‑150 (MQ=255)
                   aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag                                                                                                                                >  2:3003888/1‑150 (MQ=255)
                     cTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGcc                                                                                                                              <  2:3838234/150‑1 (MQ=255)
                      tgttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCt                                                                                                                             <  2:1634133/150‑1 (MQ=255)
                      tgttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCt                                                                                                                             <  2:5041129/150‑1 (MQ=255)
                      tgttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCt                                                                                                                             <  2:5160246/150‑1 (MQ=255)
                        ttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTATCCTGc                                                                                                                           >  1:2948292/1‑150 (MQ=255)
                         tgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCa                                                                                                                          <  2:3462044/150‑1 (MQ=255)
                            ccccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAAcc                                                                                                                       >  1:3388998/1‑150 (MQ=255)
                            ccccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAAcc                                                                                                                       >  2:3121000/1‑150 (MQ=255)
                             cccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCa                                                                                                                      <  2:3289009/150‑1 (MQ=255)
                             cccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCa                                                                                                                      <  2:4937246/150‑1 (MQ=255)
                              ccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAg                                                                                                                     >  2:1237142/1‑150 (MQ=255)
                               cAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGc                                                                                                                    <  1:4447879/150‑1 (MQ=255)
                               cAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGc                                                                                                                    >  2:901863/1‑150 (MQ=255)
                                aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCa                                                                                                                   <  1:1030122/150‑1 (MQ=255)
                                aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCa                                                                                                                   <  2:2057594/150‑1 (MQ=255)
                                   tgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGa                                                                                                                >  1:1292608/1‑150 (MQ=255)
                                    gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGAc                                                                                                               >  2:1148200/1‑150 (MQ=255)
                                    gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGAc                                                                                                               >  2:685482/1‑150 (MQ=255)
                                    gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGAc                                                                                                               <  2:2994354/150‑1 (MQ=255)
                                     tCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACt                                                                                                              >  2:2746208/1‑150 (MQ=255)
                                     tCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACt                                                                                                              >  2:129419/1‑150 (MQ=255)
                                           gTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGt                                                                                                        <  1:3621388/150‑1 (MQ=255)
                                           gTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGt                                                                                                        <  1:2829130/150‑1 (MQ=255)
                                                  ccAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCa                                                                                                 >  2:3637723/1‑150 (MQ=255)
                                                   cAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAg                                                                                                <  2:3395390/150‑1 (MQ=255)
                                                    aGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGc                                                                                               >  1:4740624/1‑150 (MQ=255)
                                                       aTCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCg                                                                                            <  1:2521765/150‑1 (MQ=255)
                                                             gAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTa                                                                                      >  1:630664/1‑150 (MQ=255)
                                                             gAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTa                                                                                      <  2:4573664/150‑1 (MQ=255)
                                                              aGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAACTCAGCCCGCCTTTAc                                                                                     >  2:2974500/1‑150 (MQ=255)
                                                              aGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAc                                                                                     >  1:2207901/1‑150 (MQ=255)
                                                                gTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCg                                                                                   >  2:2703438/1‑150 (MQ=255)
                                                                 tGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGc                                                                                  >  1:2435674/1‑150 (MQ=255)
                                                                   aaTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                                >  1:3262925/1‑150 (MQ=255)
                                                                   aaTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                                >  1:1196037/1‑150 (MQ=255)
                                                                     ttttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTgg                                                                              <  1:24767/150‑1 (MQ=255)
                                                                       ttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCa                                                                            <  2:370874/150‑1 (MQ=255)
                                                                         gCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAgg                                                                          <  2:1479883/150‑1 (MQ=255)
                                                                          cGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGc                                                                         >  1:56532/1‑150 (MQ=255)
                                                                            gCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGa                                                                       >  1:3769865/1‑150 (MQ=255)
                                                                              ttCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaa                                                                     >  1:3862712/1‑150 (MQ=255)
                                                                              ttCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaa                                                                     >  1:3843371/1‑150 (MQ=255)
                                                                               tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaaa                                                                    >  2:1607788/1‑150 (MQ=255)
                                                                                cGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAc                                                                   <  1:4053120/150‑1 (MQ=255)
                                                                                 gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  <  2:1045189/150‑1 (MQ=255)
                                                                                 gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  >  2:927638/1‑150 (MQ=255)
                                                                                  tCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCg                                                                 <  2:1888713/150‑1 (MQ=255)
                                                                                    ttCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCggg                                                               >  2:282614/1‑150 (MQ=255)
                                                                                      cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAc                                                             <  2:647541/150‑1 (MQ=255)
                                                                                       aTTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAcc                                                            <  1:2134033/150‑1 (MQ=255)
                                                                                           gggCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGccc                                                        <  1:796801/150‑1 (MQ=255)
                                                                                             gCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCAt                                                      <  2:4669548/150‑1 (MQ=255)
                                                                                               tAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAg                                                    <  2:2222528/150‑1 (MQ=255)
                                                                                               tAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAg                                                    <  1:1491933/150‑1 (MQ=255)
                                                                                                aaTTGACCGCGTTGATTGTTGTTGGGACACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGa                                                   <  1:1944329/150‑1 (MQ=255)
                                                                                                aaTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGa                                                   <  1:1565351/150‑1 (MQ=255)
                                                                                                 aTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGAt                                                  <  2:832111/150‑1 (MQ=255)
                                                                                                   tGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATgg                                                >  1:2612792/1‑150 (MQ=255)
                                                                                                    gACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGa                                               >  1:487612/1‑150 (MQ=255)
                                                                                                     aCCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGac                                              >  2:2796431/1‑150 (MQ=255)
                                                                                                      ccGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacc                                             <  1:2569049/150‑1 (MQ=255)
                                                                                                       cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacca                                            >  2:3084155/1‑150 (MQ=255)
                                                                                                                tgttgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTg                                   <  2:3904883/150‑1 (MQ=255)
                                                                                                                  ttgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGtt                                 <  1:4744369/150‑1 (MQ=255)
                                                                                                                   tgttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTg                                <  1:1146520/150‑1 (MQ=255)
                                                                                                                   tgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTg                                <  2:2889054/150‑1 (MQ=255)
                                                                                                                        ggAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCt                           <  2:2365652/150‑1 (MQ=255)
                                                                                                                        ggAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCt                           <  1:3120995/150‑1 (MQ=255)
                                                                                                                         gAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTg                          >  2:4685215/1‑150 (MQ=255)
                                                                                                                            aCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCt                       <  2:1979448/150‑1 (MQ=255)
                                                                                                                             ccGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCtt                      <  2:2597378/150‑1 (MQ=255)
                                                                                                                              cGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTc                     >  2:5077547/1‑150 (MQ=255)
                                                                                                                              cGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTc                     >  2:4108942/1‑150 (MQ=255)
                                                                                                                              cGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTc                     >  1:4348528/1‑150 (MQ=255)
                                                                                                                              cGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTc                     <  1:996494/150‑1 (MQ=255)
                                                                                                                               ggCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCa                    >  1:303980/1‑150 (MQ=255)
                                                                                                                               ggCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCa                    >  2:355527/1‑150 (MQ=255)
                                                                                                                                 cTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAAc                  <  1:3191923/150‑1 (MQ=255)
                                                                                                                                  tGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACa                 <  1:4188373/150‑1 (MQ=255)
                                                                                                                                         tcatcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCAc          <  1:2277293/150‑1 (MQ=255)
                                                                                                                                            tcaACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACttt       <  1:2583571/150‑1 (MQ=255)
                                                                                                                                              aaCCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCa     >  1:3953855/1‑150 (MQ=255)
                                                                                                                                              aaCCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTGCAACATTGTCACTTTCa     >  2:5094767/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGT  >  CP009273/2589785‑2590079

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: