Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4940TG42.0% 127.7 / 11.5 70intergenic (‑512/+91)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (7/33);  new base G (29/0);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 3.54e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTT  >  NC_000913/987352‑987640
                                                                                                                                               |                                                                                                                                                    
aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATa                                                                                                                                             <  2:583720/151‑1 (MQ=255)
 gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAg                                                                                                                                            <  1:527538/151‑1 (MQ=255)
   cATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGAt                                                                                                                                          <  2:393894/151‑1 (MQ=255)
     tCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑aaaaaaa                                                                                                                                           >  2:8887/1‑144 (MQ=255)
       cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAATaaaattct                                                                                                                                      >  2:424690/1‑144 (MQ=255)
       cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAGATTgtt                                                                                                                                      >  1:806248/1‑149 (MQ=255)
        aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGT‑AAATAAATTGCtt                                                                                                                                     >  1:842435/1‑151 (MQ=255)
         gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTg                                                                                                                                    <  2:320931/151‑1 (MQ=255)
         gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAGATTGttt                                                                                                                                     >  2:146753/1‑150 (MQ=255)
          aaTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGG‑AAATaaagttatt                                                                                                                                     >  2:768798/1‑141 (MQ=255)
              ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTa                                                                                                                               <  1:201500/151‑1 (MQ=255)
              ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTa                                                                                                                               <  1:1072964/151‑1 (MQ=255)
              ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑AAATAAATTTCttttttt                                                                                                                                >  1:821164/1‑145 (MQ=255)
                   aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑AAAAAAAATTCTTTACTTaacaa                                                                                                                           >  1:763249/1‑150 (MQ=255)
                     aTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAATAAATATCTTTACTTtacaaa                                                                                                                          >  1:386981/1‑149 (MQ=255)
                      tAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAATAAATTAGTTTTCTTtacaaaca                                                                                                                        >  2:856280/1‑150 (MQ=255)
                          aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGT‑AAAAAGAATTCTTGTCTTaaaaaaaaaaaa                                                                                                                    >  2:534600/1‑144 (MQ=255)
                          aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAAATTATTTTACTTaaaaaaaaaaaa                                                                                                                    >  1:617418/1‑144 (MQ=255)
                                aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGAc                                                                                                            <  2:37910/151‑1 (MQ=255)
                                   gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGa                                                                                                         <  2:795470/151‑1 (MQ=255)
                                    tgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTGGT‑AAAAAAAATTTTTGAATTaacaaaaaaaaaaatttaaa                                                                                                            >  1:85012/1‑134 (MQ=255)
                                    tgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAATAAATTTGTTTACTTTACAAAACA‑ATAAATTTACAg                                                                                                          >  2:773835/1‑149 (MQ=255)
                                        aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAg                                                                                                     <  1:757700/150‑1 (MQ=255)
                                              tAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGA‑AAAAAAATTGTTTTA‑TTTACAAAAAA‑AAAAATTTACAGAaaaaaaaaaa                                                                                               >  2:696355/1‑150 (MQ=255)
                                               aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATaaaaaa                                                                                             <  2:629598/151‑1 (MQ=255)
                                               aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATaaaaa                                                                                              <  2:1070083/150‑1 (MQ=255)
                                                 gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAag                                                                                           <  1:97215/151‑1 (MQ=255)
                                                 gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑GAAAAAAATGTTTGAGTTAACAAAAAA‑AAAAAATTACAGAAAAGATaaaaaaa                                                                                            >  2:358903/1‑150 (MQ=255)
                                                  aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAaga                                                                                          <  1:609993/151‑1 (MQ=255)
                                                  aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAaga                                                                                          <  2:1063494/151‑1 (MQ=255)
                                                   aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATATAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAagag                                                                                         <  2:103351/151‑1 (MQ=255)
                                                    aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAAAATGCTTGACTTTACAAAAAG‑ATAAATTGACAGaaaaaaaaaaaaaaaa                                                                                          >  2:234948/1‑147 (MQ=255)
                                                     aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                       <  2:504024/151‑1 (MQ=255)
                                                     aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                       <  2:779173/151‑1 (MQ=255)
                                                      aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                      <  2:261248/151‑1 (MQ=255)
                                                      aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                      <  2:690902/151‑1 (MQ=255)
                                                       aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑CAAAAAATTTGTTTGCTTTAAAAAAAG‑AAAAATTTACAGAAAAAAAAAAAAAAAAGACa                                                                                      >  2:1093363/1‑150 (MQ=255)
                                                        cTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGT‑AAAAAGAAAGGTTTGGTTAAAAAAAAA‑AAAAATTGACAGAAAAGATaaaaaaaaa                                                                                          >  2:956798/1‑143 (MQ=255)
                                                        cTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAATAAATATC‑TGACTTTACAAAAAA‑AAAAATTTGCAGGAAAGAaaaaaaaaaa                                                                                          >  2:317885/1‑143 (MQ=255)
                                                         tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGATGA‑GAAAAAAATGCTTGTCTTAAAAAACAG‑AAAACAAGTTAGAAAAGATaaaaaaaaaaa                                                                                        >  1:202962/1‑142 (MQ=255)
                                                          aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc                                                                                    <  1:93169/149‑1 (MQ=255)
                                                          aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑AAAAAGAATGCTTTACTTAAAAAAAAA‑AAAACTTTACAGAAAAAAAAAAAAAAAAAACAGCgg                                                                                  >  2:859143/1‑151 (MQ=255)
                                                             ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                <  1:789817/150‑1 (MQ=255)
                                                                 aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑AAATAAATTGCTTTACTTTACAAACAA‑AAAAATTTACAGAAAAAAAAAAAAAAAAACCAGCGGGGGCTgg                                                                           >  2:1016652/1‑151 (MQ=255)
                                                                     aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑AAAAAAAATGCTTTACTTAACAAAAAA‑AAAAAATTACAGAAAAAAAAAAAAAAAAAACAGCGGGGGGTgg                                                                           >  2:725538/1‑147 (MQ=255)
                                                                     aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGGGC‑TAAAAGAATGGGTTACTTTTAAAAAAA‑AAAACTTTACACCAaaaaaaaaaaaaaaaaccgcg                                                                                   >  2:102201/1‑130 (MQ=255)
                                                                          gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                   <  2:92418/150‑1 (MQ=255)
                                                                            ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAGATTGCTTTATTTTACAAAAAG‑AAAACTTTACAGAAAAAATAAAAAAAAAACCAGCGGGGGCTggggctt                                                                      >  2:683668/1‑140 (MQ=255)
                                                                              aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑AAAAAAATTGCTTTACTTTAAAAAAAA‑ATAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCGGGCTCtt                                                                      >  2:324622/1‑143 (MQ=255)
                                                                                 aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                            <  1:638479/150‑1 (MQ=255)
                                                                                 aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                            <  2:538057/150‑1 (MQ=255)
                                                                                   tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAAAATGTTTTATTTTATAAACAA‑ATAAATTGACAGAAAAAATAAAAAAAAAAACAGCGGGCGGCGGGGCTTGTAGAAGGTaaa                                                          >  1:840434/1‑150 (MQ=255)
                                                                                   tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑CAATAGATTGCTTTGCTTTACAAAAAA‑AAAAATTTACAGAAAAGAAAAAAAAAAAACCAACGGGGGCTGGCTCTTTTAGAAGCCaaa                                                          >  2:636504/1‑150 (MQ=255)
                                                                                    ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGC‑AAAAAGATTTATTTGTTTTACAAAAAGAAAAAATTTAAAAAAAAAAAAAAAAAAAAG‑CAACGGGCGCGGGCTCtttgggaagcgaaa                                                          >  2:1088455/1‑137 (MQ=255)
                                                                                      ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                       <  1:367538/150‑1 (MQ=255)
                                                                                              aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                              <  2:55304/151‑1 (MQ=255)
                                                                                                aTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGT‑AAATAGAATTCTTTACTTTAAAAACAA‑AAAAATTTACAGAAAAAAAAAAAAAAAAAACAACGGGCGGTGGGTCTTGAAAAAGCTAAATTAAAAGATgg                                               >  2:888/1‑148 (MQ=255)
                                                                                                  gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                           <  2:613019/150‑1 (MQ=255)
                                                                                                   ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAAAATGCTTTACTTTAAAAAAAA‑AAAAATTTACAGAAAAAAAAAAAAAAAAAACAGAGGGGGCGGGGTCTTGAAAACGCCAAATTAAAAACTGGCggt                                           >  2:678512/1‑147 (MQ=255)
                                                                                                          gTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGG‑CAATAAATTTCTTGACTTAACAAAAAA‑AAAACATTACAAAAAAGAAAAAAAAAAAAACAACGGGCGCCGGCTCTTGAAGACGCTAAATTAAAAGCTGGGGtttaccgg                                     >  2:569119/1‑142 (MQ=255)
                                                                                                            aaTATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTAc                                <  2:615622/151‑1 (MQ=255)
                                                                                                               atTCCGTAACAGGATGATCGTTATGGC‑GGGGTGTG‑AAAAAGATTTCTTTACTTAAAAAACAA‑ATAACTTTACAGAAAAAAAAAAAAAAAAGAAAACGGGCGCCGGCTTTTTTAGACGCTAAAATAAAAACAGGCGGTACGCGCAgg                                  >  1:763238/1‑145 (MQ=255)
                                                                                                                 tCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                            <  1:903187/150‑1 (MQ=255)
                                                                                                                      aaCAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACg                       <  1:295520/150‑1 (MQ=255)
                                                                                                                       aCAGGATGATCGTTATGGC‑GGGGGTGG‑AAAAAAAATTGTTGACTTTACAAAAAA‑AAAACATGATCGAAAAAAAAAAAAAAAAAACAACGGGGGCCGGCGCTTTGAAAAGCTAAATTAAAAACTGGGGTTATGCGGaagaccgggggacgg                      >  2:1038494/1‑136 (MQ=255)
                                                                                                                                     aTGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac       <  1:983710/151‑1 (MQ=255)
                                                                                                                                         c‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt   >  1:67014/1‑151 (MQ=255)
                                                                                                                                         c‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt   >  1:804000/1‑151 (MQ=255)
                                                                                                                                         c‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt   <  1:575101/151‑1 (MQ=255)
                                                                                                                                           ggggTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt  <  1:428410/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGC‑AAATAGATTGCTTGACTTAACAAACAG‑ATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTT  >  NC_000913/987352‑987640

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