Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,4100TG41.8% 92.6 / 11.7 56intergenic (+112/‑35)dnaJ/insL1chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base T (4/28);  new base G (23/0);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 9.40e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT  >  NC_000913/15261‑15536
                                                                                                                                                        |                                                                                                                               
aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGtt                                                                                                                                <  1:759228/151‑1 (MQ=255)
  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTtg                                                                                                                               <  2:180421/150‑1 (MQ=255)
   gtgtGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTtgtg                                                                                                                             <  2:754998/151‑1 (MQ=255)
       gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCT‑AAACG‑CCTACTT‑TAGGGGTGTgg                                                                                                                            >  2:171165/1‑148 (MQ=255)
        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTAtt                                                                                                                        <  2:1077734/151‑1 (MQ=255)
           aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACg                                                                                                                     <  2:580354/151‑1 (MQ=255)
            agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGc                                                                                                                    <  1:490325/151‑1 (MQ=255)
             gTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGcc                                                                                                                   <  1:663325/151‑1 (MQ=255)
                  ttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGCTTACCCCCAT‑AAACG‑CCGACCT‑AAAGGGTGTGGTATTACGcc                                                                                                                   >  1:1115647/1‑146 (MQ=255)
                   tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCC‑AAGCG‑CTAAATT‑AAAGGGTGTGGGATTATGCCTggt                                                                                                               >  1:566365/1‑147 (MQ=255)
                    gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATg                                                                                                            <  2:794083/151‑1 (MQ=255)
                      cgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTAATCCCCT‑AAGCG‑CCCACCT‑TAGGGGTGGGGGAGTTCGACCGATATg                                                                                                            >  2:898420/1‑149 (MQ=255)
                        acCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAACG‑CTAACTT‑AAAGGGTGTGGGATTTCGCCTGATATGAttt                                                                                                        >  1:1118499/1‑151 (MQ=255)
                         ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                       <  1:384032/151‑1 (MQ=255)
                             aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCT‑AAGAG‑CTCAATT‑AAAGGGTGGGGGGTTTCGCCTGCTATGAatt                                                                                                        >  1:182743/1‑146 (MQ=255)
                             aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATAAGATTTAAc                                                                                                     <  2:699596/149‑1 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACTCCCCT‑AAGCG‑CGAAATT‑TAGGGGTGTGGGTTTTCGCCCGGTATTAttt                                                                                                        >  2:990812/1‑145 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                                                                                                   <  1:990839/150‑1 (MQ=255)
                                cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGca                                                                                                <  2:1000454/151‑2 (MQ=255)
                                cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGca                                                                                                <  2:1000456/151‑2 (MQ=255)
                                  cTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCC‑AAACG‑CTGACTT‑TAAGGGTGTGGGGTTATCCCTGGTATGtttttaaatgcc                                                                                                >  2:1026181/1‑141 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATAAATCCCAT‑AAACG‑CTAACTT‑AAAGGGTGTGGTATTACGCCCGAAATGATTTTAAGTGcc                                                                                                >  2:561837/1‑147 (MQ=255)
                                     aCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa                                                                                           <  1:373740/151‑1 (MQ=255)
                                      ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                          <  1:807768/151‑1 (MQ=255)
                                      ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CGAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                          <  2:654213/151‑1 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                       <  1:752302/151‑1 (MQ=255)
                                          cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc                                                                                      <  1:840106/151‑1 (MQ=255)
                                              aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCT‑AAAAG‑CTGAAAT‑AAAGGGTGTGGTGTTAAGCCTGATATGGTTTTAAACGCCGATGaa                                                                                          >  1:506212/1‑143 (MQ=255)
                                              aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCC‑AAACG‑CTCAATT‑AAGGGGTGTGGGATTACGCCTGATATTTTTTAAACTGcc                                                                                                >  1:1072006/1‑137 (MQ=255)
                                              aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                                  <  2:408385/151‑1 (MQ=255)
                                                     gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa                                                                            <  2:266980/150‑1 (MQ=255)
                                                        cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAACG‑CTCACTT‑AAGGGTTGTGGTGTTACGCCTGATATGATTTAACGTGCCGATGAATTAATCTCACGATaa                                                                           >  1:837599/1‑148 (MQ=25)
                                                         gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCAA‑AAAGC‑CTAACTT‑TAAGGGTGTGGGATTAAGCCTGCTATGATTTAACAAGTCGAAGAATTTCTCTCCCCATAACTg                                                                        >  1:1073297/1‑150 (MQ=255)
                                                           gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                                                     <  1:614548/151‑1 (MQ=255)
                                                              cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAAAG‑CTAAATT‑TAAGGGTGTGGGATTTAGCCTGGTATGATTTTAAGTGCCGATGAATTAATCTCCCGATAACTgg                                                                       >  2:463900/1‑146 (MQ=255)
                                                              cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGACTCCCCT‑AAAGG‑CTAACTTCAAGGGTT‑TGGTATTTCCCCCGCTTTGATTTTACGTGCCGGAGGATTAATCTCAccaaaaatagggcg                                                                   >  1:14890/1‑137 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCT‑AAACG‑CCAACAT‑AAAGGGTGTGGGGGTATGCCCGGTTTTATTTAACCTGCCGGAGAAATAATCTCACGATAACAGGGCa                                                                    >  1:617618/1‑147 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                                                <  1:540027/151‑1 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATAAACCCCCT‑AAACG‑CTAACTT‑TAAGGGTGTGGTGTTTCACCCGGTTTTATTTAATCTGCCGATGGATTACTCTCCCGCAAAATGGTCAGCa                                                                 >  1:565225/1‑150 (MQ=255)
                                                                     gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCCCTCCCAT‑AAAAG‑CTCAATT‑AAGGGGTGTGGGGTTTCGCCTGATTTGATTTTAAGGGCCGCTGAATTAATTTCACGaaaaataggccacaaattt                                                            >  2:579739/1‑135 (MQ=255)
                                                                     gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCC‑AAAAG‑CTAACTT‑AAAGGGTGTGGGGTTACGCCTGTTTTTATTTAACGTGCCGATGAAATACTCTCACCATAAATTGTCAGCAATTCt                                                            >  2:707816/1‑150 (MQ=255)
                                                                      ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                          <  2:672461/151‑1 (MQ=255)
                                                                          aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCT‑AAACG‑CTAACTT‑AAGGGGTGTGGGATTACGCCTGATATGATTTTACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAAATCTGGccc                                                       >  1:1174788/1‑150 (MQ=21)
                                                                            aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCT‑AAAAG‑CCCACTT‑AAGGGGTGGGGGATTACGCCCGCTATGATTTAACGTGCCGGTGAATTACTCTCACCATAACTGGTCAGGAATTCTGGccc                                                       >  2:416263/1‑148 (MQ=255)
                                                                               aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTAACTCCCAT‑AAACG‑CTAAATT‑TAGGGGTGTGGGGTTACGCCTGATATGATTTAAAGTGCCGACGGATTAATTTCACAAAAACATGGCAGACATTCTTGGCCACAtt                                                  >  2:835729/1‑150 (MQ=255)
                                                                                     cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCCCTCCCCT‑AAACG‑CCAACTT‑TAGGGGTGTGGGATTTCGCCCGTTATGATTTTAAATGCCGATGAATTACTCTCACCATAAATGGTCAGCAAATCTTGGCCCTAtt                                                  >  1:677184/1‑144 (MQ=255)
                                                                                        tGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGAGACTCCCCT‑AAGGG‑GCAACTT‑AAGGGGTGTGGGGTTACGCCTGATATGATTTTACATTGCGGTGAAATACTCTCCAGAGAAAATGTCAGCAAATCTTGGCCACATTGGT‑AAGc                                          >  2:913716/1‑148 (MQ=255)
                                                                                            gATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCgaaga                                   <  1:480022/151‑1 (MQ=255)
                                                                                             atatGCGAGCACCTGTAAAGTGGCGGGGGTCCCTCCCCC‑AAACGCCTACCTA‑AAGGGTT‑TGGGATTTCTCCCGGTATGATTTTAAGTTCCGATGAATTTCTCTCCCCATAAATGGTCCGCAATTCTTGCCCCCATTTGT‑AAGccc                                        >  2:621493/1‑145 (MQ=255)
                                                                                               atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACt                                <  1:368365/151‑1 (MQ=255)
                                                                                                 gCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTgg                              <  1:407607/151‑1 (MQ=255)
                                                                                                     gCACCTGTAAAGTGGCGGGGGTCACCCCCCT‑AAACG‑CCAAAAT‑AAGGGGTGGGGGGTTAAGCCTGATATGATTTAACGTGCCGCTGGAATAATCTCACGATAACTGGTTCGCAATTTTGGCCCCCATTTGT‑AAAAACGAAGAAAAGGATa                           >  2:627366/1‑150 (MQ=9)
                                                                                                       aCCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACtt                        <  2:1161389/151‑1 (MQ=255)
                                                                                                         cTGTAAAGTGGCGGGGGTCACCCCCCTAAACCC‑CCAATTT‑AAGGG‑TGTGGGATTACGCCTGATATGATTTAAAGTGCCCCTGAATTAATCTTCCGATAAATGGGCAGCCAATTTGGGCCACATTTGTAAAACCCGAAAAACgggtaatt                         >  2:386378/1‑144 (MQ=2)
                                                                                                                   gCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATg            <  2:138699/151‑1 (MQ=34)
                                                                                                                         aTCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgggg      <  2:206196/151‑1 (MQ=33)
                                                                                                                             cTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctct  >  1:1064235/1‑151 (MQ=11)
                                                                                                                                                        |                                                                                                                               
ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT  >  NC_000913/15261‑15536

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑