Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1400AG20.0% 177.4 / 10.8 82intergenic (+123/‑75)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (50/14);  new base G (0/16);  total (50/32)
Fisher's exact test for biased strand distribution p-value = 5.39e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAA  >  NC_000913/3781041‑3781238
                                                                                                   |                                                                                                    
cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                                      <  2:567730/100‑1 (MQ=255)
    ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                                  >  2:2015089/1‑100 (MQ=255)
     tcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                                                 <  1:1460591/100‑1 (MQ=255)
      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑g                                                                                               >  2:1826496/1‑100 (MQ=255)
      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑g                                                                                               >  2:2255902/1‑100 (MQ=255)
      cccTCAGGGAGAGGGGACCGTGCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACGGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                     <  2:2174850/95‑1 (MQ=255)
         tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑gggg                                                                                            <  2:71257/100‑1 (MQ=255)
           aGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTg                                                                                          <  2:30668/100‑1 (MQ=255)
            gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGa                                                                                         >  1:1152569/1‑100 (MQ=255)
            gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGa                                                                                         >  1:1066072/1‑100 (MQ=255)
            gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGa                                                                                         >  2:1055004/1‑100 (MQ=255)
               agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAggg                                                                                      <  2:1951707/100‑1 (MQ=255)
                  ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑aa                                                                                  >  2:2136249/1‑100 (MQ=255)
                    ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑aaaa                                                                                >  1:840052/1‑100 (MQ=255)
                       ccGTTCGGCGCTGTATGTTCTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGc                                                                             <  1:1826508/100‑1 (MQ=255)
                       ccGTTCGGCGCTGTATGTCCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGc                                                                             <  1:844428/100‑1 (MQ=255)
                       ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGc                                                                             <  1:1052518/100‑1 (MQ=255)
                        cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGcc                                                                            >  2:2129567/1‑100 (MQ=255)
                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCg                                                                           >  2:746556/1‑100 (MQ=255)
                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCg                                                                           >  1:1066082/1‑100 (MQ=255)
                            cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCAc                                                                        >  2:1368981/1‑100 (MQ=255)
                            cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCAc                                                                        >  2:1369392/1‑100 (MQ=255)
                              gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTg                                                                      >  1:1094928/1‑100 (MQ=255)
                                gcTGCATGTCCTCCCCCCCTCTGAACCAACACCCCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCC‑GGGGAGAGAGA‑AAAAGGCCGCGCTGc                                                                     <  2:1634964/99‑1 (MQ=255)
                                 cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTc                                                                   >  1:683745/1‑100 (MQ=255)
                                   gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTTCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCAc                                                                 >  2:656104/1‑100 (MQ=255)
                                   gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCAc                                                                 >  1:1171862/1‑100 (MQ=255)
                                     atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCa                                                               >  1:495791/1‑100 (MQ=255)
                                     atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCa                                                               >  1:61317/1‑100 (MQ=255)
                                      tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGTGTGAGGGG‑AAAAGGCCGCACTGCTCACCaa                                                              >  2:1255364/1‑100 (MQ=255)
                                      tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCaa                                                              >  1:764555/1‑100 (MQ=255)
                                       gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAAt                                                             >  2:2114585/1‑100 (MQ=255)
                                          cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAAtttt                                                          >  1:458676/1‑100 (MQ=255)
                                          cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAAtttt                                                          >  1:1807038/1‑100 (MQ=255)
                                            cctccACCCTGAACCGACGCCGGACTTGTTTTTCCCCCCCAGCGCCGCCGGGGAGGGGGCC‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTc                                                         <  2:144065/95‑1 (MQ=255)
                                            cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCTCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  2:944006/1‑100 (MQ=255)
                                            cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  2:1105174/1‑100 (MQ=255)
                                            cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  1:1946670/1‑100 (MQ=255)
                                            cccTCACTCTGAAACGACACCGAACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  2:1314345/1‑100 (MQ=255)
                                             ccacactgggAAACGACCCCGCTCTCTTTTTTTCCCCTACACCCGTAGGGGGAGTGGTCC‑GGGGGGAGGGG‑AAAAGGCCGCACTGCGCACCAATTTTCTg                                                       <  1:410363/91‑1 (MQ=255)
                                             ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCGCTCTCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                       >  1:1587807/1‑100 (MQ=255)
                                              ccccctggaaaaCGACCCCCCCCT‑TTTTTTTCCCCCCCCCCCCCCCGGGGGGGGGGGC‑GGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                      <  1:1723584/91‑1 (MQ=255)
                                              cTCACTCTGAAACGACACCGCCCTCTTTTTTTCCCCCTGGCCCCTCCGGGGAGAGGGCG‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                      <  2:495788/100‑1 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:1308497/1‑100 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:1460728/1‑100 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:639826/1‑100 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGCGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:1102601/1‑100 (MQ=255)
                                                    ctGAAACGACACCGCACTCTTTTTTTTTTTTTTTCACCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑ATAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                >  2:804727/1‑100 (MQ=255)
                                                    ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATTCt                                                >  1:733135/1‑100 (MQ=255)
                                                       aagcGACACCGCACTCTTTTTTTCCCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                             <  2:438752/97‑1 (MQ=255)
                                                          cccccccGCACTTTTTTTTCCCCCCCCCCCCCTCGGGGGAGGGGGCG‑GGGGGGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                           <  1:2129582/94‑1 (MQ=255)
                                                             ccccGCCCTCTTTTTTCCTCCCCCGCCCCGCGGGGGGGGGGGCG‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                        <  2:1603970/97‑1 (MQ=255)
                                                             cacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                                       >  2:873691/1‑100 (MQ=255)
                                                              cccGCTCTCTTTTTTTCCCCCCCGCCCCTGCGGGGGAGGGGCC‑GGGTGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                         <  1:1314354/96‑1 (MQ=255)
                                                              cccGCCCTCTTTTTTTCTCCCCCGCCCCTCCGGGGAGGGGCCC‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                      <  1:2017037/99‑1 (MQ=255)
                                                               ccgctttttttttttCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                             <  1:1127736/85‑1 (MQ=255)
                                                               ccGCACTC‑TTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                    >  1:1076010/1‑100 (MQ=255)
                                                                cGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                    >  1:1253197/1‑100 (MQ=255)
                                                                 gCACTCTTTTTTTTTTTTTCTCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                   >  2:1571270/1‑100 (MQ=255)
                                                                  cgttgttttttgCCCCCGCCCCCCTCGGGGGGGGGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                      <  2:1308488/91‑1 (MQ=255)
                                                                    ctctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                >  2:377381/1‑100 (MQ=255)
                                                                       tttttttCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                             >  1:1716936/1‑100 (MQ=255)
                                                                       tttttttCCCCCCCCCCCCTCCGGGGGGGGGGCC‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                <  2:61315/97‑1 (MQ=255)
                                                                        ttttttCCCCCCCCCCCCCCCGGGGGAGGGGCC‑GGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                             <  1:746563/100‑1 (MQ=255)
                                                                        ctcttccgatcTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                            <  1:1803910/90‑1 (MQ=255)
                                                                         tttttCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                           >  2:1778398/1‑100 (MQ=255)
                                                                          ttttCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                          >  2:2127257/1‑100 (MQ=255)
                                                                              ctcCCTCGCCCCCCGGGGGAGAGGGCC‑GGGGGGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                      <  2:1094921/100‑1 (MQ=255)
                                                                              ccccccccccccccGGGGGAGGGGGCG‑GGGTCGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                      <  1:2255917/92‑1 (MQ=255)
                                                                                ccccccccccTCCGGGGAGGGGGCC‑GGGGTGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    <  2:1152562/94‑1 (MQ=255)
                                                                                cccTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCATTTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    >  2:1587898/1‑100 (MQ=255)
                                                                                cccTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    >  1:612635/1‑100 (MQ=255)
                                                                                cccTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    >  2:2089839/1‑100 (MQ=255)
                                                                                 ccTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                   >  1:1354663/1‑100 (MQ=255)
                                                                                 ccTCCCCCCCCGGGGGAGGGGCCG‑GGGGTGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                   <  1:1368990/100‑1 (MQ=255)
                                                                                  cTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAgatcggaaga                  >  2:1477648/1‑91 (MQ=255)
                                                                                    cgcgccgcgggggggCGGGCG‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                 <  1:2254686/90‑1 (MQ=255)
                                                                                    cGCCCCCCGGGGGGGGGGGCC‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                  <  2:111274/98‑1 (MQ=255)
                                                                                     gCCCCTGCGGGGGGGGGGCC‑GGGGGGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    <  2:1295481/95‑1 (MQ=255)
                                                                                     ctcccTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt               <  2:1047580/98‑1 (MQ=255)
                                                                                       cccTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAtag             >  1:1047587/1‑98 (MQ=255)
                                                                                            cGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt        >  1:1304543/1‑100 (MQ=255)
                                                                                             ggggggAGGGCG‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc          <  2:1102594/97‑1 (MQ=255)
                                                                                             ggggAGGGGGCG‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt       <  2:1076003/100‑1 (MQ=255)
                                                                                                 cgagGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg    >  1:772365/2‑99 (MQ=255)
                                                                                                  gagGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa  >  2:421740/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAA  >  NC_000913/3781041‑3781238

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: ‑