Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG40.9% 108.4 / 19.1 66T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (7/32);  new base G (27/0);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 2.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356054‑1356351
                                                                                                                                                      |                                                                                                                                                          
cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGtt                                                                                                                                                         >  1:90474/1‑151 (MQ=255)
    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTaaa                                                                                                                                                      <  1:924316/150‑1 (MQ=255)
         aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCa                                                                                                                                                 <  2:19833/150‑1 (MQ=255)
            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATa                                                                                                                                             <  1:293473/151‑1 (MQ=255)
             ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAg                                                                                                                                            <  1:282973/151‑1 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTtaaagac                                                                                                                                                  >  2:533783/1‑140 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTtaaattcaaaaagg                                                                                                                                           >  1:887087/1‑140 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGCACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTTCGGTGTTTAATACCAGaaa                                                                                                                                             >  1:1357983/1‑147 (MQ=255)
                 gAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAGAAAGGCg                                                                                                                                         >  1:856402/1‑150 (MQ=255)
                  aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAATGCCAGATAAGCg                                                                                                                                         >  1:614214/1‑149 (MQ=255)
                     cgcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGtt                                                                                                                                   <  1:352833/151‑1 (MQ=255)
                       cgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTc                                                                                                                                  <  1:1278589/150‑1 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATccc                                                                                                                             <  1:475139/151‑1 (MQ=255)
                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAATACCCGAGAGGCGACGGTTATCCCAgg                                                                                                                          >  1:437205/1‑151 (MQ=255)
                                   gagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCa                                                                                                                     <  1:1311413/151‑1 (MQ=255)
                                     ggCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAg                                                                                                                   <  2:846703/150‑1 (MQ=255)
                                           cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAAGGCAAAAAGGGGACGTTCATTCCCAGAAGCAGCggt                                                                                                              >  2:961707/1‑146 (MQ=255)
                                              ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAAAGCCCGAAAGGCGACGTTTCTACCCCGACGCCGCGGTaaa                                                                                                           >  2:597092/1‑148 (MQ=255)
                                               tCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGc                                                                                                        <  2:965943/150‑1 (MQ=255)
                                                cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCaa                                                                                                      <  1:114117/151‑1 (MQ=255)
                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATTCCAGAAAAGCGACGGTCATCCCCGGACGCAGCGTTaaa                                                                                                           >  1:152095/1‑145 (MQ=255)
                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGGCAGCGAGGCGGCGTTCATTCCCGGACGCAGCGGTaaaaacaa                                                                                                      >  2:515611/1‑145 (MQ=255)
                                                  gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTAAAAGACAAAAAGGCGGCGTTTATCCCCAGGCGCAGCggttaa                                                                                                           >  1:889/1‑139 (MQ=255)
                                                   ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgcg                                                                                                   <  1:79163/151‑1 (MQ=255)
                                                    cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAATGCCCAATAAGCGACGTTTATCCCCCGACGCAGCGGTAAACACCAgcg                                                                                                   >  2:526333/1‑150 (MQ=255)
                                                       aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAATGCCCGAAAGGGGGCGTTCATCCCCGGAAGCAGCGGTTAAAACAAGAGCggg                                                                                               >  2:797407/1‑151 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                              <  2:513170/151‑1 (MQ=255)
                                                         cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAATGCCCGAAAAGCGGAGGTCTTTCCCCGGCGCAGCGGTTAAAGCAAACCCggg                                                                                               >  2:1320040/1‑149 (MQ=255)
                                                         cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtc                                                                                             <  1:1666/151‑1 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAATGGCCGAGAAGGGACGGTTCTCCCACGACGCAGCGGTTAAACCAAGCGCgggggcgccc                                                                                        >  2:1256913/1‑144 (MQ=255)
                                                               aGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTTAATTCCAGAAAAGGGACGGTCATCCCCCGACGCCGGGGTTAACGCAAGCCCGGGGCGGCAc                                                                                        >  2:143156/1‑150 (MQ=255)
                                                                 aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAAGCCCAATAAGCGACGGTCATCCCAAGGCGCAGCGGTTAAAACAAGCGCgggggg                                                                                            >  1:166677/1‑141 (MQ=255)
                                                                  tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATATCAGATAAGGGGCGGTCATCCCCCGACGCAGCGGTAAAAGCAAGAGCGGggcg                                                                                            >  1:462107/1‑143 (MQ=255)
                                                                   ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAAATGCAGAAAAGGGACGGTTCATCCCCGGAGCACCGGTAAAAGCCAGCGGGGGGCGGCACGc                                                                                      >  2:674654/1‑148 (MQ=255)
                                                                      ttGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAATTACAAAAAAGGGGCGGTTCTTCCCGGACGCAGCGGTAAACGCAAGGGCGGGGGGTGACGCaaaa                                                                                  >  1:440040/1‑149 (MQ=255)
                                                                       tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAAAGGCACATATGAGAAGGTCATCCCAGGACGCAACGTTAAAAGCAAGAGCGGGGCGGCCCCCaaaaa                                                                                 >  1:1287508/1‑148 (MQ=255)
                                                                       tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATg                                                                                <  2:440519/150‑1 (MQ=255)
                                                                       tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATg                                                                                <  2:10931/150‑1 (MQ=255)
                                                                          aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATa                                                                            <  2:1382476/151‑1 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGTTAAATGCCCGATAAGCGGCGTTCATCCCAGGGCGCAGCGGTAAAAGCAAACGCGGGGCGTCACGCCAAAttggta                                                                            >  1:942697/1‑143 (MQ=255)
                                                                                 gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAAAACCAGAGAAGAGAGGGTCCTCCCCGGGAGCAGGGGTTAACAGCAACGCGGGGCGGCGCGCaaaaaaggaaacggg                                                                       >  2:986424/1‑138 (MQ=255)
                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                                <  1:854176/150‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagca                                                             <  1:475378/151‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                               <  2:1204069/150‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                               <  2:737197/150‑1 (MQ=255)
                                                                                          taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATAACAGATAAGGGACGGTTTTCCCCGGGCGCAGCGGTAAACCCAAACGCGGGGCGGCCCGCCAAATGGATAGCTGGCTGGCAGcagg                                                           >  1:205209/1‑151 (MQ=255)
                                                                                            aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcc                                                         <  1:797433/151‑1 (MQ=255)
                                                                                              cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcca                                                        <  1:866756/150‑1 (MQ=255)
                                                                                                ccAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCATCACGCAAAATGGATAGCTGGCTGGCAGcaggccaggc                                                     <  2:1150777/151‑1 (MQ=255)
                                                                                                  aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAAAGCCCGATAGGCGACGTTCATCCCCCGGCGCAGCGGTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGATAGGTGGGCGGCAGCaggccaggcc                                                    >  1:727894/1‑150 (MQ=255)
                                                                                                   gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAggccaggcca                                                   <  1:815502/150‑1 (MQ=255)
                                                                                                     gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGCCAGATAAGCGGCGTTCATCCCCAGGCGCAGCGTTTAACACAAACGCGGGTCGTCACGCCAAAAGGATAGGTTGCTGGCAGCaggccaggcc                                                    >  1:886858/1‑147 (MQ=255)
                                                                                                          cgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCg                                           <  2:1204668/151‑1 (MQ=255)
                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCAGATAGGCGACGTTCATCCCCGGACCCAGCGGTTAAAGCCAACGCGGGGCGTCACGCAAAATTGGTAACTGGCTGGCAGCAGGGCAGGCCCGAACGTCGc                                          >  1:792294/1‑149 (MQ=255)
                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCAt                                        <  2:833747/151‑1 (MQ=255)
                                                                                                              tGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGGGGGGGTAAAATCCAGAAAGGCGAGGTTCATCCCCAGACGGAGCGTTTAACGCAAACGCGGGTCGTCACGCCAAATAGGGAACTGGCTGGCGGGAGGACCGGGCAGaaagccg                                           >  1:381563/1‑142 (MQ=255)
                                                                                                               ggCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAAGGCAGGTAGGCGACGTTTATTCCAAAAGGCAGCCATAAAAAAAAGCGGGGGGCGTGAAGCAAAAAAGATAAATAGCTGGCCGCAaggccgggcataaacacgaat                                         >  2:1083222/1‑137 (MQ=255)
                                                                                                                  ggACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTcccc                                   <  2:887058/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:953810/151‑1 (MQ=255)
                                                                                                                               attTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                       <  1:1259577/150‑1 (MQ=255)
                                                                                                                                tttAGCTGCGGTTTGGCGGTGTTAAAGGCCAGATAGGCGACGTTCAGCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTc                      <  2:468761/150‑1 (MQ=255)
                                                                                                                                     cggcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa                <  2:1207807/151‑1 (MQ=255)
                                                                                                                                        cggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCCGGCCAGAACGGCGCATTCCCCGGTCAGGGGTTTCGACAg                 >  2:219990/1‑147 (MQ=255)
                                                                                                                                        cggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg              >  2:422652/1‑150 (MQ=255)
                                                                                                                                          gTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc           >  1:1231218/1‑151 (MQ=255)
                                                                                                                                             ttGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg         >  1:1233407/1‑151 (MQ=255)
                                                                                                                                                     gtgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc   >  2:1191789/1‑151 (MQ=255)
                                                                                                                                                      tgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGGCAGGAGTTTCGACAGACGTCCGGTgccgcc  >  1:568572/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                          
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356054‑1356351

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: