Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,346,661 T→A 6.1% P77P (CCA→CCT sapF ← antimicrobial peptide transport ABC system ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,346,6610TA 6.1% 1013.7 / 41.5 297P77P (CCA→CCTsapFantimicrobial peptide transport ABC system ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (137/140);  new base A (7/11);  total (146/151)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 308 total aligned reads displayed.

CGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTACCGCTTC  >  CP009273/1346512‑1346810
                                                                                                                                                     |                                                                                                                                                     
cGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGt                                                                                                                                                       <  1:2176045/150‑1 (MQ=255)
   cGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAGTACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGga                                                                                                                                                    >  2:3200420/1‑150 (MQ=255)
   cGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGga                                                                                                                                                    <  2:4807582/150‑1 (MQ=255)
     cTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGgatg                                                                                                                                                  <  1:5066463/150‑1 (MQ=255)
     cTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGgatg                                                                                                                                                  >  1:2113888/1‑150 (MQ=255)
         tCCGGGTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCg                                                                                                                                              <  1:4203491/150‑1 (MQ=255)
           cGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGtc                                                                                                                                            >  2:2113558/1‑150 (MQ=255)
            ggTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGtca                                                                                                                                           >  1:3182708/1‑150 (MQ=255)
               tCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGtcaatc                                                                                                                                        >  2:156548/1‑150 (MQ=255)
               tCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGtcaatc                                                                                                                                        <  1:1199630/150‑1 (MQ=255)
               tCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGtcaatc                                                                                                                                        <  2:4358313/150‑1 (MQ=255)
               tCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGtcaatc                                                                                                                                        <  2:117267/150‑1 (MQ=255)
                 cAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGtcaatcaa                                                                                                                                      <  2:4498730/150‑1 (MQ=255)
                  aGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTcaatcaat                                                                                                                                     <  1:4867264/150‑1 (MQ=255)
                  aGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTcaatcaat                                                                                                                                     <  2:3012278/150‑1 (MQ=255)
                  aGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTcaatcaat                                                                                                                                     >  1:3588436/1‑150 (MQ=255)
                      cGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACt                                                                                                                                 >  1:2331664/1‑150 (MQ=255)
                        gtgtTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTcg                                                                                                                               >  1:756286/1‑150 (MQ=255)
                         tgtTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTcgc                                                                                                                              <  1:1443834/150‑1 (MQ=255)
                           ttCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTcgccg                                                                                                                            >  2:823462/1‑150 (MQ=255)
                                 cgcAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTa                                                                                                                      >  2:543923/1‑150 (MQ=255)
                                 cgcAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTa                                                                                                                      <  1:3233292/150‑1 (MQ=255)
                                  gcAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAg                                                                                                                     >  1:1841101/1‑150 (MQ=255)
                                    aGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGc                                                                                                                   <  2:2648803/150‑1 (MQ=255)
                                      tGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGctc                                                                                                                 <  2:3413683/150‑1 (MQ=255)
                                      tGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGctc                                                                                                                 >  1:297287/1‑150 (MQ=255)
                                       ggAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGctct                                                                                                                <  2:2154931/150‑1 (MQ=255)
                                       ggAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGctct                                                                                                                <  2:278738/150‑1 (MQ=255)
                                        gAACATCCAGACTGTGCGAGAGACGTTGACGGGGATTCACCGAGGTCGAAGGATCCTGAAACATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTa                                                                                                               <  2:3168923/150‑1 (MQ=255)
                                        gAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTa                                                                                                               >  2:2532967/1‑150 (MQ=255)
                                        gAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTa                                                                                                               >  2:3575530/1‑150 (MQ=255)
                                         aaaaTCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTAt                                                                                                              >  1:3338658/1‑150 (MQ=255)
                                          aaaTCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATc                                                                                                             >  1:967597/1‑150 (MQ=255)
                                            aTCCAGAATGTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCAt                                                                                                           <  2:3126704/150‑1 (MQ=255)
                                               cAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTcc                                                                                                        <  1:3829263/150‑1 (MQ=255)
                                                 gAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCg                                                                                                      <  1:3453383/150‑1 (MQ=255)
                                                   aTTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGcc                                                                                                    >  2:2281104/1‑150 (MQ=255)
                                                      tGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGc                                                                                                 <  1:2060018/150‑1 (MQ=255)
                                                      tGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGc                                                                                                 >  1:5005800/1‑150 (MQ=255)
                                                         gagaTACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCAtt                                                                                              >  1:2782912/1‑150 (MQ=255)
                                                          agaTACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCAttt                                                                                             >  1:4705799/1‑150 (MQ=255)
                                                           gaTACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCAtttt                                                                                            >  1:4539463/1‑150 (MQ=255)
                                                            aTACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCAttttt                                                                                           >  2:1425029/1‑150 (MQ=255)
                                                            aTACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCAttttt                                                                                           <  1:1128974/150‑1 (MQ=255)
                                                              aCGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGc                                                                                         <  1:4098500/150‑1 (MQ=255)
                                                               cGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGcc                                                                                        >  2:200255/1‑150 (MQ=255)
                                                                gTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCa                                                                                       <  1:179385/150‑1 (MQ=255)
                                                                  tGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGc                                                                                     >  2:4058257/1‑150 (MQ=255)
                                                                   gACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCg                                                                                    <  2:4105689/150‑1 (MQ=255)
                                                                     cGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTg                                                                                  >  2:1144590/1‑150 (MQ=255)
                                                                     cGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTg                                                                                  >  2:2388883/1‑150 (MQ=255)
                                                                      gggtATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTgg                                                                                 >  2:3979611/1‑150 (MQ=255)
                                                                      ggggATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTgg                                                                                 <  1:251610/150‑1 (MQ=255)
                                                                         gATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGAtt                                                                              <  2:2908040/150‑1 (MQ=255)
                                                                                  aGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGacc                                                                     <  1:1916262/150‑3 (MQ=255)
                                                                                  aGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAAc                                                                     >  2:1023873/1‑150 (MQ=255)
                                                                                  aGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAAc                                                                     <  2:4149129/150‑1 (MQ=255)
                                                                                  aGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAAc                                                                     >  2:1851534/1‑150 (MQ=255)
                                                                                     tCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCAt                                                                  >  1:523822/1‑150 (MQ=255)
                                                                                       gAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATtc                                                                >  2:1383327/1‑150 (MQ=255)
                                                                                       gAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATtc                                                                >  2:2827976/1‑150 (MQ=255)
                                                                                        aaGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATtct                                                               >  1:2535545/1‑150 (MQ=255)
                                                                                         aGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATtctc                                                              >  2:3849394/1‑150 (MQ=255)
                                                                                         aGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATtctc                                                              <  2:691813/150‑1 (MQ=255)
                                                                                           gATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGc                                                            >  1:667011/1‑150 (MQ=255)
                                                                                               cTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCaat                                                        <  2:4438401/150‑1 (MQ=255)
                                                                                               cTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCaat                                                        <  1:2375948/150‑1 (MQ=255)
                                                                                                 gAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCaataa                                                      <  2:3064060/150‑1 (MQ=255)
                                                                                                  aaaaaTCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCaataat                                                     <  1:465404/150‑1 (MQ=255)
                                                                                                  aaaaaTCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCaataat                                                     <  2:350534/150‑1 (MQ=255)
                                                                                                      aTCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGcc                                                 <  1:1049775/150‑1 (MQ=255)
                                                                                                      aTCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGcc                                                 <  2:4414387/150‑1 (MQ=255)
                                                                                                       tCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGAGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCa                                                >  2:4274387/1‑150 (MQ=255)
                                                                                                         aTGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAgt                                              >  1:3544388/1‑150 (MQ=255)
                                                                                                                tACGCTTACTTCGGAAGGAATAATCGCCAAAAAGCAGTGGATGATCGACAATCAAAAACTCGACGCTGGTAGGCTCTAACATTCCCGCCAGCATTTTAGCCAGCGTGGATTTACCAGAACCAATCTCGCCAATAATCACCAGAGTCTgtc                                       >  1:2917337/1‑148 (MQ=255)
                                                                                                                  cGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCt                                     <  2:1297924/150‑1 (MQ=255)
                                                                                                                  cGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCt                                     <  2:1297835/150‑1 (MQ=255)
                                                                                                                     tGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCa                                  <  2:1029578/150‑1 (MQ=255)
                                                                                                                       aCTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACg                                >  1:2432943/1‑150 (MQ=255)
                                                                                                                       aCTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACg                                <  1:1293279/150‑1 (MQ=255)
                                                                                                                           cGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGc                            <  2:3603468/150‑1 (MQ=255)
                                                                                                                           cGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGc                            >  1:4447183/1‑150 (MQ=255)
                                                                                                                            ggaaggaaTAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCg                           >  2:1287257/1‑150 (MQ=255)
                                                                                                                            ggaaggaaTAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCg                           >  1:1730424/1‑150 (MQ=255)
                                                                                                                              aaggaaTAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTa                         >  2:3303972/1‑150 (MQ=255)
                                                                                                                                 gaaTAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAg                      <  2:1678532/150‑1 (MQ=255)
                                                                                                                                  aataatCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGc                     >  2:3316472/1‑150 (MQ=255)
                                                                                                                                  aataatCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGc                     <  2:2331664/150‑1 (MQ=255)
                                                                                                                                     aatCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCa                  >  1:2683627/1‑150 (MQ=255)
                                                                                                                                           cAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGtt            <  2:830867/150‑1 (MQ=255)
                                                                                                                                           cAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGtt            >  1:3764963/1‑150 (MQ=255)
                                                                                                                                            aaaaTGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGttt           <  2:3686100/150‑1 (MQ=255)
                                                                                                                                            aaaaTGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGttt           >  2:4902916/1‑150 (MQ=255)
                                                                                                                                             aaaTGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGtttt          <  2:1372205/150‑1 (MQ=255)
                                                                                                                                                tGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTAcc       >  1:4603865/1‑150 (MQ=255)
                                                                                                                                                tGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTAcc       <  2:3443020/150‑1 (MQ=255)
                                                                                                                                                 gCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTACCg      <  1:134468/150‑1 (MQ=255)
                                                                                                                                                     tGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTACCGCTTc  >  1:4004381/1‑150 (MQ=255)
                                                                                                                                                     tGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTACCGCTTc  <  1:4932335/150‑1 (MQ=255)
                                                                                                                                                     tGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTACCGCTTc  >  2:3772951/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTTGCGAGATACGTTGACGGGGATTCAACGAGGTCGAAGGATCCTGAAAAATCATGCGAATACGCTGACTGCGGAAGGAATAATCGCCAAAATGCAGTGGATGATCGTCAATCAATAACTCGCCGCTGGTAGGCTCTATCATTCCCGCCAGCATTTTTGCCAGCGTGGATTTACCCGAACCATTCTCGCCAATAATCGCCAGTGTCTGGCCTTCACGTAGCGTAAAGCTCAAGGGTTTTACCGCTTC  >  CP009273/1346512‑1346810

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: