Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1520AG26.4% 152.2 / 21.9 87intergenic (+135/‑63)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (45/19);  new base G (0/23);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.95e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAT  >  NC_000913/3781053‑3781246
                                                                                                    |                                                                                               
gggAGAGGAGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                                 >  2:238409/1‑100 (MQ=255)
  gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                               >  1:81515/1‑100 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                             <  1:2482769/100‑1 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                             <  2:2673560/100‑1 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                             <  1:1138662/100‑1 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                          >  1:257849/1‑100 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                          <  2:2045527/100‑1 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                          >  1:200350/1‑100 (MQ=255)
       gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                         <  1:2088184/100‑1 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                        >  1:2133575/1‑100 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                       <  1:1741815/100‑1 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                       <  1:1741938/100‑1 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                                     <  1:192232/100‑1 (MQ=255)
            cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                       <  1:2198222/97‑1 (MQ=255)
            cGTTCGGCCCTGTATGTACTCCCTCGCTCTGAAACGACACCGCACGCTTTTTTTCTCCCGCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                                                    <  2:607765/100‑1 (MQ=38)
              ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGc                                                                                  >  2:423217/1‑100 (MQ=255)
                ccgcgcTGTATGTACTCCCTCACTCTGAAACCACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                                >  1:1813878/3‑100 (MQ=255)
                 ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACt                                                                               >  1:306824/1‑100 (MQ=255)
                 ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACt                                                                               >  2:777541/1‑100 (MQ=255)
                   cgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGc                                                                             >  2:87275/1‑100 (MQ=255)
                    ncTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCTTCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                            >  2:1973384/2‑100 (MQ=255)
                       gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTTTTTTCTCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                         >  2:1875493/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                     >  2:2687268/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                     >  2:931777/1‑100 (MQ=255)
                             ccgcccTCACTCGAAAACAACCCCCCTCTTTTTTTTTCCCCCCCCCCCGGGGGGGGGGGGGGCGGGGGTGGGGGAAAAAGGCCGCACTGCTCAcc                                                                        <  1:1768350/92‑1 (MQ=255)
                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                                   >  2:1592138/1‑100 (MQ=255)
                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                                   >  1:289573/1‑100 (MQ=255)
                                ccctccctCGAAAGCGACCCCTCTCTTTTTTTTTCCCCCGCGCCCCTGCGGGGGGGGGCCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                 <  1:669988/99‑1 (MQ=255)
                                ccccccccctGAAGCGACCCCGCCTTTTTTTTTTCCCCCCCGCCCCCCGGGGGGGGGGGCGGGGTGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                <  1:1114753/92‑1 (MQ=255)
                                cccTCCCCCTGAAACGACACCGCTCTTTTTTTTTCCCCCTCGCCCCCCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                                   <  2:1554204/97‑1 (MQ=37)
                                cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                >  1:853067/1‑100 (MQ=255)
                                 cccccctctGAACCGCCCCGGCCCTCTTTTTTTCCCCCTCCCCCCTCCGGGGAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                               <  2:1016746/95‑1 (MQ=255)
                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                               >  2:337479/1‑100 (MQ=255)
                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                               >  1:2104139/1‑100 (MQ=255)
                                 ccTCACTCTGAAACGACACCACACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCGGGGTGAGGATAAAAGGCCGCACTGCTCACCAATTTTc                                                                 >  2:1194264/1‑98 (MQ=255)
                                 ccTCACTCGAAAACGACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                               <  1:1932069/100‑1 (MQ=255)
                                   tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                             >  2:1711085/1‑100 (MQ=255)
                                     cctcgGAACCGACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                              <  2:2018216/96‑1 (MQ=255)
                                       tctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                         >  1:518909/1‑100 (MQ=255)
                                       tataaaaCTACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCCGGGGGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                         <  1:1105424/96‑1 (MQ=255)
                                        ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                        >  2:1968622/1‑100 (MQ=255)
                                          gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                      >  2:395785/1‑100 (MQ=255)
                                           aaaCCACCCCGCACTTTTTTTTTCCCCCCCGCCCCCCCGGGGAGAGGGCCGTGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                       <  2:1495741/98‑1 (MQ=255)
                                                acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                        >  2:1704475/1‑92 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                             >  1:2487016/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                             >  1:980923/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                             >  1:508464/1‑100 (MQ=255)
                                                      cACTCTTTTTTCCCCCCTCGCCCCCCGGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                              <  2:1406017/96‑1 (MQ=255)
                                                        ctctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                        >  1:1175043/1‑100 (MQ=255)
                                                        ctctTTTTTTCTCCCCCCCCCCCCCGGGGGGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                         <  2:1433862/99‑1 (MQ=255)
                                                           tttttttCCCCCCCCCCCCCCCGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                        <  1:2157465/97‑1 (MQ=255)
                                                           tttttttCCCCCCCCCCCCCCCGGGGGGGGGCCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                              <  2:2064707/91‑1 (MQ=255)
                                                           ttttttcCGCCCCCCCCCCCCCGGGGGGAGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                      <  2:1871177/99‑1 (MQ=255)
                                                            ttttttCTCCTCCGCCCCGCGGGGGGGGGGGCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                        <  2:124733/96‑1 (MQ=255)
                                                             tttttCTCCCCCGCCCCCCCGGGGGGAGGGCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                   <  2:1380291/100‑1 (MQ=255)
                                                               tttCCCCCCCGCCCCTCCGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAgga                                 <  2:306821/100‑3 (MQ=255)
                                                               tcccccccTCCCCCCGCGGGGGAGGGGGCGGGGTGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                                 <  1:2130933/95‑1 (MQ=255)
                                                                tcccccccccccccGCCGGGGGGAGGGCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                <  1:605424/90‑1 (MQ=255)
                                                                  ctcCCTCGCCCCCCCGGGGGGAGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                              <  2:2069749/100‑1 (MQ=255)
                                                                   tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                             >  1:2458212/1‑100 (MQ=255)
                                                                   tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                             >  1:2413242/1‑100 (MQ=255)
                                                                    cccTCTCTCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACGAATTTTCTGCTATCCTGCCCCCGCACTAAGGTGGCAGTAttc                              >  2:2282255/1‑96 (MQ=255)
                                                                    cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                            >  2:944564/1‑100 (MQ=255)
                                                                     ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                           >  1:610666/1‑100 (MQ=255)
                                                                     ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                          >  1:934117/1‑100 (MQ=255)
                                                                        cGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                        >  2:1939190/1‑100 (MQ=255)
                                                                         gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                       >  1:90773/1‑100 (MQ=255)
                                                                         gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                       >  1:625332/1‑100 (MQ=255)
                                                                          cccgactggggAGGGGGCGGGTCGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                             <  1:2583137/86‑1 (MQ=255)
                                                                          ccccggaggggCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACCTc                      >  1:1223776/8‑100 (MQ=255)
                                                                          ccccGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                          <  1:2445145/96‑1 (MQ=255)
                                                                           cccccgggggggggggCGGGGGGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                     <  1:1875515/94‑1 (MQ=255)
                                                                            ccccgcgggAGGGGGCCGGGGTCGGGGTGAAAAGGCCGCACTTCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                     <  1:2078513/94‑1 (MQ=255)
                                                                            ccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                    >  1:402264/1‑100 (MQ=255)
                                                                             ccccGGGGAGGGTGCGGGGGGCGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                    <  1:139552/97‑1 (MQ=255)
                                                                              tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                  >  1:2638336/1‑100 (MQ=255)
                                                                                   ggAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc             <  1:568275/100‑1 (MQ=255)
                                                                                   ggAGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc             <  1:238411/100‑1 (MQ=255)
                                                                                    gggggggCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg            <  1:1951967/96‑1 (MQ=255)
                                                                                       agGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc         >  1:426630/1‑100 (MQ=255)
                                                                                        gggCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa           <  2:200349/97‑1 (MQ=255)
                                                                                        gggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc        >  2:2524803/1‑100 (MQ=255)
                                                                                         ggCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa           <  2:295984/96‑1 (MQ=255)
                                                                                          gCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg      >  1:1989224/1‑100 (MQ=255)
                                                                                          gCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg      >  1:2368483/1‑100 (MQ=255)
                                                                                          gCCGGGGGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg      <  1:777548/100‑1 (MQ=255)
                                                                                            cGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaa    >  1:2436147/1‑100 (MQ=255)
                                                                                              gggTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt  <  2:731628/100‑1 (MQ=255)
                                                                                                    |                                                                                               
GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAT  >  NC_000913/3781053‑3781246

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