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

GTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15266‑15534
                                                                                                                                                   |                                                                                                                               
gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTgg                                                                                                                            <  1:66979/150‑1 (MQ=255)
gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGt                                                                                                                           <  2:133978/151‑1 (MQ=255)
 tgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGA‑CTCCCAAAAGCG‑GC‑AAATTTAGGGGTGTgg                                                                                                                            >  1:14511/1‑149 (MQ=255)
      aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACg                                                                                                                     <  2:39966/151‑1 (MQ=255)
      aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACg                                                                                                                     <  1:182841/151‑1 (MQ=255)
       agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑c                                                                                                                   <  1:570257/151‑1 (MQ=255)
           ttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACA‑CACCCCAAAGAA‑CT‑CACTTATGGGGTGTGGGATTatg                                                                                                                     >  2:279307/1‑144 (MQ=255)
            tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCCAAGCG‑CC‑CACTTAAAGGGTTTGGGGTTATG‑CCTgc                                                                                                               >  1:922967/1‑149 (MQ=255)
            tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGa                                                                                                               <  2:92419/150‑1 (MQ=255)
              tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CCCCCCTAAGCG‑CT‑CAATTTAGGGGTggggggttttgcctggt                                                                                                               >  1:562993/1‑136 (MQ=255)
                  gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTCCCATAAAAG‑CC‑AACTTTAGGGGTGTGGGGTTATG‑CCTGCTATg                                                                                                           >  1:680411/1‑148 (MQ=255)
                  gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCCTTAAAG‑CC‑AACTTTAGGGGTGTGGGATTAAG‑CCCGATATGG‑tt                                                                                                       >  2:799985/1‑151 (MQ=255)
                    ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CT‑AACTTTAAGGGTGTGGGATTACG‑CCTGATATTA‑ttt                                                                                                      >  2:70093/1‑150 (MQ=255)
                      tGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑ATCCCCCAAACG‑CC‑AAATTTAGGGGTGTGGGATTTCG‑CCCGATATTA‑ttt                                                                                                      >  2:50998/1‑148 (MQ=255)
                      tGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTCCCCTAAAAG‑GT‑AAATTTAAGGGTGTGGGATTTTG‑GCTGATATGA‑ttt                                                                                                      >  2:85461/1‑148 (MQ=255)
                       gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACg                                                                                                  <  1:90447/151‑1 (MQ=255)
                        aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGt                                                                                                 <  1:161803/151‑1 (MQ=255)
                         cccGCTAACCTTCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAACG‑CT‑AACTTAAAGGGTGTGGGGTTTCG‑CCTGATATTA‑ttt                                                                                                      >  2:659679/1‑145 (MQ=255)
                          ccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CCCCCCTAAACG‑CT‑AACTTAAAGGGTGTGGGATTTAG‑CCCGATATGA‑TTTAAAGt                                                                                                 >  1:48197/1‑149 (MQ=255)
                          ccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGc                                                                                               <  2:351006/151‑1 (MQ=255)
                           cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCCTTAAGGC‑CA‑AATTTAAGGG‑TGGGGGATTACG‑cccggctttgttt                                                                                                       >  1:182816/1‑132 (MQ=255)
                            gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTA‑CTCCCCTAAACG‑CC‑AACCTAAAGGGTGTGGGTTTATG‑CCCCGTATGA‑ttt                                                                                                      >  1:970161/1‑142 (MQ=255)
                             cTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTCCCCTAAACG‑CC‑AACTTTAGGGGTGTGGGATTACG‑CCTGATATGA‑TTTAAAGt                                                                                                 >  2:237602/1‑146 (MQ=255)
                                 ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGA‑CCCCCCTAAGCG‑CC‑AACTTAAAGGGTGTGGGATTAAG‑CCTGATATGA‑ttttaagggcc                                                                                              >  1:467820/1‑137 (MQ=255)
                                        aaaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CC‑AACTTTAAGGGTGTGGGATTATG‑CCTGGTATTA‑TTTAAAGTGCCCATGAAtt                                                                                      >  2:413812/1‑146 (MQ=255)
                                          aaGCCTGTCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑AACCCCTTAACGCCC‑AACTAAAAGG‑TGTGGTATTTCG‑CCTGATATTT‑TTTaatcgggcggtg                                                                                          >  1:682245/1‑130 (MQ=255)
                                            gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CT‑AACTTAAGGGGTGTGGGATTACG‑CCTGATATGA‑TTTTAAATGCCGATGAAtt                                                                                      >  2:894173/1‑142 (MQ=255)
                                                gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACA‑CTCCCATAAGCG‑CT‑CACTTAAAGGTTGTGGGATTTTG‑CCCGATATTA‑TTTTAAGTGCCGATGAATTTCTCTCCCGATa                                                                          >  1:273020/1‑150 (MQ=255)
                                                  ccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCATCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATa                                                                          <  1:813608/148‑1 (MQ=255)
                                                   cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                      <  1:80035/151‑1 (MQ=255)
                                                   cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                      <  2:285737/151‑1 (MQ=255)
                                                    gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTgg                                                                     <  2:221084/151‑1 (MQ=255)
                                                    gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCA‑CTCCCATAAGCG‑CT‑AACCTAAGGGGTGTGGGGTTAAG‑CCTGATATTA‑TTTTACGTGCCGCTGAAATTTTCTCCCCATAAATgg                                                                     >  1:773469/1‑151 (MQ=255)
                                                         cAGGCCTGGGTAAAAATAGGGGGCGTTGAAGATATGCGTGCAGCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                                 <  2:889876/150‑1 (MQ=255)
                                                            gCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATTA‑CTCCCCTAAGAG‑GT‑AACTTTAAGGTTGTGGTTTTACGACCTG‑TATTA‑GTTAATGTGCCGCCGAATTAATCCCATGCTAACTGGGCAgga                                                               >  1:894026/1‑147 (MQ=255)
                                                             ccTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAACG‑CT‑AACTTTAGGGGTGTGGGATTTCG‑CCTGGTATTA‑TTTTAAGTGCCGGTGAAATACTCTCACGATAAATGGGCAGCaaat                                                            >  2:991667/1‑149 (MQ=255)
                                                               tGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTA‑ATCCCCTAAAAG‑GTAAAATTAAGGG‑TGTGGGGTTTCG‑CCTGATATGG‑TTTTAAGTGCCGGTGAATTACTCTCTCGATAACTGGTCAGCAAtt                                                            >  1:951817/1‑149 (MQ=255)
                                                                gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt                                                          <  1:498502/150‑1 (MQ=255)
                                                                  gTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGA‑CCCCCACAAGAG‑CT‑TACTTAAGAGGTGTGGGATTAGG‑CCTGTTATGATTTTAAAGT‑CGGATGAATTATTCTCTCCATTACttagtagccattcttg                                                        >  1:648249/1‑135 (MQ=255)
                                                                   tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CC‑AACTTAAGGGGTGTGGGATTACG‑CCTGGTATGA‑TTTTACGTGGCGATGAAATACTCTCACGATAACTGGGCAGCAATTTTTGcc                                                      >  2:73462/1‑151 (MQ=255)
                                                                       aaTAGGGGGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑ATCCCCTAAAAG‑CG‑AAATTTAGGGGTTGGGGATTACC‑GCCGATATGA‑TTTTAAGTGCCCATGAATTAATCTCACCGAAAATTGTTTGCAATTCTTGccc                                                     >  2:272594/1‑148 (MQ=255)
                                                                         tAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTA‑CTCCCATAAACG‑CT‑AAATTTAGGGGTGGGGGATTAAG‑CCCGATATTA‑TTTTAAGTGCCGATGAAATTATTTCCCGATAACTTGGTAGCAATTTTTGCCCCTA‑tt                                               >  1:952474/1‑151 (MQ=255)
                                                                           gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTgg                                             <  2:150152/151‑1 (MQ=255)
                                                                             gTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTa                                           <  2:201183/151‑1 (MQ=255)
                                                                                    gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCg                                     <  1:728554/150‑1 (MQ=255)
                                                                                      aGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCgaag                                  <  2:567626/151‑1 (MQ=255)
                                                                                      aGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCgaag                                  <  1:854663/151‑1 (MQ=255)
                                                                                        atatGCGAGCACCTGTAAAGTGGCGGGGGGGA‑CTCCCCTAAGCG‑CT‑TAATCAATGGGTTTGGTATTTAG‑CCCTATATGA‑TTTTACGTGCCGCTGAATTATTTTTACGATAACTTGTTAGCAATTTCTGCCCATATTTGGTAAcccca                                     >  1:128497/1‑146 (MQ=255)
                                                                                        atatGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CCCCCCTAAACG‑CT‑AAATTAAAGGGTGTGGGATTTCG‑CCCGATATTA‑TTTTACGTGCCGATGAAATTCTCTCACCCTAACTGGTCAGCAATTCTTGCCCCTA‑TTTGTAAAccc                                      >  2:936608/1‑145 (MQ=255)
                                                                                           tGCGAGCACCTGTAAAGTGGCGGGGGGCA‑ATCCCCTCAACG‑CC‑AAATTTAGGGGTGTGGTATTACG‑CCTGGTATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACCGGGCCGCAATTCTGGCCCCTA‑TTGGTAAGCCCcaagaa                                >  1:896565/1‑148 (MQ=14)
                                                                                            gCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTg                             <  1:409357/150‑1 (MQ=255)
                                                                                                   ccTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTc                     <  1:365679/151‑1 (MQ=255)
                                                                                                      gTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCg                    <  2:273081/149‑1 (MQ=40)
                                                                                                            tGGCGGGGATCA‑CTCCCATAAGCG‑CG‑AACTTAAGGGGTGTGGTATGACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCACACGAAAACTGGTCAGCAATTCTGGCCCATC‑CTGGTAAGGCCGAAGAACTGGATAATTCGGccagta             >  1:432725/1‑145 (MQ=11)
                                                                                                             ggCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt           >  2:48349/1‑151 (MQ=35)
                                                                                                                 gggATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg       <  2:469817/151‑1 (MQ=33)
                                                                                                                 gggATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg       <  1:263477/151‑1 (MQ=33)
                                                                                                                      cA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct  <  2:123099/151‑1 (MQ=31)
                                                                                                                                                   |                                                                                                                               
GTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15266‑15534

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