BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,211,610	A→G	15.3%	intergenic (+33/+70)	icdC → / ← iraM	pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,211,610	0	A	G	15.3%	371.7 / 54.2	156	intergenic (+33/+70)	icdC/iraM	pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor
Reads supporting (aligned to +/- strand): ref base A (63/69); new base G (14/10); total (77/79)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 166 total aligned reads displayed.

TCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTAT  >  NC_000913/1211511‑1211709

 cTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtat                                                                                                    >  2:1232560/1‑100 (MQ=255)

  tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                   >  2:423530/1‑100 (MQ=255)

  tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                   <  2:354660/100‑1 (MQ=255)

  tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                                                                                   >  2:3096231/1‑100 (MQ=255)

   gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGGGTTTGGTGAAGCGATCATCGAAAACATGTAATCTTTCCATGTGTTAATTATTGAAACGGGTGTATaa                                                                                                  >  2:2534681/1‑100 (MQ=21)

   gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATaa                                                                                                  >  1:2692533/1‑100 (MQ=255)

    aTGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAac                                                                                                 <  2:172704/100‑1 (MQ=255)

     tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                >  2:1352303/1‑100 (MQ=255)

      ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                                                                               <  1:2519178/100‑1 (MQ=255)

        aaGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGc                                                                                             >  2:979577/1‑100 (MQ=255)

          ggCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGccc                                                                                           >  2:620002/1‑100 (MQ=255)

           gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCg                                                                                          >  1:1633518/1‑100 (MQ=255)

            cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                                                                         >  1:3389724/1‑100 (MQ=255)

            cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                                                                         >  1:2864860/1‑100 (MQ=255)

             gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtt                                                                                        <  1:489000/100‑1 (MQ=25)

              cTAAGCTGCTGAAATTTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                       <  2:890217/100‑1 (MQ=255)

              cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                       <  2:1576700/100‑1 (MQ=255)

              cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                       <  1:758814/100‑1 (MQ=255)

              cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                       <  1:505826/100‑1 (MQ=255)

              cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                       <  1:3357339/100‑1 (MQ=255)

              cTAAGCTGCGGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                       <  2:2017942/100‑1 (MQ=255)

               tAATCTGCTGAAATGTTCAGAGTTTAGCGGAGTGCTAATCGAAAACATGTAATCACTCCATGTGTTAGATATTGAAACGTGGGTACAACACGTCCg                                                                                          <  2:753303/96‑1 (MQ=16)

               tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                      <  2:2271028/100‑3 (MQ=25)

               tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                      <  2:1649372/100‑3 (MQ=25)

               tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                      <  1:1080517/100‑3 (MQ=25)

               tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                      <  1:1924148/100‑3 (MQ=25)

               tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                                                                                      <  1:807572/100‑1 (MQ=255)

               tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATAGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgtttc                                                                                      <  2:2900757/100‑3 (MQ=25)

                 aGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                                                                                    <  1:1447114/100‑1 (MQ=255)

                  gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                   <  1:2540074/100‑1 (MQ=255)

                  gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                   >  2:2430238/1‑100 (MQ=255)

                  gctgctGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGtttt                                                                                   <  1:1328286/100‑1 (MQ=255)

                     gctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttatt                                                                                <  2:1339062/100‑1 (MQ=255)

                      ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                               >  1:2269996/1‑100 (MQ=255)

                        gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                             >  2:2075855/1‑100 (MQ=255)

                        gAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                             >  1:2054636/1‑100 (MQ=255)

                         aaaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtg                                                                            <  2:917693/100‑1 (MQ=255)

                           aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                          <  1:3045869/100‑1 (MQ=255)

                           aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                          <  1:1022542/100‑1 (MQ=255)

                           aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                          <  1:3528089/100‑1 (MQ=255)

                           aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAAAATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                          >  2:530696/1‑100 (MQ=255)

                             gTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGa                                                                        >  1:2479089/1‑100 (MQ=255)

                               tCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtgatat                                                                      >  2:1743116/1‑96 (MQ=255)

                                cAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGatat                                                                     >  2:936567/1‑100 (MQ=255)

                                     tttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATAttatta                                                                >  1:2757240/1‑100 (MQ=255)

                                      ttGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTaa                                                               >  1:1270817/1‑100 (MQ=255)

                                        ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAATTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                             <  2:2333498/100‑1 (MQ=255)

                                         gTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                            >  1:2353628/1‑100 (MQ=255)

                                           gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCa                                                          >  1:94088/1‑100 (MQ=255)

                                               cGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATc                                                      <  2:1997988/100‑1 (MQ=255)

                                                gATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCg                                                     <  2:2028350/100‑1 (MQ=255)

                                                   catcGAAAACATGTAATCTCTCCATTTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                  <  2:1525537/100‑1 (MQ=255)

                                                      atctaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                               <  2:1453385/96‑1 (MQ=255)

                                                      atctaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                               <  2:1454157/96‑1 (MQ=255)

                                                        aaaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATAtt                                             <  2:2531189/100‑1 (MQ=255)

                                                        aaaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATAtt                                             >  1:1559592/1‑100 (MQ=255)

                                                          aaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCataagat                                           >  1:1453393/1‑96 (MQ=255)

                                                             aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                        >  2:3385889/1‑100 (MQ=255)

                                                             aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                        >  2:2665072/1‑100 (MQ=255)

                                                             aTGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                        >  2:61300/1‑100 (MQ=255)

                                                               gacaTCTCTCCATGTGTTAAATATTGAAACGGGCTTATACCACGCCCGTTGTTTTATTTATGTGGATATGATTAACAGCATATCGGGCATATTTATATGa                                      <  2:672988/97‑1 (MQ=255)

                                                                 aaTCTCTCCATGTGTTAAATATTTAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAg                                    <  2:2617630/100‑1 (MQ=255)

                                                                  aTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                   <  1:2431949/100‑1 (MQ=255)

                                                                    ctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGccc                                 >  1:2196543/1‑100 (MQ=255)

                                                                     tctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCa                                <  1:2990748/100‑1 (MQ=255)

                                                                      ctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCAt                               >  1:1735108/1‑100 (MQ=255)

                                                                         cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTAc                            <  2:3440349/100‑1 (MQ=255)

                                                                          aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACt                           <  1:1221103/100‑1 (MQ=255)

                                                                            gtgtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTg                         <  2:1670221/100‑1 (MQ=255)

                                                                              gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAg                       >  1:311512/1‑100 (MQ=255)

                                                                                tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                     >  1:2940662/1‑100 (MQ=255)

                                                                                tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                     <  2:3041717/100‑1 (MQ=255)

                                                                                tAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGcc                     <  2:2479073/100‑1 (MQ=255)

                                                                                 aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGccc                    >  1:2080500/1‑100 (MQ=255)

                                                                                 aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGccc                    <  2:811463/100‑1 (MQ=255)

                                                                                  aataaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                   >  1:386060/9‑100 (MQ=37)

                                                                                  aataaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                   >  1:3635429/9‑100 (MQ=37)

                                                                                    taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                 >  1:564345/7‑100 (MQ=38)

                                                                                    taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                 >  1:1661659/7‑100 (MQ=38)

                                                                                     atTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCTATTACTTGAGCCCatat                >  2:547935/1‑100 (MQ=255)

                                                                                     atTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                >  1:879426/1‑100 (MQ=255)

                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                >  2:245961/1‑99 (MQ=255)

                                                                                      ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg               <  2:490283/100‑1 (MQ=255)

                                                                                      aaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg               >  1:114563/5‑100 (MQ=38)

                                                                                        gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCAATTACTTTAGCCCATATg               >  2:1559374/1‑98 (MQ=255)

                                                                                         caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc            >  2:453590/2‑100 (MQ=255)

                                                                                         caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc            >  2:2507559/2‑100 (MQ=255)

                                                                                         caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc            >  1:417177/2‑100 (MQ=255)

                                                                                         aaaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc            <  1:981738/100‑1 (MQ=255)

                                                                                          aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa           >  2:3277866/1‑100 (MQ=255)

                                                                                          aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa           >  1:2339959/1‑100 (MQ=255)

                                                                                          aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa           <  1:1331167/100‑1 (MQ=255)

                                                                                          aaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCa           <  1:3032095/100‑1 (MQ=255)

                                                                                              ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAtt       >  1:1809444/1‑100 (MQ=255)

                                                                                               gCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttt      >  1:1382286/1‑100 (MQ=255)

                                                                                                cGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATGTGAAGCCCATTACTTGAGCCCATATGGGCATAtttt     >  1:3469252/1‑100 (MQ=255)

                                                                                                cGTATAACACGCCCGTTGTTTTAATTATGTGGATATTATTAATAGCATATCGCGCATATTTATATGAAGCCAATTACTTGAGCCCATATGGGGATAtttt     >  2:1951411/1‑100 (MQ=255)

                                                                                                 gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt    <  2:585558/100‑1 (MQ=255)

                                                                                                   ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat  >  2:2675605/1‑100 (MQ=255)

                                                                                                   ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat  >  2:1886470/1‑100 (MQ=255)

TCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTAT  >  NC_000913/1211511‑1211709

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