Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913847,4440GT43.8% 93.4 / 32.5 64I187I (ATC→ATA) glnHglutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (6/30);  new base T (28/0);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.33e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTACC  >  NC_000913/847298‑847590
                                                                                                                                                  |                                                                                                                                                     
ttGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTttat                                                                                                                                                     >  2:11290/1‑149 (MQ=255)
 tGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTttat                                                                                                                                                     >  2:963477/1‑148 (MQ=255)
     aGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGaa                                                                                                                                                  >  1:309079/1‑147 (MQ=255)
     aGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTa                                                                                                                                               >  2:485987/1‑150 (MQ=255)
            ttCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGTAGTAAAGAATGt                                                                                                                                       >  1:1138263/1‑151 (MQ=255)
             tctcGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGtt                                                                                                                                      >  1:814259/1‑151 (MQ=255)
                 gcgcAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGGTTTTTTGAAGGAACGAAAGttt                                                                                                                                     >  2:249003/1‑148 (MQ=255)
                  cgcAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCg                                                                                                                                 <  1:643520/151‑1 (MQ=255)
                   gcAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTAAAGAATGttt                                                                                                                                     >  1:1156617/1‑146 (MQ=255)
                     aGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGGACAGAAAGTTTgggg                                                                                                                                 >  1:829908/1‑146 (MQ=255)
                       ggTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGTTGAAGGACAGAATGTTTGGCGTATc                                                                                                                             >  2:56138/1‑150 (MQ=255)
                          tttCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGGAGTAAAGAATGTTTTGGGGATCGt                                                                                                                           >  2:820999/1‑149 (MQ=255)
                             cAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGa                                                                                                                      <  2:820545/151‑1 (MQ=255)
                              aaCGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGaa                                                                                                                     <  1:1184922/151‑1 (MQ=255)
                              aaCGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGaa                                                                                                                     <  1:1141610/151‑1 (MQ=255)
                               aCGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAAc                                                                                                                    <  1:333894/151‑1 (MQ=255)
                                   gcCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTACAGGATGATTgggggat                                                                                                                              >  2:1206598/1‑132 (MQ=255)
                                          aCTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGGTTTTTTGGAGGACAAAAAGTTTTGCGGATCGTggcgaaagggg                                                                                                                >  2:825592/1‑134 (MQ=255)
                                              tGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGt                                                                                                     <  2:190925/150‑1 (MQ=255)
                                                     cAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑Gccc                                                                                            <  1:475827/151‑1 (MQ=255)
                                                          cgtcgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTGGATTTGGCACCAAAAGGGTGGCGTATTCTGCAGAACGGGGACTGAGG‑GGCTGGG‑Gccc                                                                                            >  2:906125/1‑146 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGGAGGACAGAAAGTTTTGGGTATCGTGCAGAACAGGGGCTGCGC‑GGGTGGT‑GCCCCGtt                                                                                        >  2:479305/1‑147 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCat                                                                                    <  1:193146/151‑1 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCa                                                                                     <  2:41368/150‑1 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCa                                                                                     <  1:178552/150‑1 (MQ=255)
                                                                tACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGTAGTACAAAAAGTTTTGCGTATCGGGCCGAACGGCGGCTGCGC‑GGGTGGT‑GCCCAGTTcc                                                                                      >  2:47932/1‑146 (MQ=255)
                                                                 aCCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGGACAGAATTTTTTGCGGATCGTGCAGAAAGGCGGCTGCGC‑GGGTGGT‑GGCCAGtttc                                                                                      >  1:759580/1‑143 (MQ=255)
                                                                 aCCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAg                                                                                <  1:537378/151‑1 (MQ=255)
                                                                   cTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTGTGGAGGACAGAATATTTTGCGGATTGTGCAGGAAGGCGGCTTCGC‑GGGTGGG‑GGCCCGTTTCATATAgg                                                                               >  2:1068913/1‑150 (MQ=255)
                                                                    tttCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGTAATAAAGAAAGTTTTGGGTATCGTGCAGAAAGGGGTCTTCGC‑GGGTGGT‑GCCCAGGTTCCTATAAGCg                                                                             >  1:411088/1‑151 (MQ=255)
                                                                    tttCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTACAGAATGTTTTGCGTTTCGTGCAGAACGGCGTCTTCGC‑GGGTGGTGGCCCGGTTCC‑TATAGGCg                                                                             >  2:159196/1‑151 (MQ=255)
                                                                       cGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGtt                                                                           <  2:1116617/150‑1 (MQ=255)
                                                                         ggAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGAAAGAATGTTTGGCGTATCGTGCCGAAAGGCGTCTGCGC‑GGGTGGT‑GGCCCGCTCCCTATATGCGTTAt                                                                         >  2:705731/1‑150 (MQ=255)
                                                                          gAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTACAGAAAGTTTGGGGTATTGTGCAGAACGGGGTCTGCGC‑GGGTGGT‑GCCCCGTTTCATATAGGCGGTAt                                                                         >  1:422222/1‑149 (MQ=255)
                                                                           aaCGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTACAGAATGTTTTGGGGATCGTGCCGAACGGCGTCTTCGC‑GGGTGGT‑GCCCAGGTCCCTATAGGGGGTAt                                                                         >  1:157020/1‑148 (MQ=255)
                                                                             cGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATg                                                                    <  1:287767/151‑1 (MQ=255)
                                                                               cAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTTGAGTACAGAAAGTTTTGCGGATCGTGGAGAAAGGCGTCTGCGC‑GGGTGGT‑GGCCCGTTCCATATAGGGGGTTTTGATGtt                                                                  >  1:13614/1‑151 (MQ=255)
                                                                                 aTACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑g                                                               <  1:189051/151‑1 (MQ=255)
                                                                                 aTACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑g                                                               <  1:250973/151‑1 (MQ=255)
                                                                                 aTACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑g                                                               <  2:300282/151‑1 (MQ=255)
                                                                                  tACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAATAAAAAAAGTTTTGCGTGTCCTGCAGAACGGCGGCTTCGC‑GGGTGGG‑GCCCCGTTCCCTATAGGCGGTTTCGATGTTT‑gg                                                              >  2:657715/1‑151 (MQ=255)
                                                                                     cGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGaa                                                           <  2:743555/151‑1 (MQ=255)
                                                                                     cGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGaa                                                           <  2:99577/151‑1 (MQ=255)
                                                                                       tATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGAAAGAAAGTTTTGCGTATCGTGCCGAACGGCGGCTGCGC‑GGTTGGT‑GGCCAGTTCCATATAGGCGTTATCGCTGTTC‑ggg                                                             >  1:842918/1‑147 (MQ=255)
                                                                                       tATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGTTGAAGAACAAAATGTTTTGCGTGTCTTGCAGAACGGCGTCTTCGC‑GGGTTGT‑GCCCCGTTTCCTATAAGCGTTTTTGCTGttttgg                                                              >  2:330083/1‑143 (MQ=255)
                                                                                        aTTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTACAGAATGTTTTGGGGGTCGTGTAGAACGGGGGCTGCGC‑GGGTGGG‑GCCCCGTTCCATATAGGCGGTTTCGATGTTT‑ggggaatg                                                        >  1:142806/1‑146 (MQ=255)
                                                                                        aTTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTg                                                        <  2:175762/151‑1 (MQ=255)
                                                                                           gctgcGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGAc                                                      <  2:971330/150‑1 (MQ=255)
                                                                                           gctgcGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACg                                                     <  1:131368/151‑1 (MQ=255)
                                                                                              gcgcTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGGAAAACAGAAAGTTTTGGGGGTAGTGCAGAAAGGCGGCTGCGC‑GGGTGGG‑GGCCCGTTCCCTATAAGGGTTATCGATGTTC‑GGGGAATGACGc                                                    >  1:711958/1‑149 (MQ=255)
                                                                                                  ttCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAAGACAGAATGTTTGGCGGATCCGGCAGAACGGGGGCTGCGC‑GGGTGGT‑GGCCCGTTCCCTATAGGCGTTCTCGGTGTTT‑GGGGACTGACGCAGATCt                                              >  2:788590/1‑151 (MQ=255)
                                                                                                      agagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAg                                          <  2:386410/151‑1 (MQ=255)
                                                                                                       gagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTTGAGTACAGAAAGGTTTGGGTATCGTGCAGAAAGGGGGCTTCGC‑GGGTGGG‑GGCCCGTTTCCTATAAGCGGTTTCGAGGGTT‑GGGGACTGACGGCGATCTTTTGt                                         >  2:1117234/1‑151 (MQ=255)
                                                                                                       gagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGt                                         <  1:936934/151‑1 (MQ=255)
                                                                                                         gagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAATACAGAAAGTTTTGGGTATTGTGCAGAACCGCGGCTGCGC‑GGGTGGT‑GCCCAGTTCCCTATAAGGGGTATTGATGTTT‑GGGGACAGACGCCGATCttt                                            >  1:311577/1‑146 (MQ=255)
                                                                                                           gTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTg                                     <  2:198490/151‑1 (MQ=255)
                                                                                                             cACCTACCGCTTTGAACTGACCGTTACCGGCGGGTTTTGTGAAGTACAGAAAGTTTTGCGGGTCGGGCAGAACGGCGGCTTCGCGGGTTGGG‑GCCC‑GTTCCCTTTAAGGGGTTTTGCGGTTCTGGGAAC‑GACGCGGATCTTTTGtttttgt                                   >  2:1213939/1‑148 (MQ=255)
                                                                                                                  aCCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGtt                                <  2:358289/149‑1 (MQ=255)
                                                                                                                      cTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCtt                           <  1:95939/150‑1 (MQ=255)
                                                                                                                        ttGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTc                         <  1:404630/150‑1 (MQ=255)
                                                                                                                          gAACTGACCGTTACCGGCGGTTTTTATGAAGGACAAAAAGGTTTGGGTATCGTGGAGAACCGCGGCTGCGC‑GGGTGGG‑GGCCAGTTCCATATAAGCGGTATCGATGTTC‑TGGGAATGAAGCAGATCTTTAGTTTTGTTTTTTGCttt                          >  1:843883/1‑147 (MQ=255)
                                                                                                                          gAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTcgc                       <  1:909376/150‑1 (MQ=255)
                                                                                                                           aaCTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGt                     <  1:106/151‑1 (MQ=255)
                                                                                                                            aCTGACCGTTACCGGCGGTTTTGTTTTAGGACAGAATGTTTGGCGTAGCGGGCGGGAAGGCGGCTGCGC‑GGGTTGT‑TCCCAGTGCCATATAGGCGGTATAGATGTTT‑GGGGAATGTGGCCGATCTTTTGTTTTGTTGTTTTATTTcgcg                      >  1:396349/1‑149 (MQ=255)
                                                                                                                                         ggcggTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCa       <  1:884413/151‑1 (MQ=255)
                                                                                                                                              ttttGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTAcc  >  1:285537/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGC‑GGTTGGT‑GCCCAGTTCCATATAGGCGTTATCGATGTTC‑GGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTACC  >  NC_000913/847298‑847590

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