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

CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAA  >  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)
      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                               >  2:1826496/1‑100 (MQ=255)
      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                               >  2:2255902/1‑100 (MQ=255)
      cccTCAGGGAGAGGGGACCGTGCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACGGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                     <  2:2174850/95‑1 (MQ=255)
         tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                                            <  2:71257/100‑1 (MQ=255)
           aGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                          <  2:30668/100‑1 (MQ=255)
            gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                         >  1:1152569/1‑100 (MQ=255)
            gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                         >  1:1066072/1‑100 (MQ=255)
            gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                         >  2:1055004/1‑100 (MQ=255)
               agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                      <  2:1951707/100‑1 (MQ=255)
                  ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                  >  2:2136249/1‑100 (MQ=255)
                    ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                >  1:840052/1‑100 (MQ=255)
                       ccGTTCGGCGCTGTATGTTCTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                             <  1:1826508/100‑1 (MQ=255)
                       ccGTTCGGCGCTGTATGTCCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                             <  1:844428/100‑1 (MQ=255)
                       ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                             <  1:1052518/100‑1 (MQ=255)
                        cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                                            >  2:2129567/1‑100 (MQ=255)
                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                                                           >  2:746556/1‑100 (MQ=255)
                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                                                           >  1:1066082/1‑100 (MQ=255)
                            cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                        >  2:1368981/1‑100 (MQ=255)
                            cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                        >  2:1369392/1‑100 (MQ=255)
                              gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTg                                                                      >  1:1094928/1‑100 (MQ=255)
                                gcTGCATGTCCTCCCCCCCTCTGAACCAACACCCCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCGGGGAGAGAGAAAAAGGCCGCGCTGc                                                                     <  2:1634964/99‑1 (MQ=255)
                                 cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                   >  1:683745/1‑100 (MQ=255)
                                   gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTTCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                 >  2:656104/1‑100 (MQ=255)
                                   gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                 >  1:1171862/1‑100 (MQ=255)
                                     atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                                               >  1:495791/1‑100 (MQ=255)
                                     atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                                               >  1:61317/1‑100 (MQ=255)
                                      tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGTGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                              >  2:1255364/1‑100 (MQ=255)
                                      tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                              >  1:764555/1‑100 (MQ=255)
                                       gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                             >  2:2114585/1‑100 (MQ=255)
                                          cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                          >  1:458676/1‑100 (MQ=255)
                                          cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                          >  1:1807038/1‑100 (MQ=255)
                                            cctccACCCTGAACCGACGCCGGACTTGTTTTTCCCCCCCAGCGCCGCCGGGGAGGGGGCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                         <  2:144065/95‑1 (MQ=255)
                                            cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCTCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  2:944006/1‑100 (MQ=255)
                                            cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  2:1105174/1‑100 (MQ=255)
                                            cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  1:1946670/1‑100 (MQ=255)
                                            cccTCACTCTGAAACGACACCGAACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                        >  2:1314345/1‑100 (MQ=255)
                                             ccacactgggAAACGACCCCGCTCTCTTTTTTTCCCCTACACCCGTAGGGGGAGTGGTCCGGGGGGAGGGGAAAAGGCCGCACTGCGCACCAATTTTCTg                                                       <  1:410363/91‑1 (MQ=255)
                                             ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCGCTCTCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                       >  1:1587807/1‑100 (MQ=255)
                                              ccccctggaaaaCGACCCCCCCCTTTTTTTTCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                      <  1:1723584/91‑1 (MQ=255)
                                              cTCACTCTGAAACGACACCGCCCTCTTTTTTTCCCCCTGGCCCCTCCGGGGAGAGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                      <  2:495788/100‑1 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:1308497/1‑100 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:1460728/1‑100 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:639826/1‑100 (MQ=255)
                                                  ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGCGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                  >  1:1102601/1‑100 (MQ=255)
                                                    ctGAAACGACACCGCACTCTTTTTTTTTTTTTTTCACCTCCGGGGAGAGGGCCGGGGTGAGGGGATAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                >  2:804727/1‑100 (MQ=255)
                                                    ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATTCt                                                >  1:733135/1‑100 (MQ=255)
                                                       aagcGACACCGCACTCTTTTTTTCCCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                             <  2:438752/97‑1 (MQ=255)
                                                          cccccccGCACTTTTTTTTCCCCCCCCCCCCCTCGGGGGAGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                           <  1:2129582/94‑1 (MQ=255)
                                                             ccccGCCCTCTTTTTTCCTCCCCCGCCCCGCGGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                        <  2:1603970/97‑1 (MQ=255)
                                                             cacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                                       >  2:873691/1‑100 (MQ=255)
                                                              cccGCTCTCTTTTTTTCCCCCCCGCCCCTGCGGGGGAGGGGCCGGGTGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                         <  1:1314354/96‑1 (MQ=255)
                                                              cccGCCCTCTTTTTTTCTCCCCCGCCCCTCCGGGGAGGGGCCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                      <  1:2017037/99‑1 (MQ=255)
                                                               ccgctttttttttttCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                             <  1:1127736/85‑1 (MQ=255)
                                                               ccGCACTCTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                    >  1:1076010/1‑100 (MQ=255)
                                                                cGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                    >  1:1253197/1‑100 (MQ=255)
                                                                 gCACTCTTTTTTTTTTTTTCTCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                   >  2:1571270/1‑100 (MQ=255)
                                                                  cgttgttttttgCCCCCGCCCCCCTCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                      <  2:1308488/91‑1 (MQ=255)
                                                                    ctctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                >  2:377381/1‑100 (MQ=255)
                                                                       tttttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                             >  1:1716936/1‑100 (MQ=255)
                                                                       tttttttCCCCCCCCCCCCTCCGGGGGGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                <  2:61315/97‑1 (MQ=255)
                                                                        ttttttCCCCCCCCCCCCCCCGGGGGAGGGGCCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                             <  1:746563/100‑1 (MQ=255)
                                                                        ctcttccgatcTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                            <  1:1803910/90‑1 (MQ=255)
                                                                         tttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                           >  2:1778398/1‑100 (MQ=255)
                                                                          ttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                          >  2:2127257/1‑100 (MQ=255)
                                                                              ctcCCTCGCCCCCCGGGGGAGAGGGCCGGGGGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                      <  2:1094921/100‑1 (MQ=255)
                                                                              ccccccccccccccGGGGGAGGGGGCGGGGTCGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                      <  1:2255917/92‑1 (MQ=255)
                                                                                ccccccccccTCCGGGGAGGGGGCCGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    <  2:1152562/94‑1 (MQ=255)
                                                                                cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCATTTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    >  2:1587898/1‑100 (MQ=255)
                                                                                cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    >  1:612635/1‑100 (MQ=255)
                                                                                cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    >  2:2089839/1‑100 (MQ=255)
                                                                                 ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                   >  1:1354663/1‑100 (MQ=255)
                                                                                 ccTCCCCCCCCGGGGGAGGGGCCGGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                   <  1:1368990/100‑1 (MQ=255)
                                                                                  cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAgatcggaaga                  >  2:1477648/1‑91 (MQ=255)
                                                                                    cgcgccgcgggggggCGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                 <  1:2254686/90‑1 (MQ=255)
                                                                                    cGCCCCCCGGGGGGGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                  <  2:111274/98‑1 (MQ=255)
                                                                                     gCCCCTGCGGGGGGGGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                    <  2:1295481/95‑1 (MQ=255)
                                                                                     ctcccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt               <  2:1047580/98‑1 (MQ=255)
                                                                                       cccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAtag             >  1:1047587/1‑98 (MQ=255)
                                                                                            cGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt        >  1:1304543/1‑100 (MQ=255)
                                                                                             ggggggAGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc          <  2:1102594/97‑1 (MQ=255)
                                                                                             ggggAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt       <  2:1076003/100‑1 (MQ=255)
                                                                                                 cgagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg    >  1:772365/2‑99 (MQ=255)
                                                                                                  gagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa  >  2:421740/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAA  >  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: