Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913345,3040GA43.3% 74.0 / 35.4 60intergenic (+313/‑100)yahL/yahMuncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (22/12);  new base A (0/26);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 4.80e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTAT  >  NC_000913/345163‑345452
                                                                                                                                             |                                                                                                                                                      
cAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGc                                                                                                                                              >  2:605139/1‑151 (MQ=255)
       gCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAg                                                                                                                                       <  2:209357/150‑1 (MQ=255)
            gtATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCat                                                                                                                                 <  2:506597/151‑1 (MQ=255)
            gtATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCat                                                                                                                                 >  2:163795/1‑151 (MQ=255)
              ataAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCatat                                                                                                                               >  1:288864/1‑151 (MQ=255)
              ataAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCata                                                                                                                                <  2:135755/150‑1 (MQ=255)
                aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCa                                                                                                                             <  1:42321/151‑1 (MQ=255)
                aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCa                                                                                                                             <  2:161771/151‑1 (MQ=255)
                 aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCaa                                                                                                                            <  1:82614/151‑1 (MQ=255)
                  aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCaaa                                                                                                                           <  1:76865/151‑1 (MQ=255)
                     aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAAc                                                                                                                        <  2:502312/151‑1 (MQ=255)
                       gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCg                                                                                                                      >  1:437033/1‑151 (MQ=255)
                       gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGCGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAAcc                                                                                                                       <  1:217754/150‑1 (MQ=255)
                            cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGc                                                                                                                     <  1:691411/147‑1 (MQ=255)
                            cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCAtt                                                                                                                 >  1:205501/1‑151 (MQ=255)
                               gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGc                                                                                                              >  2:172144/1‑151 (MQ=255)
                               gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGc                                                                                                              >  1:22536/1‑151 (MQ=255)
                                         tgtgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa                                                                                                     >  1:440432/1‑150 (MQ=255)
                                          gtgAATGTGTTGGTTGACATTCATAAAAAAAAAATCAAAATCCCATCATGTTTGTGGGGGTTTTTTTCTTATTCTTACCCACCGCACCCAAGCCTAACAAAAAA‑TATGC‑TGTAGGTCTCATATAAAAAACCGCCATTAGCTCATCAGGaa                                                                                                    <  2:159521/150‑1 (MQ=255)
                                             aaTGTGTTGGTTGCCTTTCATAAAAAAAAAATCATAATCCCATCATTTTTGTGTGGGGTTTTTTCTGTATCTTACGCTCCGCACTCAACCCTGAAAAAAAA‑TATGC‑TGTAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAag                                                                                                   <  1:248725/148‑1 (MQ=255)
                                                   ttggttgACATTCATATAAAAAAAACCATAATCCCATCTTGTTTGTGGGGGTTTTTTTCTTTATCTTACGCACCGCACTCAGGCCTGACAAAAAA‑TATGC‑TGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacg                                                                                           <  1:296885/150‑1 (MQ=255)
                                                     gtttgACATTCATAAAAAAAAAATCATAATCCCATCATGTTTGTGTGGGGTTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAA‑TATGC‑TGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                            <  1:11421/145‑1 (MQ=255)
                                                       ttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                      >  1:139643/1‑151 (MQ=255)
                                                       ttgACATTAATATAAAAAAAACCTAAATCCCATCATTTTTGTGGGGGTTTTTTTCTTTATCTTACCCACCCCACACAAGCCTAAAAAAAAA‑TATGC‑TGAAAGCCTCAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACaa                                                                                       <  2:34230/150‑1 (MQ=255)
                                                            aTTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                  >  1:458046/1‑150 (MQ=255)
                                                                atatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                             >  2:506429/1‑151 (MQ=255)
                                                                  aaaaaaaaaaTCAAAATCCCATCTTGTTTGGGGGGGGTTTTTTCTGTATTTTCCGCATCGCCCTCAAGCCGGAAAAAAAA‑TATGC‑TGAAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttg                                                                              <  1:271546/145‑1 (MQ=255)
                                                                   tGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                           >  1:244135/1‑150 (MQ=255)
                                                                     aaaaaaaaCAAAATCCCATATTGTTTGGGGGGGTTTTTTTCTTTCTCTACCCAACCGCCCCCAACCCAAAAAAAAAA‑TATGC‑GGTAGGGCTCATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                           <  2:359487/148‑1 (MQ=255)
                                                                     aaaaaaaTCATAATCCCATCATTTTTGTGTGGGTTTTTTTCTGTATCTTCCGCACCGCACTAAAGCCGGAAAAAAAA‑TATGC‑TGTAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                        <  1:458052/151‑1 (MQ=255)
                                                                     aaaaaaaTCAAAATTCCATCTTTTTTGTGTGTGGTTTTTTCTTTATCTTACGCACCCCACACAACCCTAACAAAAAA‑TATGC‑TGTAAGACTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                          <  2:491525/149‑1 (MQ=255)
                                                                         aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                     >  1:284804/1‑150 (MQ=255)
                                                                             cATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACggg                                                                >  1:164069/1‑151 (MQ=255)
                                                                              aaaaTCCCATAATTTTTGGGGGGGTTTTTTTTTTTTTCTCCCCCACCGCACTCAAGCCTAAAAAAAAA‑TATGG‑TGAAAGCCTAATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg                                                                  <  2:186138/146‑1 (MQ=255)
                                                                                atgcctatcatGTTTGTGGGGGTTTTTTTCTTTTTCTTCCCCCTCGCACTAAACCCTGA‑AAAAAA‑TATTCTTGAAAGACCGACTTCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                               <  1:528433/143‑1 (MQ=255)
                                                                                  ttCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                           >  2:289285/1‑151 (MQ=255)
                                                                                    ccaccttttttGGGTGGGGTCTTTTCGGTATATTACACACCCCCCCCAAGCTGGACAGAAAA‑AATGC‑TGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGa                                                          <  1:335336/142‑1 (MQ=255)
                                                                                     catcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                        >  2:230326/1‑151 (MQ=255)
                                                                                          tttttGTGGGGGTTTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAA‑TATGC‑TGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                     <  1:34308/147‑1 (MQ=255)
                                                                                          tGTTGGTGGGGGGTCTTTTCTGTTTCTTCCGCACCGCACACAACCCGAAAAAAAAA‑GATGC‑TGAAGGGCAAATATCAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGcc                                                      <  2:356646/148‑1 (MQ=255)
                                                                                            tttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGa                                                  >  1:588725/1‑150 (MQ=255)
                                                                                            tttGTGGGGGTTTTTTTCTGTATCTTACCCATCGCACTAAACCCTGACAAAAAA‑TATGT‑TGTAAGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                    <  1:432498/148‑1 (MQ=255)
                                                                                              tgtgtgGGGTTTTTTCTGTTTCTTACACACCGCACACAACCCTGACAAAAAA‑TATGC‑TGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                <  2:237222/150‑1 (MQ=255)
                                                                                                   gggtttttttCTGTTTCTAACCCCCCGCCCTAAACCGTGAAAAAAAA‑TATGC‑TGAAGCGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                          <  2:54305/145‑1 (MQ=255)
                                                                                                   ggggTTTTTTTTTTATTTTACCCATCCCACTAAAGCCTGAAAAAAAA‑TATGT‑TGTAGGGCAAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                <  1:541779/145‑1 (MQ=255)
                                                                                                    gttttttttttttttCTTCCGCATCTCACCAAAGCCTAAAAAAAAA‑AATGT‑TATAGTTCAAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcc                                         <  2:681185/137‑1 (MQ=255)
                                                                                                       tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                      >  2:279676/1‑151 (MQ=255)
                                                                                                         ttttCTATATCTCACGCACCCCAC‑CAAACCAGAAAAAAAATTATGA‑GGTAGGCCTCATTTAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                               <  2:407204/140‑1 (MQ=255)
                                                                                                           ttCTTTATCTTCCGCATCGCACTCAAGCCTAACAAAAAA‑TATGC‑TGTAAGGCTAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                   <  1:515870/150‑1 (MQ=255)
                                                                                                            tCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTc                                 >  1:263479/1‑151 (MQ=255)
                                                                                                                  tCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                           >  2:64820/1‑151 (MQ=255)
                                                                                                                            cGCACTCAAGCCAGAAAAAAAA‑TATGC‑TGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGaa                   <  1:556776/149‑1 (MQ=255)
                                                                                                                            cGCACTAAAGCCAGACAAAAAA‑TATGC‑TGTAAGTCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt                  <  1:603422/150‑1 (MQ=255)
                                                                                                                               aCTCAACCTTGAAAAAAAA‑TATGC‑TGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                <  1:224850/149‑1 (MQ=255)
                                                                                                                                cccAAGCCTGACAAAAAA‑TAGGC‑TGTAAGCCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa              <  1:718359/148‑1 (MQ=255)
                                                                                                                                 tCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                 >  1:360867/1‑146 (MQ=255)
                                                                                                                                   aaGCCTGAAAAAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt             <  1:29423/148‑1 (MQ=255)
                                                                                                                                    accgtaaCAAAAAA‑TATGT‑TTTAGGGTTAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt             <  2:690350/141‑1 (MQ=255)
                                                                                                                                           cAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat  >  2:29877/1‑151 (MQ=255)
                                                                                                                                           cAGAAAA‑TATGC‑TGTAAGCCTCATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatt    <  2:602859/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
CAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTAT  >  NC_000913/345163‑345452

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