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

GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAT  >  NC_000913/3781053‑3781246
                                                                                                    |                                                                                               
gggAGAGGAGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑Ga                                                                                                 >  2:238409/1‑100 (MQ=255)
  gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAgg                                                                                               >  1:81515/1‑100 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAgggg                                                                                             <  1:2482769/100‑1 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAgggg                                                                                             <  2:2673560/100‑1 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAgggg                                                                                             <  1:1138662/100‑1 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑aa                                                                                          >  1:257849/1‑100 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑aa                                                                                          <  2:2045527/100‑1 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑aa                                                                                          >  1:200350/1‑100 (MQ=255)
       gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑aaa                                                                                         <  1:2088184/100‑1 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑aaaa                                                                                        >  1:2133575/1‑100 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAg                                                                                       <  1:1741815/100‑1 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAg                                                                                       <  1:1741938/100‑1 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGc                                                                                     <  1:192232/100‑1 (MQ=255)
            cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAg                                                                                       <  1:2198222/97‑1 (MQ=255)
            cGTTCGGCCCTGTATGTACTCCCTCGCTCTGAAACGACACCGCACGCTTTTTTTCTCCCGCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGcc                                                                                    <  2:607765/100‑1 (MQ=38)
              ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGc                                                                                  >  2:423217/1‑100 (MQ=255)
                ccgcgcTGTATGTACTCCCTCACTCTGAAACCACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCAc                                                                                >  1:1813878/3‑100 (MQ=255)
                 ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACt                                                                               >  1:306824/1‑100 (MQ=255)
                 ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACt                                                                               >  2:777541/1‑100 (MQ=255)
                   cgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGc                                                                             >  2:87275/1‑100 (MQ=255)
                    ncTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCTTCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCt                                                                            >  2:1973384/2‑100 (MQ=255)
                       gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTTTTTTCTCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCAc                                                                         >  2:1875493/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAAt                                                                     >  2:2687268/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAAt                                                                     >  2:931777/1‑100 (MQ=255)
                             ccgcccTCACTCGAAAACAACCCCCCTCTTTTTTTTTCCCCCCCCCCCGGGGGGGGGGGGGGCGGGGGT‑GGGGGA‑AAAAGGCCGCACTGCTCAcc                                                                        <  1:1768350/92‑1 (MQ=255)
                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAAttt                                                                   >  2:1592138/1‑100 (MQ=255)
                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAAttt                                                                   >  1:289573/1‑100 (MQ=255)
                                ccctccctCGAAAGCGACCCCTCTCTTTTTTTTTCCCCCGCGCCCCTGCGGGGGGGGGCCCGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTc                                                                 <  1:669988/99‑1 (MQ=255)
                                ccccccccctGAAGCGACCCCGCCTTTTTTTTTTCCCCCCCGCCCC‑CCGGGGGGGGGGGCGGGGTGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                                <  1:1114753/92‑1 (MQ=255)
                                cccTCCCCCTGAAACGACACCGCTCTTTTTTTTTCCCCCTCGCCCCCCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAAttt                                                                   <  2:1554204/97‑1 (MQ=37)
                                cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                                >  1:853067/1‑100 (MQ=255)
                                 cccccctctGAACCGCCCCGGCCCTCTTTTTTTCCCCCTCCCCCCTCCGGGGAGGGGGCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                               <  2:1016746/95‑1 (MQ=255)
                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                               >  2:337479/1‑100 (MQ=255)
                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                               >  1:2104139/1‑100 (MQ=255)
                                 ccTCACTCTGAAACGACACCACACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCGGGGT‑GAGGAT‑AAAAGGCCGCACTGCTCACCAATTTTc                                                                 >  2:1194264/1‑98 (MQ=255)
                                 ccTCACTCGAAAACGACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                               <  1:1932069/100‑1 (MQ=255)
                                   tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                             >  2:1711085/1‑100 (MQ=255)
                                     cctcgGAACCGACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                              <  2:2018216/96‑1 (MQ=255)
                                       tctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                         >  1:518909/1‑100 (MQ=255)
                                       tataaaaCTACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCCGGGGG‑GAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                         <  1:1105424/96‑1 (MQ=255)
                                        ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                        >  2:1968622/1‑100 (MQ=255)
                                          gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                      >  2:395785/1‑100 (MQ=255)
                                           aaaCCACCCCGCACTTTTTTTTTCCCCCCCGCCCCCCCGGGGAGAGGGCCGTGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                       <  2:1495741/98‑1 (MQ=255)
                                                acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCGGGGT‑GAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                        >  2:1704475/1‑92 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                             >  1:2487016/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                             >  1:980923/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                             >  1:508464/1‑100 (MQ=255)
                                                      cACTCTTTTTTCCCCCCTCGCCCCCCGGGGGGGGGGGCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                              <  2:1406017/96‑1 (MQ=255)
                                                        ctctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                        >  1:1175043/1‑100 (MQ=255)
                                                        ctctTTTTTTCTCCCCCCCCCCCCCGGGGGGGGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                         <  2:1433862/99‑1 (MQ=255)
                                                           tttttttCCCCCCCCCCCCCCCGGGGGGGGGGCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                        <  1:2157465/97‑1 (MQ=255)
                                                           tttttttCCCCCCCCCCCCCCCGGGGGGGGGCCCGGGGT‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                              <  2:2064707/91‑1 (MQ=255)
                                                           ttttttcCGCCCCCCCCCCCCCGGGGGGAGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                      <  2:1871177/99‑1 (MQ=255)
                                                            ttttttCTCCTCCGCCCCGCGGGGGGGGGGGCCGGGGT‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                        <  2:124733/96‑1 (MQ=255)
                                                             tttttCTCCCCCGCCCCCCCGGGGGGAGGGCCGGGGT‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                   <  2:1380291/100‑1 (MQ=255)
                                                               tttCCCCCCCGCCCCTCCGGGGAGGGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAgga                                 <  2:306821/100‑3 (MQ=255)
                                                               tcccccccTCCCCCCGCGGGGGAG‑GGGGCGGGGT‑GGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                                 <  1:2130933/95‑1 (MQ=255)
                                                                tcccccccccccccGCCGGGGGGAGGGCCGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                <  1:605424/90‑1 (MQ=255)
                                                                  ctcCCTCGCCCCCCCGGGGGGAGGGCCGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                              <  2:2069749/100‑1 (MQ=255)
                                                                   tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                             >  1:2458212/1‑100 (MQ=255)
                                                                   tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                             >  1:2413242/1‑100 (MQ=255)
                                                                    cccTCTCTCCTCCGGGGAGAGGGCCGGGGG‑GAGGGG‑AAAAGGCCGCACTGCTCACGAATTTTCTGCTATCCTGCCCCCGCACTAAGGTGGCAGTAttc                              >  2:2282255/1‑96 (MQ=255)
                                                                    cccTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                            >  2:944564/1‑100 (MQ=255)
                                                                     ccTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                           >  1:610666/1‑100 (MQ=255)
                                                                     ccTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                          >  1:934117/1‑100 (MQ=255)
                                                                        cGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                        >  2:1939190/1‑100 (MQ=255)
                                                                         gCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                       >  1:90773/1‑100 (MQ=255)
                                                                         gCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                       >  1:625332/1‑100 (MQ=255)
                                                                          cccgactggggAGGGGGCGGGTCG‑GAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                             <  1:2583137/86‑1 (MQ=255)
                                                                          ccccggaggggCGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACCTc                      >  1:1223776/8‑100 (MQ=255)
                                                                          ccccGCGGGGGGGGGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                          <  1:2445145/96‑1 (MQ=255)
                                                                           cccccgggggggggggCGGGGGG‑GAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                     <  1:1875515/94‑1 (MQ=255)
                                                                            ccccgcgggAGGGGGCCGGGGTCGGGGTG‑AAAAGGCCGCACTTCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                     <  1:2078513/94‑1 (MQ=255)
                                                                            ccTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                    >  1:402264/1‑100 (MQ=255)
                                                                             ccccGGGGAGGGTGCGGGGGG‑CGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                    <  1:139552/97‑1 (MQ=255)
                                                                              tCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                  >  1:2638336/1‑100 (MQ=255)
                                                                                   ggAGGGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc             <  1:568275/100‑1 (MQ=255)
                                                                                   ggAGGGGGCGGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc             <  1:238411/100‑1 (MQ=255)
                                                                                    gggggggCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg            <  1:1951967/96‑1 (MQ=255)
                                                                                       agGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc         >  1:426630/1‑100 (MQ=255)
                                                                                        gggCCGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa           <  2:200349/97‑1 (MQ=255)
                                                                                        gggCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc        >  2:2524803/1‑100 (MQ=255)
                                                                                         ggCCGGGGG‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa           <  2:295984/96‑1 (MQ=255)
                                                                                          gCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg      >  1:1989224/1‑100 (MQ=255)
                                                                                          gCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg      >  1:2368483/1‑100 (MQ=255)
                                                                                          gCCGGGGG‑GAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg      <  1:777548/100‑1 (MQ=255)
                                                                                            cGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaa    >  1:2436147/1‑100 (MQ=255)
                                                                                              gggT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt  <  2:731628/100‑1 (MQ=255)
                                                                                                    |                                                                                               
GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAT  >  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: ‑