Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,487,1630CT41.5% 64.1 / 24.1 53intergenic (+200/‑72)hrpA/ydcFputative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (8/23);  new base T (22/0);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.43e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCA  >  NC_000913/1487016‑1487311
                                                                                                                                                    |                                                                                                                                                          
tttttATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTtcct                                                                                                                                                          <  1:345941/150‑1 (MQ=255)
 ttttATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTtcctc                                                                                                                                                         <  2:329440/150‑1 (MQ=255)
   ttattaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGATATAAGATTtcctcc                                                                                                                                                        >  1:249505/1‑149 (MQ=255)
            cTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAATAAGATTTTCTCCtttttcaggg                                                                                                                                              >  1:95200/1‑143 (MQ=255)
            cTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAGATTTTCTCCTTTTACAGTg                                                                                                                                              >  1:637805/1‑150 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAGGAAAAAACATTTTCTCCTTTTAGAGGGCt                                                                                                                                            >  1:665641/1‑150 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGc                                                                                                                                             <  2:635014/149‑1 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCt                                                                                                                                            <  1:11585/150‑1 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGACAAACGAATTTCTCCTTTTTGAGTg                                                                                                                                              >  1:429098/1‑148 (MQ=255)
               aaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTTTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTcc                                                                                                                                          <  1:27708/151‑1 (MQ=255)
               aaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAATACAATTTTtttctttcttggggcgcc                                                                                                                                          >  2:370526/1‑141 (MQ=255)
                 aCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGAAAAACAATTttctcc                                                                                                                                                        >  2:762224/1‑135 (MQ=255)
                      aaGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTtct                                                                                                                                    <  1:92672/150‑1 (MQ=255)
                         ggTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑ATTAGGAAATACAATTTTCTTCTTTCTGAGGGCCCCGcgtttccc                                                                                                                                 >  2:512076/1‑143 (MQ=255)
                          gTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGACAAAACAATTCCTCCTTTTAAAGTGCTccccct                                                                                                                                      >  1:118419/1‑140 (MQ=255)
                                aTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTAtct                                                                                                                         <  2:511636/151‑1 (MQ=255)
                                aTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACAAAAGAATTTCTCCTTTTTGAGGGCTCCGCTTCTCCCtattat                                                                                                                           >  2:829528/1‑149 (MQ=255)
                                 ttCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAAGAAATACAATTTTCTCCTTTCTGAGGTGTCCGCCTCTCACtatt                                                                                                                             >  2:698646/1‑146 (MQ=255)
                                   cAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCAc                                                                                                                      <  1:298522/151‑1 (MQ=255)
                                      gAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAAATTTTCTTCCTTCTAAGGTGGCCGCTTCTTACGATTATCTCCCGc                                                                                                                    >  1:116082/1‑150 (MQ=255)
                                       aaTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑g                                                                                                                 <  2:680764/151‑1 (MQ=255)
                                           cAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑t                                                                                                             <  2:386029/150‑1 (MQ=255)
                                            aGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAAAATTTCTCCTTTTTTAGTGGTTCGCCTTTCCCTTTTTTTTCCCGCC‑GTA‑tt                                                                                                            >  2:664126/1‑150 (MQ=255)
                                                  aTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAgg                                                                                                    <  1:902703/151‑1 (MQ=255)
                                                   tACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGGA‑ATAAGAATTTCTCCTTTTAGAGGGCTCCCCCTCTCCCTATTTTTTCACCCA‑GGA‑Ttt‑tttatg                                                                                                    >  1:437865/1‑146 (MQ=255)
                                                     cAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGaa                                                                                                 <  2:541005/151‑1 (MQ=255)
                                                       aTTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGACATAAAAATTCCTTCTTTCTGAGTGCGCCGCCTCTCACTATTTTTTTACCCA‑CTATTTT‑TTAAGGGa                                                                                                  >  2:282244/1‑149 (MQ=255)
                                                                        gAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATg                                                                              <  2:75486/150‑1 (MQ=255)
                                                                           taatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAAAATTTCTTCTTTTTGAGTGGGCCGCTTCTTCCTATTATTTTACCCC‑GTA‑TTT‑TTAAAGGGACGAAAAAGA‑GGAACAATGaaa                                                                           >  1:439281/1‑149 (MQ=255)
                                                                           taatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑GTAATGACATAAGAATTCTTTCTTTCATAGTGGTCCGCCTCTCCCTTTTTTCTTACCCC‑ATA‑TTT‑TTTATGGGACGACAAAGA‑GGGAACATAAAc                                                                           >  2:943307/1‑150 (MQ=255)
                                                                            aatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGACAAAAGATTTTCTCCCTTTTTAGTGCTCCCCCTCTCCCTATTATTTTACCCA‑GTA‑TTT‑TTTAAGGAAAGAAAAAGAGGGAACCagaaa                                                                            >  1:824416/1‑145 (MQ=255)
                                                                             atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑At                                                                        <  1:935870/150‑1 (MQ=255)
                                                                                       aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTAACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑Gttt                                                             <  2:519850/150‑1 (MQ=255)
                                                                                        cAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAATAAAAATTTCTTCTTTTATAGGGCTCCGCTTCTCTCTATTTTTTCTCGCA‑GTATTTT‑TTAAGGGAA‑GATAAAGA‑GGAAACATGAAC‑ATTAACCC‑Cttt                                                             >  2:127087/1‑149 (MQ=255)
                                                                                         aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCg                                                          <  1:631392/151‑1 (MQ=255)
                                                                                         aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑GTAAAGAAATAAGAATTTCTTCTTTTTGGGGGGTCCCCGTCTCACTATTATCTCACCCC‑GGA‑TTC‑TTTAGGGAAAGAGAAAGA‑GGGAACATAAAC‑CATAACCC‑Cttt                                                             >  2:785426/1‑148 (MQ=255)
                                                                                         aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑CTTAAGACAAACTATTTTCTCCTTTTCGAACGGTTCCCTTTTCTCTTTTATTTCAAGCC‑GAATTTC‑TTAAGGG‑ACGAAAAGGA‑GGGACCATGAACAATAACCCC‑Gtt                                                              >  1:452495/1‑148 (MQ=255)
                                                                                           cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAAGAAAAACCATTTTCTCCCTTTTTAGGGCTCCGCCTCTCCCTATTATTTTACCCC‑GTATTTT‑TTAAGGGAA‑GAAAAAGA‑GGGACCATAAAA‑AATAACCC‑CCTTTCGAc                                                        >  1:494763/1‑151 (MQ=255)
                                                                                                ccTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCttt                                                   <  2:60252/151‑1 (MQ=255)
                                                                                                 cTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAAGAAAAAAGATTTTCTCCCTTTCGAGTGGTCCCCCTCTCCCTATTTTTTCCCCCA‑GTATTTT‑TTAAGGGAA‑GATAAAGA‑GGGAACAATAAC‑AATAACCC‑CCTTcccgcgcctt                                                   >  1:353353/1‑142 (MQ=255)
                                                                                                 cTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAACGATTTTCTCCTTTTATAGTGCTCCCCCTCTCCCTATTTTTTCCCCCC‑GTATTTT‑TTAAGGGAA‑GAAAAAGA‑GGGAACATGAAA‑ATTAcccccctttccacgcctt                                                    >  1:946156/1‑136 (MQ=255)
                                                                                                   tAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTc                                                  <  2:258522/149‑1 (MQ=255)
                                                                                                    aaaaGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGc                                                <  1:69167/150‑1 (MQ=255)
                                                                                                     aaaGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAACGATTTTCTCCCTTTAAAGGGCTCCCCTTCTCCCTATTATTTCACGCC‑GGA‑ATCTTTAAGGGAA‑GAAAAGGA‑GGAAACATGGAC‑AATACCCCGGTTTccccgcctttccc                                                >  1:44042/1‑139 (MQ=255)
                                                                                                        gATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGAAAAACAATTTTCTCCCTTTTGAGTGCTCCCCCTCTCACTCTTATCTCCCCCC‑GCG‑TTT‑TTTAGGGAGAGCTAAAGG‑GGGACAAAGGAC‑ATTAACCC‑CCTTTCGCCGCTTTCCCCCGc                                            >  2:746789/1‑150 (MQ=255)
                                                                                                             gTCAGGCTTGCGGAGTGGCGG‑TTATGGAAATACAAATTTCTCCCTTTTGAGTGGTCCCCCTTTCCCTATTTTTTCCCCCAGGAA‑TTT‑TTAAGGGAAGACAAAGGG‑GGAACAATGAA‑‑ATTATCCC‑CTTTcccccgccttttgc                                               >  2:354509/1‑130 (MQ=255)
                                                                                                                  gCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATg                                 <  1:614913/151‑1 (MQ=255)
                                                                                                                      gCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCa                              <  1:729364/150‑1 (MQ=255)
                                                                                                                               cGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTAt                    >  1:525242/1‑151 (MQ=255)
                                                                                                                                      aGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAc               >  1:530349/1‑150 (MQ=255)
                                                                                                                                      aGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACa              <  2:335097/151‑1 (MQ=255)
                                                                                                                                         aCATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAg             <  2:480941/149‑1 (MQ=255)
                                                                                                                                           aTACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggc         >  2:550074/1‑151 (MQ=255)
                                                                                                                                             aCGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctg       >  1:705321/1‑151 (MQ=255)
                                                                                                                                             aCGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggct        <  1:827535/150‑1 (MQ=255)
                                                                                                                                               gATTTCCTCCTTTCAGAGTGCTCCGCTTCTCCCTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctggc     <  2:249469/151‑1 (MQ=255)
                                                                                                                                                aTTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGgcg    >  2:374970/1‑151 (MQ=255)
                                                                                                                                                   tcctccTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCa  <  1:370587/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                          
TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCA‑GTA‑TTC‑TTAAGGGAACGATAAGGA‑GGAACCATGAAC‑ATTACCCC‑GTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCA  >  NC_000913/1487016‑1487311

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