BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,140	0	A	G	25.0%	207.5 / 28.1	104	intergenic (+123/‑75)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base A (54/24); new base G (0/26); total (54/50)
Fisher's exact test for biased strand distribution p-value = 5.39e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 113 total aligned reads displayed.

CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAC  >  NC_000913/3781041‑3781239

cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGgaga                                                                                                        <  2:2561843/100‑1 (MQ=255)

 ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGagag                                                                                                       <  2:718195/100‑1 (MQ=255)

  ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAgg                                                                                                      <  2:1725844/100‑1 (MQ=255)

    ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGc                                                                                                    <  1:3369442/100‑1 (MQ=255)

      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCTGGGAGAGGGCC‑g                                                                                                 >  2:2163764/1‑100 (MQ=255)

      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑g                                                                                                 >  1:3141025/1‑100 (MQ=255)

      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑g                                                                                                 >  1:2791800/1‑100 (MQ=255)

      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑g                                                                                                 <  1:1905586/100‑1 (MQ=255)

       ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑gg                                                                                                >  1:1027559/1‑100 (MQ=255)

        cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑ggg                                                                                               >  2:3322391/1‑100 (MQ=255)

        cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑ggg                                                                                               <  1:1805343/100‑1 (MQ=255)

         tCAGGGAGAGGGGACCTTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑gggg                                                                                              <  2:1991843/100‑1 (MQ=255)

         tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑gggg                                                                                              <  1:1574454/100‑1 (MQ=255)

         tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑gggg                                                                                              <  1:47987/100‑1 (MQ=255)

          cAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGt                                                                                             >  1:1240665/1‑100 (MQ=255)

           aGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTg                                                                                            <  2:3194865/100‑1 (MQ=255)

           aGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTg                                                                                            <  2:1890654/100‑1 (MQ=255)

               agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAggg                                                                                        <  1:2525156/100‑1 (MQ=255)

                gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAggg‑g                                                                                      <  2:32771/100‑1 (MQ=255)

                gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAggg‑g                                                                                      >  2:1547431/1‑100 (MQ=255)

                    ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑aaaa                                                                                 <  2:3184565/100‑1 (MQ=255)

                    ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑aaaa                                                                                 >  1:2189997/1‑100 (MQ=255)

                     gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGCGAGGGCC‑AGGGTGAGGG‑G‑AAAAg                                                                                >  1:2830124/1‑100 (MQ=255)

                     gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAg                                                                                <  1:3111586/100‑1 (MQ=255)

                     gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAg                                                                                >  2:2729489/1‑100 (MQ=255)

                     gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAg                                                                                >  1:1697859/1‑100 (MQ=255)

                       ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGc                                                                              >  2:45768/1‑100 (MQ=255)

                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAggga                                                                             >  2:587992/1‑97 (MQ=255)

                           tCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCa                                                                          >  2:833720/1‑100 (MQ=255)

                             ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACt                                                                        >  1:721934/1‑100 (MQ=255)

                                gcTGTATTTCCCCCCCCCCCCTAAACCGACCCCGCACTC‑TTTTTTCCCC‑CCCCCCCCCCCGGGGGGGGGGCCGGGGGGGGGGG‑G‑AAAAGGCCGc                                                                           <  2:1222325/94‑1 (MQ=255)

                                gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCt                                                                     >  1:1410975/1‑100 (MQ=255)

                                 cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAACGGCCGCACTGCTc                                                                    >  2:1953272/1‑100 (MQ=255)

                                  tgtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCa                                                                   >  2:1369640/1‑100 (MQ=255)

                                     ttccgatctcTCACTCGGAAACGACCCCGCACTCTTTTTTTCTCC‑CGCGCCCCCCCGGGGGGGGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTc                                                                    <  2:2300780/87‑1 (MQ=255)

                                     atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCa                                                                >  1:1705791/1‑100 (MQ=255)

                                     atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCa                                                                >  1:1706318/1‑100 (MQ=255)

                                      tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCaa                                                               >  2:2248076/1‑100 (MQ=255)

                                       gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAAt                                                              >  2:1864433/1‑100 (MQ=255)

                                             ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                        >  1:3183117/1‑100 (MQ=255)

                                             ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                        >  2:1347117/1‑100 (MQ=255)

                                              ccccctctAAAACGACCCCCCCCTCTTTTTTCCCCC‑CCCCCCCCCCCGGGGAGGGGGCG‑GGGGGGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                       <  1:480527/96‑1 (MQ=255)

                                              cTCCCTCTAAAGCGACACCGCACTCTTTTTTTCCCC‑CCCGCCCCCCCGGGGAGAGGCCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                       <  2:1326902/100‑1 (MQ=255)

                                              cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑CAAAGGCCGCCCTGCTCACCAATTTTCTGc                                                       >  2:3077151/1‑100 (MQ=255)

                                               tCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                      >  1:2265495/1‑100 (MQ=255)

                                                cACTCTAAAACGACACCGCACTCTTTTTTTCGCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑A‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                                     <  2:1068236/100‑1 (MQ=255)

                                                   tctGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCTGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                  <  1:2163769/100‑1 (MQ=255)

                                                    ctGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                 >  1:2005585/1‑100 (MQ=255)

                                                    ctGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                 >  1:1336380/1‑100 (MQ=255)

                                                        aaCGACCCCGCCCTCTTTTTTCCCCC‑CTCGCCCCCCGGGGGAGGGGCCG‑GGGGTGGGGG‑A‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                             <  2:2490489/100‑1 (MQ=255)

                                                          cGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                           >  1:1333874/1‑100 (MQ=255)

                                                           gACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                          >  1:2705697/1‑100 (MQ=255)

                                                           gACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                          >  1:2461463/1‑100 (MQ=255)

                                                              nccGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                       >  2:1101061/2‑100 (MQ=255)

                                                               cgccACTCTTTTTTTCTCC‑CGCCCCCGCCGGGGGGGGGGGCC‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                  <  2:1240662/85‑1 (MQ=255)

                                                               ccGCCCTTTTTTTTTCCCC‑CCCGCCCCCCGGGGGGGGGGGCG‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                       <  2:1336377/99‑1 (MQ=255)

                                                               ccGCCCTCTTTTTTTCCCC‑CTCGCCCCCCGGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                       <  2:2265489/99‑1 (MQ=38)

                                                                  cgctctTTTTTTCCCC‑CTCTCCCCCCCGGGGAGATGCCG‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                              <  2:1306190/87‑1 (MQ=255)

                                                                  cACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                   >  1:1892018/1‑100 (MQ=255)

                                                                   ccgttttttgtCTCCGCTCCCCCCCCGGAGGGGAGGGCG‑GGGGTGGGGG‑A‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGTCCCCGCAc                                       <  1:2175069/92‑1 (MQ=255)

                                                                   ccctcttttttCCCC‑CCCCCCCCTCCGGGGGGGGGGGC‑GGGGGGGGGG‑GAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                    <  1:3322404/94‑1 (MQ=255)

                                                                     tctTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                >  2:2793058/1‑100 (MQ=255)

                                                                      ctTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                               >  2:1150833/1‑100 (MQ=255)

                                                                      ctTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                               >  2:3149337/1‑100 (MQ=255)

                                                                       tttttttCTCC‑CCCCCCCCGCCGGGGGGGGGGCC‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                              <  1:2466641/100‑1 (MQ=255)

                                                                       tttttttC‑CC‑CCCCCCCCCCCGGGGGGGGGGCCGGGGGTGAGGG‑A‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                 <  2:3189450/97‑1 (MQ=255)

                                                                        ttttttCCCC‑CGGGCCCCCCCGGGGGGGGGGCG‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                              <  2:2791790/99‑1 (MQ=255)

                                                                        ttttttCACC‑CAACCCCCCCGGGAGAGGCGCCC‑GGGGTGGGGG‑A‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                              <  2:2203588/83‑1 (MQ=255)

                                                                         tttttCCCC‑CCCCCCCCCCGGGGGGGGGGGCG‑GGGGGGGGGG‑A‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                            <  1:435293/100‑1 (MQ=255)

                                                                          ttttCCCC‑CCCCCCCCGCGGGGGGGGGGGCG‑GGGGTGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                           <  1:1876722/100‑1 (MQ=255)

                                                                           tttCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                          >  1:1704695/1‑100 (MQ=255)

                                                                           tttCTCC‑CTCCCCCTTCGGGGGAGGGGGCG‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                          <  1:1369643/100‑1 (MQ=255)

                                                                            ttCTCC‑CCCGCCCCCCGGGGGGGGGGGCG‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                         <  1:1101063/100‑1 (MQ=255)

                                                                            tccccc‑cTCCCCCCCGGGGGGGGGGGGCG‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                          <  1:3149349/95‑1 (MQ=255)

                                                                             tctcC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggagatcggaag                        >  2:736080/1‑90 (MQ=255)

                                                                              ctcC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                       >  2:582918/1‑100 (MQ=255)

                                                                              cccc‑cTCGCCCCCCGGGGGAGGGGGCG‑GGGGTGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                       <  1:1547434/98‑1 (MQ=255)

                                                                               ccccccccccccccCGGGGGAGGGGCC‑GGGGGGGGGG‑GAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                        <  1:2729498/92‑1 (MQ=255)

                                                                                cc‑cTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                     >  1:2577297/1‑100 (MQ=255)

                                                                                cc‑cTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                     >  1:33360/1‑100 (MQ=255)

                                                                                cc‑cTCGCCCCCCGGGGGAGGGGGCC‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                        <  2:1755970/97‑1 (MQ=255)

                                                                                 c‑cTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                    >  2:2577672/1‑100 (MQ=255)

                                                                                   cTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                   >  1:1721217/1‑100 (MQ=255)

                                                                                   cTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                   >  2:1721530/1‑100 (MQ=255)

                                                                                    cgcccccTCCGGGGGGGAGGCC‑GGGGGGGGGGTT‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                   <  1:2648896/97‑1 (MQ=255)

                                                                                     cGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                 >  2:2865503/1‑100 (MQ=255)

                                                                                       ccccggaggggCGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc               >  2:845701/8‑100 (MQ=255)

                                                                                        cccTCGGGGGGGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg              <  1:579783/100‑1 (MQ=255)

                                                                                        cccTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg              >  1:676210/1‑100 (MQ=255)

                                                                                         ccccgggggggggggCG‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg              <  1:294077/94‑1 (MQ=255)

                                                                                          cTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa            >  1:607656/1‑100 (MQ=255)

                                                                                          cTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGAACTAAGGGGGCAGTATGCTAAACATCGTa            >  2:1361687/1‑100 (MQ=255)

                                                                                           ccgggggggggggCG‑GGGGTGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc           <  1:149017/97‑1 (MQ=255)

                                                                                            cgggggAGAGGGCG‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc           <  2:3058551/97‑1 (MQ=255)

                                                                                            ccGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt          >  2:1576419/1‑100 (MQ=255)

                                                                                              ggggAGGGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt             <  2:1892014/95‑1 (MQ=255)

                                                                                              ggggAGGGGCCG‑GGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt        <  1:350548/100‑1 (MQ=255)

                                                                                              ggggAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt        >  1:3071619/1‑100 (MQ=255)

                                                                                                 gcgagGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg     <  1:845701/98‑1 (MQ=255)

                                                                                                    agGGGC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc  <  2:360609/100‑1 (MQ=255)

CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCC‑CTCGCCCCTCCGGGGAGAGGGCC‑GGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAC  >  NC_000913/3781041‑3781239

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