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

ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15261‑15534
                                                                                                                                                        |                                                                                                                             
aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGtt                                                                                                                              <  2:826349/151‑1 (MQ=255)
  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTtgt                                                                                                                            <  4:433968/151‑1 (MQ=255)
  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTtg                                                                                                                             <  2:227565/150‑1 (MQ=255)
    tgtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑GCACCCCCCCAAAA‑GCT‑TAATTAAAGGGTGTgg                                                                                                                          >  1:1393641/1‑151 (MQ=255)
     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TCACTCCCCCAAAC‑GCC‑CACTTAAGGGGTGTgg                                                                                                                          >  3:7125/1‑150 (MQ=255)
     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTgg                                                                                                                          <  2:740992/150‑1 (MQ=255)
     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGt                                                                                                                         <  2:748916/151‑1 (MQ=255)
     ctgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TGACCCCCCTGAAC‑GCC‑GACTTTAGGGGTgggg                                                                                                                          >  1:59586/2‑150 (MQ=255)
       gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑t                                                                                                                      <  2:310808/151‑1 (MQ=255)
        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑tt                                                                                                                     <  1:184213/151‑1 (MQ=255)
        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑GTACTCCCCCAAAC‑GCC‑AACTTAAGGGGtgtg                                                                                                                           >  2:851479/1‑146 (MQ=255)
         agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TCACCCCCCTAAGC‑GCG‑AACTTTAGGGGTGTGGtt‑ttt                                                                                                                    >  1:790372/1‑150 (MQ=255)
          gaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTAc                                                                                                                   <  3:536973/151‑1 (MQ=255)
            agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGc                                                                                                                 <  1:772287/151‑1 (MQ=255)
                ttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑GTACTCCCCTAAAA‑GCC‑AACTTTATGGGTGGGGGA‑ttttg                                                                                                                  >  4:54815/1‑143 (MQ=255)
                 tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGa                                                                                                             <  1:75323/150‑1 (MQ=255)
                 tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑AAACTCCCCCAAGC‑GGT‑AACATAAAGGGTGTGGTA‑TTACGcccggt                                                                                                            >  2:1131683/1‑147 (MQ=255)
                   tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGata                                                                                                           <  4:186590/150‑1 (MQ=255)
                    gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATg                                                                                                         <  2:837355/151‑1 (MQ=255)
                       gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TCACTCCCCCAAAC‑GCT‑AACTTTAGGGGTGTGGGA‑TTTAGCCCGATattatt                                                                                                      >  1:1417399/1‑151 (MQ=255)
                        acCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑GCACCCCCCTAAAC‑GCT‑AACTTTAAGGGTGTGGGA‑TTAAGCCTGGTATGAttt                                                                                                     >  2:1419102/1‑151 (MQ=255)
                         ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑TCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTa                                                                                                    <  1:1174638/151‑1 (MQ=255)
                          cTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑GCACCCCCCCAAGA‑GCG‑AACTTAAGGGGTGTGGGA‑TTTTGCCTGATATGAttt                                                                                                     >  1:543263/1‑149 (MQ=255)
                            gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACg                                                                                                 <  3:108990/151‑1 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTg                                                                                               <  2:1089350/151‑1 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTg                                                                                               <  1:796463/151‑1 (MQ=255)
                                 gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCg                                                                                            <  3:54816/151‑1 (MQ=255)
                                   tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATTGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TCACTCCCCTAAGC‑GCT‑AACTTAAGGGGTGTGGTA‑TTAAGCCTGATATGATTTAACGTGCCg                                                                                            >  2:862164/1‑149 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGTGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCATGATATGATTTAACGTGCCGAt                                                                                          <  1:1370731/150‑1 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑GTACTCCCATAAGC‑GCT‑AACTTAAGGGGTGTGGGA‑TTTAGCCTGATATTAttttaaaag                                                                                               >  1:394026/1‑139 (MQ=37)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGa                                                                                           <  1:1358460/149‑1 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGAt                                                                                          <  1:139707/150‑1 (MQ=255)
                                     aCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGa                                                                                        <  2:18388/151‑1 (MQ=255)
                                       cTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TCACTCCCATAAAC‑GCT‑AACTTAAAGGGTGTGGGA‑TTTAACCCGATATGATTTAAAGTGCCGATGaa                                                                                       >  2:243166/1‑150 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAAtt                                                                                     <  3:341845/150‑1 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAAtt                                                                                     <  1:1353172/150‑1 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                    <  1:1285692/151‑1 (MQ=255)
                                          cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                      <  2:803789/148‑1 (MQ=255)
                                                aGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                                             <  2:481734/151‑1 (MQ=255)
                                                 gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACg                                                                            <  2:1417466/151‑1 (MQ=255)
                                                 gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACg                                                                            <  1:924714/151‑1 (MQ=255)
                                                   cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TGTCTCCCCATAGG‑GGT‑TACTTAAGGGGTGTGGGG‑TTAAGCCTGGTATGATTTTAAGCGCCGACGGAATATTCCCACGAt                                                                          >  1:557347/1‑151 (MQ=255)
                                                           gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                                                  <  1:1455863/151‑1 (MQ=255)
                                                            agCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGAGA‑TTAATCCCCCAAGA‑GGTAAACATAAAGG‑TGCGGGG‑TTACGCGTGTAAATATTTAAAAAGTCGATGGATTACACTCACGCTaaactgca                                                                  >  4:106284/2‑144 (MQ=255)
                                                               aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TGACTGCCATAAAA‑GCG‑AACTTAAGGGGTGTGGGG‑TTACGCCCGGGTTGGTTTTAAGTGCCGATGAATTACTCTTACGGTAAATGGTCAGc                                                               >  2:404734/1‑150 (MQ=255)
                                                               aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑GGAACCCCATAAGG‑GCC‑AACCTAAGGGGTGGGGGATTTACGCC‑GGTATGGTTTAAACTGGCGGAGAAAATCTTCCCCGATAACTGGGCAGc                                                               >  2:669239/1‑150 (MQ=255)
                                                               aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑ACACTCCCATAAAC‑GCC‑CAATTAAGGGGTGTGGTG‑TTAAACCTGATAAGATTTTACGGGGCGATGAAATAATCTCACggaaaattgtcgg                                                                >  2:948458/1‑137 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACCCCCCTAAAC‑GCT‑AACTTAAGGGGTGTGGGA‑TTAAGCCTGGTATGATTTAAAGTGCCGATGAATTAATCTCACGATAAATgg                                                                    >  3:193646/1‑144 (MQ=255)
                                                                     gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTc                                                          <  2:334826/149‑1 (MQ=255)
                                                                      ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                       <  2:1543421/151‑1 (MQ=255)
                                                                      ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                       <  1:409683/151‑1 (MQ=255)
                                                                          aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TGACTCCCATAAAC‑GCT‑AACTTTAAGGGTGGGGGA‑TTACGCCTGATATGATTTAAAGTGGCGATGAATTATTCTCACGATAAATGGTTAGGAATTCTGGcccc                                                   >  2:924768/1‑150 (MQ=255)
                                                                           aaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACATAAAGGGTGTGGTA‑TTACGCCTGGTATTATTTAAAGTGCCGATGAAATACTCTCACGATAAATGGTCAGCAATTTTGGccc                                                    >  4:243918/1‑149 (MQ=255)
                                                                             aTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑GCACTCCCATAAAC‑GCT‑AACTTAAGGGGTGTGGGG‑TTAAGCCCGGTATTATTTTAACTGCCGATGAATTAATTTCACGATAAATGGGCAGGAATTCttggcc                                                    >  2:1256521/1‑142 (MQ=255)
                                                                               aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                              <  2:135452/151‑1 (MQ=255)
                                                                                tggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑TCACTCCCATAAGCGGCA‑AACTT‑AGGGGTGTGGGA‑TTAAGCCCGATATGATTTAAAGTGCCGATGAATTACTTTCACGATAAATGGTTAGCAATTTTGGCCCCTAtt                                               >  4:76461/2‑149 (MQ=255)
                                                                                gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg                                             <  1:433312/151‑1 (MQ=255)
                                                                                 ggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt                                            <  1:1694154/151‑1 (MQ=255)
                                                                                  gTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGAGATTCCCTCCCA‑AAAC‑GCC‑AAATTTACAGGTGTGGGA‑TTTCGTACGGTATGATTTAACGTGACGGCGAAATAATCTTCCGATAAATTGTAAGCAATTCTgtaccctcttgg                                             >  1:757748/1‑138 (MQ=255)
                                                                                   tGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑ACCATCCCCTAAAC‑GCT‑AAATTAAAGGGTTTGGGA‑TTTACCCTGGTATTATTTTAAGTGCCGGTGAATTACTCTCACCCTAAATGGGCCGCAAATCTAGCCCCTATTTGt                                            >  1:1589066/1‑149 (MQ=255)
                                                                                     cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGc                                        <  2:184134/151‑1 (MQ=255)
                                                                                         gAAGATATGCGAGCACCTGTAAAGTGGCGGGGG‑ACACCCCCCTAAAA‑GCC‑AACTTAAGGGGTGTGGTA‑TTACGCCTGATATTATTTAAAGTGCCGATGAATTACTCTCACGAAAAATAGTCAGCAAATTTTGCCCATTTTTGTAAGacc                                      >  2:1543482/1‑149 (MQ=17)
                                                                                         gAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCCTAAAC‑GCT‑AACTTAAAGGGTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAATTTTGGCCCCTATTTGTAAGccc                                      >  2:365086/1‑149 (MQ=17)
                                                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCACAAAC‑GCT‑AACTTAAAGGGTGTGGGG‑TTTCGCCTGATATGATTTAACGTGCCGCCGAAATAATCTCCCCATAAATGGTCAGCAATTCTTGCCCCTATTTGTAAAccc                                      >  4:593940/1‑148 (MQ=17)
                                                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCAATCCCCTAAGC‑CCC‑AAATTTAGGGGTGTGGGA‑TTTTGCCTGATAAGATTTAAAGTGCCGGTGAATTAATCTTACGATAAATAGTCCGCAATTCTTGGCCATATTGGTaaagcc                                      >  3:530016/1‑144 (MQ=255)
                                                                                           aGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTTCCCTAAGC‑ACG‑CACATAAGGGGGGTGGGA‑TTAAGGCTGTTATGATTTTACAAGCCGGTGAATTACTCTCCCTATAAAAGGTGAGCAATTTTTTCCCCTAATTGTAAGacc                                      >  3:121798/1‑147 (MQ=255)
                                                                                           aGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaag                                  <  4:334766/151‑1 (MQ=255)
                                                                                            gATATGCGAGCACCTGTAAAGTGGCGGGGA‑GTACACCCCCAAGC‑GCC‑TAATTCAGGGGTGGGGGG‑TTATCCCTTGTATGAGTTAACGCGCCGAAGAATTTATCTCACTACAACACGTCTGCAAGTCTGGTCCCTATTgcgtagtgcacca                                  >  3:581843/1‑138 (MQ=255)
                                                                                             atatGCGAGCACCTGTAAAGTGGCGGGGG‑TTACTCCCCTAAAC‑GCC‑AACTTAAAGGGTGTGGGG‑TTTCCCCTGATATTATTTTAAGTGCCGGTGAAATAATCTCACGATAACTTGTCAGCAAATCTTGCCCCTATTTGGAAGACCgaa                                   >  2:1686576/1‑148 (MQ=255)
                                                                                             atatGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaagaa                                <  1:964629/151‑1 (MQ=255)
                                                                                               atGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAAc                               <  4:384172/150‑1 (MQ=255)
                                                                                               atGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                              <  1:1561499/151‑1 (MQ=255)
                                                                                                tGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCCTAAGC‑GCT‑TAATTAAGGGGTGTGGGA‑TTATGCCTGGTATTATTTTACATGCCGATGAATTACTTTCACGATAAATTGTCAGCAATTCTTGCCCATATTTGTAAGCCCggagaa                                >  2:1452242/1‑148 (MQ=255)
                                                                                                tGCGAGCACCTGTAAAGTGGCGGGGA‑TCACCCCCCTAAGC‑GCT‑AACTTTAAGGGTGGGGGG‑TTAAGCCTGATATTAATTTAAGTGCCGGTGAATTAATCTCACGATAAATTGGCAGCAAATTTGGGCCCTATTTGTAAAACCgaagaa                                >  4:451512/1‑148 (MQ=255)
                                                                                                 gCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCCTAAAC‑GCT‑AACTTAAAGGGTGTGGGA‑TTACGCCTGATATGATTTTACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAATTCTTGCCCCTATTTGTAAGCCCGAAGAACTgg                            >  2:1564771/1‑151 (MQ=17)
                                                                                                 gCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                            <  2:170319/151‑1 (MQ=255)
                                                                                                     gCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATAc                        <  1:1609414/151‑1 (MQ=255)
                                                                                                        ccTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACtt                      <  2:459354/150‑1 (MQ=255)
                                                                                                        ccTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc                     <  3:209243/151‑1 (MQ=255)
                                                                                                             aaaGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc                <  2:1296883/151‑1 (MQ=39)
                                                                                                                gTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa             <  1:1527185/151‑1 (MQ=37)
                                                                                                                gTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa             <  2:1273469/151‑1 (MQ=37)
                                                                                                                   gCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg          <  2:1208264/151‑1 (MQ=34)
                                                                                                                   gCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg          <  4:517447/151‑1 (MQ=34)
                                                                                                                    cGGTGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTAAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGc         >  1:66559/1‑151 (MQ=18)
                                                                                                                      gggA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg       <  4:521751/151‑1 (MQ=33)
                                                                                                                       ggA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgg      >  1:1453868/1‑151 (MQ=33)
                                                                                                                           tCACTCCCATAAGC‑GCT‑AACTTAAGGGGTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgggg    >  1:908042/1‑150 (MQ=17)
                                                                                                                            cACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct  <  1:243126/151‑1 (MQ=31)
                                                                                                                            cACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct  >  2:1105355/1‑151 (MQ=31)
                                                                                                                                                        |                                                                                                                             
ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGA‑TCACTCCCATAAGC‑GCT‑AACTTAAGGGTTGTGGTA‑TTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15261‑15534

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: ‑