Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG41.3% 144.6 / 11.9 76T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (6/38);  new base G (31/0);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 2.05e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCC  >  NC_000913/1356054‑1356342
                                                                                                                                                    |                                                                                                                                               
cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGtt                                                                                                                                              <  2:204826/151‑1 (MQ=255)
   tGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTaa                                                                                                                                            >  1:242137/1‑150 (MQ=255)
    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTaaa                                                                                                                                           <  2:858010/150‑1 (MQ=255)
       aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAAAtcc                                                                                                                                        >  2:975053/1‑149 (MQ=255)
       aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGcc                                                                                                                                       <  1:1237264/151‑1 (MQ=255)
          ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGa                                                                                                                                    <  2:961499/151‑1 (MQ=255)
            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGcc                                                                                                                                       >  1:1125433/1‑146 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTAAATGCCCGa                                                                                                                                    >  1:1339434/1‑147 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAgg                                                                                                                                <  2:619768/151‑1 (MQ=255)
               aCGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGATaag                                                                                                                                >  2:1276889/1‑148 (MQ=255)
                  aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGAc                                                                                                                            <  2:362487/151‑1 (MQ=255)
                  aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGAc                                                                                                                            <  2:294668/151‑1 (MQ=255)
                       cgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGCGGTAAAttacaaaaaggggg                                                                                                                             >  2:822110/1‑132 (MQ=255)
                         cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATc                                                                                                                     <  2:69094/151‑1 (MQ=255)
                             aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGTGTTTAAAGCCCGATAGGGGGCGCTCATccc                                                                                                                   >  2:1228475/1‑149 (MQ=255)
                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAGGCGGCGTTCATCCCCgg                                                                                                               >  1:115858/1‑151 (MQ=255)
                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGATAGGCGACGTTCATCCCCgg                                                                                                               >  1:1196647/1‑151 (MQ=255)
                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAAATAGGCGACGGTCATtcc                                                                                                                   >  1:1108513/1‑148 (MQ=255)
                                  agagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGc                                                                                                           <  2:568446/151‑1 (MQ=255)
                                       ccgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGt                                                                                                      <  1:889617/151‑1 (MQ=255)
                                         gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTa                                                                                                    <  1:498122/151‑1 (MQ=255)
                                         gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTa                                                                                                    <  2:1070807/151‑1 (MQ=255)
                                           cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTaaa                                                                                                  <  1:714449/151‑1 (MQ=255)
                                            gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCGGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAAc                                                                                                 <  2:632014/151‑1 (MQ=255)
                                            gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGCCAAATAGGCGGCGTTTATTCCCGGACGCAGCggttaa                                                                                                  >  1:241894/1‑145 (MQ=255)
                                               tCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAAGGGGCGTTTATCCCCGGACGCAGCGGTAAAAGc                                                                                               >  2:135815/1‑150 (MQ=255)
                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAATGCCCGAAAGGCGGCGTTCATCCCCGGACGCCGCGTTAAACGCaa                                                                                             >  2:934192/1‑150 (MQ=255)
                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAg                                                                                            <  1:1115804/151‑1 (MQ=255)
                                                  gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgc                                                                                           <  1:1382384/151‑1 (MQ=255)
                                                   ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgcg                                                                                          <  1:531300/151‑1 (MQ=255)
                                                    cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgcgc                                                                                         <  1:490522/151‑1 (MQ=255)
                                                     gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGcgcg                                                                                        <  1:842425/151‑1 (MQ=255)
                                                       aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAATGCCCGATAAGCGGCGCTCATCCCCGGACGCAGCGGTTAAAGCAAACGCggg                                                                                      >  2:1404273/1‑151 (MQ=255)
                                                         cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtc                                                                                    <  2:83980/151‑1 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGCCAGAAAGGCGCCGATCATCCCCCGACGCCGCGCTTAAAGCAAACGCgggggg                                                                                   >  1:362374/1‑148 (MQ=255)
                                                            aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtcgtc                                                                                 <  2:847015/151‑1 (MQ=255)
                                                            aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAAAGCCAGATAGGCGGCGGTCATCCCCGGACGCAGCGTTAAACGCAAACGCGGggcgt                                                                                  >  1:1157381/1‑150 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAAGCGACGTTTATCCCCCGACGCAACGTTTAACCCAAACCGGGGGCGTCAc                                                                               >  2:725471/1‑151 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATACCCGAAAAGGGGCGGTTATCCCCGGGCGCAGGGGTAAAAGCAAGCGCgggg                                                                                     >  2:881940/1‑144 (MQ=255)
                                                                  tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAATGCCCGATAGGCGACGTTCATCCCAAGACGCAGCGGTTAAAGCAAACGCggg                                                                                      >  1:691954/1‑140 (MQ=255)
                                                                    gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAATCCAGATAAGCGGCGGTCATCCCCGGGAGCCGCGGTTAAAGCCAACGCGGGGCGTCCCGCcaa                                                                          >  2:156146/1‑150 (MQ=255)
                                                                     aTTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGCCAGAGAAGCGACGATCAATCCCCGATGCAACGTTTTACGCCAGCGCGGGTCGGCACggcaaa                                                                         >  2:90747/1‑145 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGa                                                                     <  2:803598/151‑1 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAGAAAGGGGACGTTCCTCCCAGGGCGCAGCGGTTAAAGCAAACGCGGGGCGGCACGCAAAAAgg                                                                      >  1:60790/1‑150 (MQ=255)
                                                                            cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAg                                                                  <  1:1388362/150‑1 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGCCAGAAAAGGGGCGGTCATCCCCCGGCGCAGCGGTTAAAGCAAACGCGGGGGGGCCCGCaaaaaaga                                                                     >  2:1068653/1‑142 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAATCCAGAAAAGCGACGGTCCTCCCCAGGCGCAGCGTTAAAAGCAAACGCGGGGCGTCCCCCCAAATgg                                                                      >  2:1143210/1‑145 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgct                                                                <  2:192055/151‑1 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTCGGGGGGTTTAAAGCCAGATAAGGGGCGTTCATCCCCCGACGCAGCGGTTAAAACAAACCCGGGGCGGCCCGCCAAAAGGGTAggt                                                                >  2:1381985/1‑149 (MQ=255)
                                                                                 gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggc                                                             <  1:1374001/150‑1 (MQ=255)
                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                       <  2:1231591/150‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                      <  2:1396851/150‑1 (MQ=255)
                                                                                               gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggccagg                                             <  2:141951/151‑1 (MQ=255)
                                                                                                  aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTAAATGGCCGATAGGGGACGGTTCTCCCAGGACGCAGCGGTAAAAGCAAGCGCGGGGCGGAACGCAAAATGGATAGCTGGGTGgcagcaag                                                   >  1:930492/1‑141 (MQ=255)
                                                                                                  aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAggccaggcca                                          <  2:1118330/151‑1 (MQ=255)
                                                                                                    tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGCGACGGTCATCCCCGGGCGCAGCGTTAAACGCAAACGCGGGGCGTCACGCCAAATGGGTAGCTGGCTGGCAGCaagccaggcc                                           >  1:1379636/1‑148 (MQ=255)
                                                                                                     gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGaa                                       <  2:87197/151‑1 (MQ=255)
                                                                                                      ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGCTTGGCGGGGGTAAATATCCGAAAGGCGGCGTTTAACCCCCAAAGCAGCGGTAAACGCAAGCGCGGGGCGGCCAGCAAAATGGAAAACCGGGTGGGCGCacgccaggcccg                                         >  2:860823/1‑146 (MQ=255)
                                                                                                       cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGCTTGGGGGGGTTTAAAGGCAAAGAGGCGACGGTCCTCCCACGGACCACCGGTTAACGCAAACGGGGGGTGCCCCGCACAATGGGGAGCTGGGGGGGAACaggccaggcc                                           >  1:514657/1‑145 (MQ=255)
                                                                                                           gcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAATGCCAGAGAGGGGGCGTTTATCCCCGGGCGCAGCGGTAAAAGCAAGCGCGGGGCGGCACGCAAAATAGATAGCTGGCTGGGAGCAGGGCCGGCCAGAACg                                     >  2:392471/1‑147 (MQ=255)
                                                                                                           gcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGc                                 <  2:159320/151‑1 (MQ=255)
                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCa                                <  1:1118280/150‑1 (MQ=255)
                                                                                                                 cGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAATGGCACATAAGCGGCGATCATCCCAGGACGCAGCGGTAAAAGCCAACGCGGGGGGGCACGCAAAAAGGATAGCCGGCTGGGAGCAGGGCCGGCCCGAAAGTCGCATTtcc                           >  1:1146491/1‑151 (MQ=255)
                                                                                                                  ggACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTcccc                          <  2:56799/151‑1 (MQ=255)
                                                                                                                  ggACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTcccc                          <  1:592420/151‑1 (MQ=255)
                                                                                                                   gACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAATGACAGAGAGGGGACGGTCATCCCCGGACGCAGCGTTAAAAGCAAGCGGGGGGCGTGACGCCAAAAGGAGAACTGGCCGGCGGCAGGCCCGGCGCGAACGCCGCATTcccc                          >  2:286789/1‑150 (MQ=255)
                                                                                                                    aCAGCGGGATTATTTAGCGGCGGTTTGGCGGTGGTAAAAATCAGATAGGGGACGTTTATCCCAGGGCGCAGCGGTAAAAGCAAACGCGGGGCGTCACGCAAAAAGGATAACTGGCTGGCAGCAGGGCAGGGCCGAACGGCGCAATcccc                          >  2:1147035/1‑149 (MQ=255)
                                                                                                                      aGCGGGATTATTTAGCGGCGGTTTTGCGGTGGTAAATGCCAGATAGGCGACGTTTATCCCCGGACGCAGCGTTAAAAGCAAACGCGGGTCGTCACGCAAAATGGAAAGCTGGGTGGGAGCAGGGCCGGGCCAAACGTCGCATTtccc                          >  2:495145/1‑147 (MQ=255)
                                                                                                                       gCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                     <  2:322843/151‑1 (MQ=255)
                                                                                                                        cGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGAAAAGCGGCGGTTATCCCCGGGCGCAGCGTTAAACGCAAGCGCGGGGCGGCACGCAAAAAGGATAGCTGGGTGGGCGCCGGGCCGGCCCGAACGGCGCCTTtccc                          >  2:1233320/1‑145 (MQ=255)
                                                                                                                         gggATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAgg                   <  1:573218/151‑1 (MQ=255)
                                                                                                                                  tAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa           <  2:241887/150‑1 (MQ=255)
                                                                                                                                   aGCTGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc          >  2:666941/1‑150 (MQ=255)
                                                                                                                                   aGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc          <  1:332675/150‑1 (MQ=255)
                                                                                                                                         ggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc   >  1:162544/1‑151 (MQ=255)
                                                                                                                                          gTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc  >  1:293270/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCC  >  NC_000913/1356054‑1356342

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: