BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,211,614	C→T	18.3%	intergenic (+37/+66)	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,614	0	C	T	18.3%	457.9 / 80.2	195	intergenic (+37/+66)	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 C (84/75); new base T (20/16); total (104/91)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 207 total aligned reads displayed.

ATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATG  >  NC_000913/1211515‑1211713

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

 tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                    <  1:3139838/100‑1 (MQ=255)

 tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTCGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAaca                                                                                                    <  1:2319095/100‑1 (MQ=255)

  ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                                                                                   >  2:3690007/1‑100 (MQ=255)

     aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                                <  2:2654631/100‑1 (MQ=255)

     aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                                <  1:1861456/100‑1 (MQ=255)

       gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAATCGATCATCGAAAACATGGAATCACCCCCCCTGTTAAATATTGAAACGGGCGTATAACACGcc                                                                                                <  1:1371064/98‑1 (MQ=21)

       gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCg                                                                                              <  2:466101/100‑1 (MQ=25)

       gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGccc                                                                                               <  1:3147599/99‑1 (MQ=255)

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

        cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                                                                             <  1:2063856/100‑1 (MQ=255)

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

           tAAGCTGCTGAAATGTTCAGAGTTTGTTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGACACGGGCGTGTTATACGCCCgtttc                                                                                          <  1:205134/100‑3 (MQ=25)

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

             aGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                                                                                        >  2:2218115/1‑100 (MQ=255)

               ctgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTttta                                                                                      <  2:2290918/100‑1 (MQ=255)

               ctgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTttta                                                                                      <  1:1754017/100‑1 (MQ=255)

                tgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttat                                                                                     <  2:1174172/100‑1 (MQ=255)

                  ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                   <  2:840675/100‑1 (MQ=255)

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

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

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

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

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

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

                        tGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTgg                                                                             >  2:1855857/1‑100 (MQ=255)

                         gTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCACGTTGTTTTATTTATATGGa                                                                            >  1:4248547/1‑100 (MQ=255)

                          ttCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTgggt                                                                           >  2:2377271/1‑98 (MQ=255)

                           tCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGata                                                                          <  2:2841481/100‑1 (MQ=255)

                             agagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtatt                                                                        >  2:1498717/1‑100 (MQ=255)

                              gagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACCCGCCCGTTGTTTTATTTATGTGGAtatta                                                                       >  1:1222097/1‑100 (MQ=255)

                               agTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtattat                                                                      >  2:60886/1‑100 (MQ=255)

                                gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                     <  2:2168585/100‑1 (MQ=255)

                                gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                     >  2:556217/1‑100 (MQ=255)

                                    ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                                 <  1:4139830/100‑1 (MQ=255)

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

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

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

                                      tGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGc                                                               <  1:544929/100‑1 (MQ=255)

                                       gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGc                                                               <  2:3510631/99‑1 (MQ=255)

                                       gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCa                                                              <  2:3431810/100‑1 (MQ=255)

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

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

                                       gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTGAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCa                                                              <  2:4392122/100‑1 (MQ=255)

                                         aGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCata                                                            <  2:3470392/100‑1 (MQ=255)

                                          gCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCatat                                                           >  2:1918872/1‑100 (MQ=255)

                                              tcatcGAAAACATGTAATCTCTCCATGTGTTAAATACTGAATCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAg                                                       <  1:2426197/100‑1 (MQ=255)

                                               catcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                      <  1:3335926/100‑1 (MQ=255)

                                                atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                      <  1:1312459/99‑1 (MQ=255)

                                                 tcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCat                                                    >  2:1068459/1‑100 (MQ=255)

                                                  cGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                                   <  2:591985/100‑1 (MQ=255)

                                                   gAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGCGCatat                                                  >  1:1125220/1‑100 (MQ=255)

                                                   gAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCatat                                                  >  2:510098/1‑100 (MQ=255)

                                                   gAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCatat                                                  >  1:4311275/1‑100 (MQ=255)

                                                     aaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATAttt                                                <  1:288355/100‑1 (MQ=255)

                                                      aaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTta                                               >  2:434392/1‑100 (MQ=255)

                                                        cATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                             <  1:1051106/100‑1 (MQ=255)

                                                        cATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                             <  1:3683927/100‑1 (MQ=255)

                                                         atttAATCTCTCCATGTGTTAAATATTTAAACTGGCTTATAACACTCCCGTTGTTTTATTTATTTTGATATTATTAATAGCATATCGAGCATATTTatat                                            <  1:2532533/97‑1 (MQ=255)

                                                          tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                           >  2:3056233/1‑100 (MQ=255)

                                                          tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                           >  2:1930221/1‑100 (MQ=255)

                                                           gTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGa                                          >  1:3003656/1‑100 (MQ=255)

                                                           gTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGa                                          >  2:3793944/1‑100 (MQ=255)

                                                            tAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGaa                                         >  2:3814413/1‑100 (MQ=255)

                                                               tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTGTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGc                                       >  2:3632691/1‑99 (MQ=255)

                                                                   tcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCAtt                                  >  2:1376257/1‑100 (MQ=255)

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

                                                                          tttAAATAGTGAAACGGGCGTGTAACACGCCCGTTGGTTTATTAATGTGGATATTATTAATAGCATATCTAGCATTTTTATATGATGCCCAGTAGTTGa                            >  1:2622811/2‑99 (MQ=255)

                                                                          gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAg                           <  1:1792275/100‑1 (MQ=255)

                                                                           ttAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGc                          <  2:3272879/100‑1 (MQ=255)

                                                                            tAAATANTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTTTTAATAGCATATCGAGCGTATTTATATGAAGCCNATTACTTGAGc                          >  1:2565080/1‑99 (MQ=255)

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

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

                                                                              aaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                       >  2:1980099/1‑100 (MQ=255)

                                                                               atatTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCat                      >  2:702480/1‑100 (MQ=255)

                                                                                taaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                     >  2:3403144/7‑100 (MQ=38)

                                                                                 aaaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGAATATTTATATGAAACCCATTACTTTAGCCtata                     >  1:4384949/6‑99 (MQ=25)

                                                                                  ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg                   >  2:204242/1‑100 (MQ=255)

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

                                                                                  aaacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATg                   >  2:1352574/5‑100 (MQ=38)

                                                                                    gAAACGTGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                 <  1:2811961/100‑1 (MQ=255)

                                                                                    acaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATggg                 >  2:2755130/3‑100 (MQ=255)

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

                                                                                     aaaCGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc                >  2:918176/1‑100 (MQ=255)

                                                                                       aCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCat              >  2:1279236/1‑100 (MQ=255)

                                                                                       aCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCat              >  2:2470042/1‑100 (MQ=255)

                                                                                        cGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCata             >  1:3020349/1‑100 (MQ=255)

                                                                                        cGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCata             <  1:3811226/100‑1 (MQ=255)

                                                                                          ggCGTGTAATACGCCCGTTTTGTGATTTATGTGGATATTATTAATAGGATAACGAGCATATGTATATGAAGCCCAGTATTTGAGCCCATATGGGCata             <  1:4107932/98‑1 (MQ=255)

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

                                                                                            cGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAtttt         >  1:2760090/1‑100 (MQ=255)

                                                                                             gTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt        >  1:4170598/1‑100 (MQ=255)

                                                                                              tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta       <  2:3992650/100‑1 (MQ=255)

                                                                                               ataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtat      <  2:3003211/100‑1 (MQ=255)

                                                                                                taACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtata     >  2:4294874/1‑100 (MQ=255)

                                                                                                 aaCACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttagatc    >  1:3475244/1‑95 (MQ=255)

                                                                                                  acacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt   <  2:2037077/100‑1 (MQ=255)

                                                                                                  acacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt   >  2:1381182/1‑100 (MQ=255)

                                                                                                  acacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt   >  2:2905796/1‑100 (MQ=255)

                                                                                                   cacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg  >  2:1765725/1‑100 (MQ=255)

ATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATG  >  NC_000913/1211515‑1211713

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