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→ATAglnHglutamine 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

TTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTACC  >  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)
                                              tGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGt                                                                                                     <  2:190925/150‑1 (MQ=255)
                                                     cAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGccc                                                                                            <  1:475827/151‑1 (MQ=255)
                                                          cgtcgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTGGATTTGGCACCAAAAGGGTGGCGTATTCTGCAGAACGGGGACTGAGGGGCTGGGGccc                                                                                            >  2:906125/1‑146 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGGAGGACAGAAAGTTTTGGGTATCGTGCAGAACAGGGGCTGCGCGGGTGGTGCCCCGtt                                                                                        >  2:479305/1‑147 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCat                                                                                    <  1:193146/151‑1 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCa                                                                                     <  2:41368/150‑1 (MQ=255)
                                                             cgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCa                                                                                     <  1:178552/150‑1 (MQ=255)
                                                                tACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGTAGTACAAAAAGTTTTGCGTATCGGGCCGAACGGCGGCTGCGCGGGTGGTGCCCAGTTcc                                                                                      >  2:47932/1‑146 (MQ=255)
                                                                 aCCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGGACAGAATTTTTTGCGGATCGTGCAGAAAGGCGGCTGCGCGGGTGGTGGCCAGtttc                                                                                      >  1:759580/1‑143 (MQ=255)
                                                                 aCCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAg                                                                                <  1:537378/151‑1 (MQ=255)
                                                                   cTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTGTGGAGGACAGAATATTTTGCGGATTGTGCAGGAAGGCGGCTTCGCGGGTGGGGGCCCGTTTCATATAgg                                                                               >  2:1068913/1‑150 (MQ=255)
                                                                    tttCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGTAATAAAGAAAGTTTTGGGTATCGTGCAGAAAGGGGTCTTCGCGGGTGGTGCCCAGGTTCCTATAAGCg                                                                             >  1:411088/1‑151 (MQ=255)
                                                                    tttCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTACAGAATGTTTTGCGTTTCGTGCAGAACGGCGTCTTCGCGGGTGGTGGCCCGGTTCCTATAGGCg                                                                             >  2:159196/1‑151 (MQ=255)
                                                                       cGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGtt                                                                           <  2:1116617/150‑1 (MQ=255)
                                                                         ggAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGAAAGAATGTTTGGCGTATCGTGCCGAAAGGCGTCTGCGCGGGTGGTGGCCCGCTCCCTATATGCGTTAt                                                                         >  2:705731/1‑150 (MQ=255)
                                                                          gAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTACAGAAAGTTTGGGGTATTGTGCAGAACGGGGTCTGCGCGGGTGGTGCCCCGTTTCATATAGGCGGTAt                                                                         >  1:422222/1‑149 (MQ=255)
                                                                           aaCGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGTACAGAATGTTTTGGGGATCGTGCCGAACGGCGTCTTCGCGGGTGGTGCCCAGGTCCCTATAGGGGGTAt                                                                         >  1:157020/1‑148 (MQ=255)
                                                                             cGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATg                                                                    <  1:287767/151‑1 (MQ=255)
                                                                               cAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTTGAGTACAGAAAGTTTTGCGGATCGTGGAGAAAGGCGTCTGCGCGGGTGGTGGCCCGTTCCATATAGGGGGTTTTGATGtt                                                                  >  1:13614/1‑151 (MQ=255)
                                                                                 aTACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCg                                                               <  1:189051/151‑1 (MQ=255)
                                                                                 aTACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCg                                                               <  1:250973/151‑1 (MQ=255)
                                                                                 aTACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCg                                                               <  2:300282/151‑1 (MQ=255)
                                                                                  tACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAATAAAAAAAGTTTTGCGTGTCCTGCAGAACGGCGGCTTCGCGGGTGGGGCCCCGTTCCCTATAGGCGGTTTCGATGTTTgg                                                              >  2:657715/1‑151 (MQ=255)
                                                                                     cGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGaa                                                           <  2:743555/151‑1 (MQ=255)
                                                                                     cGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGaa                                                           <  2:99577/151‑1 (MQ=255)
                                                                                       tATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGAAAGAAAGTTTTGCGTATCGTGCCGAACGGCGGCTGCGCGGTTGGTGGCCAGTTCCATATAGGCGTTATCGCTGTTCggg                                                             >  1:842918/1‑147 (MQ=255)
                                                                                       tATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGTTGAAGAACAAAATGTTTTGCGTGTCTTGCAGAACGGCGTCTTCGCGGGTTGTGCCCCGTTTCCTATAAGCGTTTTTGCTGttttgg                                                              >  2:330083/1‑143 (MQ=255)
                                                                                        aTTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTACAGAATGTTTTGGGGGTCGTGTAGAACGGGGGCTGCGCGGGTGGGGCCCCGTTCCATATAGGCGGTTTCGATGTTTggggaatg                                                        >  1:142806/1‑146 (MQ=255)
                                                                                        aTTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTg                                                        <  2:175762/151‑1 (MQ=255)
                                                                                           gctgcGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGAc                                                      <  2:971330/150‑1 (MQ=255)
                                                                                           gctgcGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACg                                                     <  1:131368/151‑1 (MQ=255)
                                                                                              gcgcTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGGAAAACAGAAAGTTTTGGGGGTAGTGCAGAAAGGCGGCTGCGCGGGTGGGGGCCCGTTCCCTATAAGGGTTATCGATGTTCGGGGAATGACGc                                                    >  1:711958/1‑149 (MQ=255)
                                                                                                  ttCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAAGACAGAATGTTTGGCGGATCCGGCAGAACGGGGGCTGCGCGGGTGGTGGCCCGTTCCCTATAGGCGTTCTCGGTGTTTGGGGACTGACGCAGATCt                                              >  2:788590/1‑151 (MQ=255)
                                                                                                      agagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAg                                          <  2:386410/151‑1 (MQ=255)
                                                                                                       gagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTTGAGTACAGAAAGGTTTGGGTATCGTGCAGAAAGGGGGCTTCGCGGGTGGGGGCCCGTTTCCTATAAGCGGTTTCGAGGGTTGGGGACTGACGGCGATCTTTTGt                                         >  2:1117234/1‑151 (MQ=255)
                                                                                                       gagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGt                                         <  1:936934/151‑1 (MQ=255)
                                                                                                         gagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAATACAGAAAGTTTTGGGTATTGTGCAGAACCGCGGCTGCGCGGGTGGTGCCCAGTTCCCTATAAGGGGTATTGATGTTTGGGGACAGACGCCGATCttt                                            >  1:311577/1‑146 (MQ=255)
                                                                                                           gTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTg                                     <  2:198490/151‑1 (MQ=255)
                                                                                                             cACCTACCGCTTTGAACTGACCGTTACCGGCGGGTTTTGTGAAGTACAGAAAGTTTTGCGGGTCGGGCAGAACGGCGGCTTCGCGGGTTGGGGCCCGTTCCCTTTAAGGGGTTTTGCGGTTCTGGGAACGACGCGGATCTTTTGtttttgt                                   >  2:1213939/1‑148 (MQ=255)
                                                                                                                  aCCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGtt                                <  2:358289/149‑1 (MQ=255)
                                                                                                                      cTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCtt                           <  1:95939/150‑1 (MQ=255)
                                                                                                                        ttGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTc                         <  1:404630/150‑1 (MQ=255)
                                                                                                                          gAACTGACCGTTACCGGCGGTTTTTATGAAGGACAAAAAGGTTTGGGTATCGTGGAGAACCGCGGCTGCGCGGGTGGGGGCCAGTTCCATATAAGCGGTATCGATGTTCTGGGAATGAAGCAGATCTTTAGTTTTGTTTTTTGCttt                          >  1:843883/1‑147 (MQ=255)
                                                                                                                          gAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTcgc                       <  1:909376/150‑1 (MQ=255)
                                                                                                                           aaCTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGt                     <  1:106/151‑1 (MQ=255)
                                                                                                                            aCTGACCGTTACCGGCGGTTTTGTTTTAGGACAGAATGTTTGGCGTAGCGGGCGGGAAGGCGGCTGCGCGGGTTGTTCCCAGTGCCATATAGGCGGTATAGATGTTTGGGGAATGTGGCCGATCTTTTGTTTTGTTGTTTTATTTcgcg                      >  1:396349/1‑149 (MQ=255)
                                                                                                                                         ggcggTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCa       <  1:884413/151‑1 (MQ=255)
                                                                                                                                              ttttGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTAcc  >  1:285537/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTACC  >  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: