Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,120,3680TG40.6% 109.9 / 15.8 64E73D (GAA→GACsolAN‑methyltryptophan oxidase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (8/30);  new base G (26/0);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 3.83e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCCCTTGCGCATCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGTTTAAACCGGC  >  NC_000913/1120221‑1120517
                                                                                                                                                   |                                                                                                                                                     
aaTCCCTTGCGCATCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTc                                                                                                                                                     <  2:652093/150‑1 (MQ=255)
   cccTTGCGCATCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGGTCATcc                                                                                                                                                 >  2:581412/1‑151 (MQ=255)
        gcgcATCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCcaca                                                                                                                                             <  1:297827/150‑1 (MQ=255)
        gcgcATCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAg                                                                                                                                            <  2:33580/151‑1 (MQ=255)
           cATCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTCATCCCCCAGc                                                                                                                                           >  1:932169/1‑149 (MQ=255)
              cGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGGTCATCCCCCAACGttt                                                                                                                                       >  1:451387/1‑150 (MQ=255)
              cGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTg                                                                                                                                      <  1:837869/151‑1 (MQ=255)
               gAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTCATTCCCCAACGgtt                                                                                                                                       >  2:496987/1‑149 (MQ=255)
                aGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGGTCTTCCcccaac                                                                                                                                           >  2:620849/1‑142 (MQ=255)
                   ttttCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTCATCCCCCAGCGTTTgcg                                                                                                                                    >  1:310651/1‑148 (MQ=255)
                       cAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCgga                                                                                                                              <  1:404648/150‑1 (MQ=255)
                        aaCGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCggagg                                                                                                                            <  2:359133/151‑1 (MQ=255)
                          cGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTCATTCCACAGCGGTTTTGCGCggg                                                                                                                              >  2:567565/1‑146 (MQ=255)
                             tGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCa                                                                                                                        <  2:147385/150‑1 (MQ=255)
                             tGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCa                                                                                                                        <  1:350491/150‑1 (MQ=255)
                              gAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTTATCCCACAACGTTTGCGCGCggggg                                                                                                                            >  1:1022254/1‑145 (MQ=255)
                              gAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTTAACCCCCCCCGCTTGTGCGCggggg                                                                                                                            >  2:446005/1‑145 (MQ=255)
                               aGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCg                                                                                                                     <  1:1128939/151‑1 (MQ=255)
                                   gCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGTGCATTCCACCGCGCTTGCGGGCGGGGGGGCCGCggg                                                                                                                   >  2:757999/1‑149 (MQ=255)
                                       ttGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGAGGTCCTCCCCCAACGTTTTCGCGCGGGGGGGCAGCggg                                                                                                                   >  2:128882/1‑145 (MQ=255)
                                           tCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGGTTATCCCCCACAGTTTGTGCGCGGGGGGACCGCGGGGCATAAtttt                                                                                                         >  2:1076835/1‑151 (MQ=255)
                                            cGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTc                                                                                                        <  1:328303/151‑1 (MQ=255)
                                               cgcTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTTATCCCACAGCGGTTGTGCGCGGGGGGGCCGCGGGGGATAAttttttcc                                                                                                     >  1:311701/1‑147 (MQ=255)
                                                   gtgGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCtt                                                                                                   <  1:918759/149‑1 (MQ=255)
                                                   gtgGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTc                                                                                                  <  2:937351/150‑1 (MQ=255)
                                                          aCGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGGTTATCCCACAGCGGTTTGGCGCGGGGGGGCAGAGGGACATACTTTTTGCCtttccc                                                                                               >  2:793535/1‑142 (MQ=255)
                                                               ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATgg                                                                                     <  1:458146/151‑1 (MQ=255)
                                                                gCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGAGTTTATCCCACAGCGCTTTCGCGCGGGGGACCAGCGGGGCATACATTTCTCCTTCCCCCTaacc                                                                                         >  2:644120/1‑144 (MQ=255)
                                                                gCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGc                                                                                    <  2:817813/151‑1 (MQ=255)
                                                                  gagaAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGa                                                                                  <  2:29417/151‑1 (MQ=255)
                                                                     aaaTGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAAtt                                                                               <  2:1077017/151‑1 (MQ=255)
                                                                       aTGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTaa                                                                             <  2:500887/151‑1 (MQ=255)
                                                                         gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAAt                                                                            <  2:650584/150‑1 (MQ=255)
                                                                         gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATc                                                                           <  2:77802/151‑1 (MQ=255)
                                                                         gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATc                                                                           <  1:2210/151‑1 (MQ=255)
                                                                          tgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGGTCATCCCCCAGCGCTTTTGCGCGGGGGGGCAACGGGGCATACTTTTTTCCTTCCCCACAAGCATTGCGAAt                                                                                >  1:575047/1‑145 (MQ=255)
                                                                           ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATcgc                                                                         <  2:139974/151‑1 (MQ=255)
                                                                            gAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGGGGTTATCCCCCAACGTTTGCGCGCGGGGGAGCAGCGGGGCAAACTTTTTTCCTTCACCCTAAACATGGCGAATTAATcgc                                                                         >  1:382743/1‑150 (MQ=255)
                                                                              aTCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGAGGTTCTTCCACAGAGGTTTTGCGCGGGGGGACCGCGGGGGAGAATTTTTTCCTTTTCCCTAAACATGGGGAAt                                                                                >  2:980232/1‑141 (MQ=255)
                                                                               tCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGAGGTCATCCCCCAGCGGTTTTGCGCGGGGGGCGGGCGGGACAAAATTTTTGCCCTCCCCCTCAACACGGCGGGTTAATccc                                                                         >  2:478131/1‑145 (MQ=255)
                                                                                cAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATc                                                                    <  1:177253/151‑1 (MQ=255)
                                                                                     gggCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTTATCCCACCGCCTTTGTGCGCGGGGGGGCCGCGGGGAAAAATTTTTGCCTTCACCCTAAACATGGCGGATTAATCCCGTATCGcc                                                                 >  1:1152094/1‑149 (MQ=255)
                                                                                        ccAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTCATCCCACAGCGGTTTCGCGCGGGGGGCCAGCGGGGCATACTTTTTTCCCTCACCATAAGCAAGGCGAATTAAACGCGTATCGCCgggg                                                             >  1:335873/1‑150 (MQ=255)
                                                                                          aaGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGGGAGTTTATCCCCCACCGGTTTTGCGCGGGGGGCCAGCGGGGCATACTTTTCTCCTTCACCCTAAACAAGGGGGATTAATCGCGGATCGCCCtggggg                                                          >  1:865039/1‑151 (MQ=255)
                                                                                              ttAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTTATCCCCCCCCGCTTTCGCGCGGGGGGCCAGCGGGGGCAACTTTTTGCCCTCACCCTAAACATGGCGGATTAAACTCGGATCGCCGTGGGGGGTg                                                       >  1:63089/1‑150 (MQ=255)
                                                                                                        cAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATg                                            <  1:320078/151‑1 (MQ=255)
                                                                                                        cAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTTATTCCACAACGGTTTCGCGCGGGGGGGCAGCGGGGCAAAATTTTTTCCCTTACCATAAACAAGGGGGATTAAATGCGGATCCCCCGGGGGGGCGCCgtggtggt                                             >  2:666381/1‑148 (MQ=255)
                                                                                                        cAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTCATCCCACAACGGTTTCGGGCGGGGGGCCAGCGGGGGATACTTTTTCCCCTCCCCCTAAACAAGGCGGAATTATTGCCTATCGTCGGGGTGGCGGcccggtttgt                                             >  2:751875/1‑142 (MQ=255)
                                                                                                         agagCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATTCCACAGCGGTTTCGCGCGGGGGAGCCGCGGGACATACTTTTTGCCTTCACCCTAAACATGGCGGAATAAAAGCGGATCGCCGTGGTGGCTGCCGGGTTGGTg                                            >  1:890115/1‑150 (MQ=255)
                                                                                                           agCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGGTCATCCCACAGCGTTTTTGCGCGGGGGGCCAGCGGGGCATAATTTTTTCCTTCTCCCTAAAAAAAGCGAATTAATTGCGTATCCCCCGGGGGGGTGCCGCGTTGGTGCgg                                         >  1:997881/1‑151 (MQ=255)
                                                                                                             cGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTg                                       <  1:240569/151‑1 (MQ=255)
                                                                                                                 cAAAAATGGGATCTTCTTCGTTGTGGCGGGGGAGGTCATTCCCCCGCGGTTTCGCGCGGGGGGGCAGCGGGGCATAATTTTTGCCTTCACCATAAACAAGGCGGATTAATCGCGTATCGCCGCGGGGGCTGCCGTGTTGAttgggggggat                                   >  1:76578/1‑141 (MQ=255)
                                                                                                                    aaaTGGGATCTTCTTCGTTGTGGCGGGGGAGTTTATCCCACACCGGTTGTGCGTCGGGGGCCAGCGGGGGATACCTTTTGCCTTCACCACAACAACGGCGAATTAAACTGGTACCACCGTGGGGGGTGGCGCGTTGGTGccggggg                                     >  2:316133/1‑140 (MQ=255)
                                                                                                                    aaaTGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATg                                <  2:248552/151‑1 (MQ=255)
                                                                                                                        gggATCTTCTTCGTTGTGGCGGGAGAGGTTATCCCACAACGGTTGCGGGCGGGGGGGCAGCGGGGCATACTTTTTGCCTTCCCCCCAAACAAGGCGAAATAATTGCGGAACGCCGCGGTGGCGGCCGGGGTGGTGGGGGGGGATATtgggg                            >  1:593919/1‑149 (MQ=255)
                                                                                                                          gATCTTCTTCGTTGTGGCGGGAGAGGTCCTCCCCCCGCCGGTGCGCGCGGGGGGCCCGGGGGGCGTAATTTTTTCCCTCCCCACAAACAAGGGGGAATTATATCGTATCGCCCCGGGGGCTGCCCGGTTTGTGGGGGGGGATATGGGCg                            >  2:370661/1‑149 (MQ=255)
                                                                                                                                   cGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTa                  <  1:947852/150‑1 (MQ=255)
                                                                                                                                       gtgGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACg             <  1:4379/151‑1 (MQ=255)
                                                                                                                                         ggCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGtt           >  1:731373/1‑151 (MQ=255)
                                                                                                                                            gggAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGTTTaa        <  2:1047865/151‑1 (MQ=255)
                                                                                                                                            gggAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGTTTa         >  2:1106933/1‑150 (MQ=255)
                                                                                                                                               agagTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGTTTAAAcc     >  1:153338/1‑151 (MQ=255)
                                                                                                                                                  gTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGTTTAAACCGgc  <  2:126260/151‑1 (MQ=255)
                                                                                                                                                  gTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGTTTAAACCGgc  <  1:920982/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATCCCTTGCGCATCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGGGAGAGTTCATCCCACAGCGTTTGCGCGCGGAGGACCAGCGGGACATACTTTTCGCCTTCACCATAAGCATGGCGAATTAATCGCGTATCGCCGTGGTGGCTGCCGTGTTGATGCGGTGGCATATGGGCGTCGGTCATTAGCACGTTTAAACCGGC  >  NC_000913/1120221‑1120517

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