Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,933 G→A 85.3% 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,9330GA85.3% 693.8 / 100.5 251T55I (ACC→ATC) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (23/14);  new base A (109/105);  total (132/119)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 253 total aligned reads displayed.

AGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTC  >  CP009273/2589784‑2590080
                                                                                                                                                     |                                                                                                                                                   
aGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGGCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCgg                                                                                                                                                     >  2:755733/1‑150 (MQ=255)
  gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                                                                                                                                                   <  1:3840958/150‑1 (MQ=255)
  gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                                                                                                                                                   >  2:432651/1‑150 (MQ=255)
      tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                                                                                                                                               <  2:4456230/150‑1 (MQ=255)
       aCCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCCTCAACCCGGTTAAATc                                                                                                                                              >  2:4181483/1‑150 (MQ=255)
        cccATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCa                                                                                                                                             <  1:3724980/150‑1 (MQ=255)
          cATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAac                                                                                                                                           >  2:2467814/1‑150 (MQ=255)
           aTCTTCTCGACCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAaca                                                                                                                                          <  2:2639349/150‑1 (MQ=255)
           aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAaca                                                                                                                                          >  1:1275601/1‑150 (MQ=255)
               tctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAt                                                                                                                                      >  1:763446/1‑150 (MQ=255)
                ctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAta                                                                                                                                     <  2:1350064/150‑1 (MQ=255)
                ctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAta                                                                                                                                     <  2:1227818/150‑1 (MQ=255)
                 tcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag                                                                                                                                    <  2:1351250/150‑1 (MQ=255)
                  cGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt                                                                                                                                   <  1:3961270/150‑1 (MQ=255)
                   gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagta                                                                                                                                  >  1:1536270/1‑150 (MQ=255)
                   gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagta                                                                                                                                  >  1:2974626/1‑150 (MQ=255)
                    aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag                                                                                                                                 >  1:4108730/1‑150 (MQ=255)
                       tgttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCt                                                                                                                              >  1:993389/1‑150 (MQ=255)
                        gttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTg                                                                                                                             <  1:3478399/150‑1 (MQ=255)
                        gttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTg                                                                                                                             <  1:3174879/150‑1 (MQ=255)
                        gttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTg                                                                                                                             >  2:4275680/1‑150 (MQ=255)
                         ttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGc                                                                                                                            >  1:985721/1‑150 (MQ=255)
                          tgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCa                                                                                                                           <  2:928641/150‑1 (MQ=255)
                            tCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAAc                                                                                                                         >  2:2136502/1‑150 (MQ=255)
                              cccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCa                                                                                                                       <  2:1092604/150‑1 (MQ=255)
                                 aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCa                                                                                                                    <  1:59601/150‑1 (MQ=255)
                                  ggTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCAc                                                                                                                   <  2:3200956/150‑1 (MQ=255)
                                     gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGAc                                                                                                                >  2:3855468/1‑150 (MQ=255)
                                     gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGAc                                                                                                                >  1:3947566/1‑150 (MQ=255)
                                      tCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACt                                                                                                               >  2:3010992/1‑150 (MQ=255)
                                        ttCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTAt                                                                                                             <  1:2724101/150‑1 (MQ=255)
                                         tCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATc                                                                                                            <  1:4276561/150‑1 (MQ=255)
                                          ccGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCa                                                                                                           <  2:379389/150‑1 (MQ=255)
                                           cGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAg                                                                                                          >  1:4057126/1‑150 (MQ=255)
                                           cGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGAAGCCTGCAACCAGCACGACAATCAg                                                                                                          >  1:1219452/1‑150 (MQ=255)
                                            gTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGt                                                                                                         >  1:2194961/1‑150 (MQ=255)
                                              ggTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTaa                                                                                                       <  1:859978/150‑1 (MQ=255)
                                              ggTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTaa                                                                                                       <  2:4185158/150‑1 (MQ=255)
                                                 tGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAAt                                                                                                    <  1:3036719/150‑1 (MQ=255)
                                                    cAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAg                                                                                                 >  2:500673/1‑150 (MQ=255)
                                                      gAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGcc                                                                                               >  1:1675464/1‑150 (MQ=255)
                                                       aaTCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGccc                                                                                              >  1:1165967/1‑150 (MQ=255)
                                                          cACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGcc                                                                                           >  2:2832434/1‑150 (MQ=255)
                                                           aCCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCt                                                                                          <  1:1655570/150‑1 (MQ=255)
                                                            ccGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCtt                                                                                         <  2:2653420/150‑1 (MQ=255)
                                                            ccGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCtt                                                                                         >  1:2233335/1‑150 (MQ=255)
                                                            ccGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCtt                                                                                         >  1:235854/1‑150 (MQ=255)
                                                             cGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCttt                                                                                        <  2:2893738/150‑1 (MQ=255)
                                                             cGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCttt                                                                                        <  2:1116987/150‑1 (MQ=255)
                                                               aGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAc                                                                                      >  1:677438/1‑150 (MQ=255)
                                                               aGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAc                                                                                      >  1:176451/1‑150 (MQ=255)
                                                               aGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAGTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAc                                                                                      >  1:4368286/1‑150 (MQ=255)
                                                                ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCATTAAAATCATCCCGCCTTTAcc                                                                                     <  2:3644503/150‑1 (MQ=255)
                                                                ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAcc                                                                                     >  2:1671672/1‑150 (MQ=255)
                                                                ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAcc                                                                                     <  1:89929/150‑1 (MQ=255)
                                                                 gTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCg                                                                                    >  1:3516091/1‑150 (MQ=255)
                                                                  tGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGc                                                                                   >  1:3796815/1‑150 (MQ=255)
                                                                   gAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCt                                                                                  <  1:3754226/150‑1 (MQ=255)
                                                                    aaTTTTGCGGCGTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                                 <  2:1661231/150‑1 (MQ=255)
                                                                     aTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTg                                                                                >  1:2543030/1‑150 (MQ=255)
                                                                         tGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAg                                                                            >  2:1965660/1‑150 (MQ=255)
                                                                            ggCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCg                                                                         <  1:322896/150‑1 (MQ=255)
                                                                            ggCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCg                                                                         >  1:3578221/1‑150 (MQ=255)
                                                                                tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaaa                                                                     >  1:4316653/1‑150 (MQ=255)
                                                                                 cGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAc                                                                    <  1:1911358/150‑1 (MQ=255)
                                                                                 cGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAc                                                                    <  2:1378085/150‑1 (MQ=255)
                                                                                 cGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAc                                                                    <  2:3457374/150‑1 (MQ=255)
                                                                                  gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                   >  1:3037266/1‑150 (MQ=255)
                                                                                    cTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCgg                                                                 <  2:3168810/150‑1 (MQ=255)
                                                                                    cTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCgg                                                                 <  1:36775/150‑1 (MQ=255)
                                                                                    cTTCATTTGGGCTAAGTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCgg                                                                 <  1:4479064/150‑1 (MQ=255)
                                                                                        aTTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAcc                                                             >  2:3819641/1‑150 (MQ=255)
                                                                                               cTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATa                                                      >  1:512220/1‑150 (MQ=255)
                                                                                                tAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAg                                                     >  2:3282505/1‑150 (MQ=255)
                                                                                                    tGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATgg                                                 >  1:2521143/1‑150 (MQ=255)
                                                                                                    tGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATgg                                                 >  1:1469538/1‑150 (MQ=255)
                                                                                                        cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacca                                             <  1:1769141/150‑1 (MQ=255)
                                                                                                         gcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGaccac                                            >  1:733409/1‑150 (MQ=255)
                                                                                                              gATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCagag                                       >  1:3404764/1‑150 (MQ=255)
                                                                                                              gATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCagag                                       <  1:101814/150‑1 (MQ=255)
                                                                                                              gATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCagag                                       >  1:1118399/1‑150 (MQ=255)
                                                                                                               aTTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGc                                      <  1:2650269/150‑1 (MQ=255)
                                                                                                                  gttgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGt                                   <  1:3352222/150‑1 (MQ=255)
                                                                                                                    tgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTg                                 >  2:1936992/1‑150 (MQ=255)
                                                                                                                     gttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgc                                >  1:777757/1‑150 (MQ=255)
                                                                                                                     gttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgc                                >  2:442459/1‑150 (MQ=255)
                                                                                                                         ggAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCt                            <  2:1193523/150‑1 (MQ=255)
                                                                                                                           aaaCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGt                          >  1:3024974/1‑150 (MQ=255)
                                                                                                                           aaaCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGt                          >  1:3025025/1‑150 (MQ=255)
                                                                                                                             aCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCt                        >  2:2191314/1‑150 (MQ=255)
                                                                                                                             aCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCt                        >  2:4364210/1‑150 (MQ=255)
                                                                                                                                     cccGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACAtt                >  2:3593112/1‑150 (MQ=255)
                                                                                                                                        ggTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTc             <  1:4522674/150‑1 (MQ=255)
                                                                                                                                         gTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCa            <  2:777629/150‑1 (MQ=255)
                                                                                                                                          tcatcaACCCGATTAAATCAACACCATAGTAGTCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCAc           >  1:1111025/1‑150 (MQ=255)
                                                                                                                                            atcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACtt         <  2:1506382/150‑1 (MQ=255)
                                                                                                                                            atcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACtt         <  2:1420372/150‑1 (MQ=255)
                                                                                                                                                  ccGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGt   >  1:3906673/1‑150 (MQ=255)
                                                                                                                                                   cGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTc  <  1:2613661/150‑1 (MQ=255)
                                                                                                                                                   aggTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTc  <  1:198534/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTC  >  CP009273/2589784‑2590080

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: