BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	146 [0]	[0] 142	insB1–insA	insB1, insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	122 (1.170)	72/182	0.2	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
*	?	NC_000913	1979279 =	0 (0.000)	122 (1.170)	72/182	0.2	100%	intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

Only 100 of 146 total aligned reads displayed.

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601

cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  1:2101958‑M1/100‑1 (MQ=255)

cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  1:656338‑M1/100‑1 (MQ=255)

cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  2:1482146‑M1/100‑1 (MQ=255)

 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                    >  1:2414036‑M1/1‑99 (MQ=255)

  gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                   >  1:357701‑M1/1‑98 (MQ=255)

   tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                  <  2:1288730‑M1/100‑4 (MQ=255)

      tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTCATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                               <  1:1605815‑M1/100‑7 (MQ=255)

         tCACTACACGCACATACAACGGAGGGTGGCTGAGATTTTCAATAATGCGTAATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                             <  1:1729525‑M1/99‑9 (MQ=255)

         tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                            >  1:930439‑M1/1‑91 (MQ=255)

          cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                           >  1:1875993‑M1/1‑90 (MQ=255)

              acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                       >  2:2357908‑M1/1‑86 (MQ=255)

                 cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    >  2:458530‑M1/1‑83 (MQ=255)

                 cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    <  2:2252219‑M1/100‑18 (MQ=255)

                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                               >  2:2844354‑M1/1‑78 (MQ=255)

                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                 <  2:2446405‑M1/98‑21 (MQ=255)

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  1:2309864‑M1/1‑77 (MQ=255)

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              <  2:1028107‑M1/100‑24 (MQ=255)

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  1:1530556‑M1/1‑77 (MQ=255)

                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:1138437‑M1/100‑25 (MQ=255)

                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  1:857973‑M1/100‑25 (MQ=255)

                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  2:328733‑M1/1‑75 (MQ=255)

                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                           <  1:1697087‑M1/100‑27 (MQ=255)

                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                           <  2:2489927‑M1/100‑27 (MQ=255)

                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     >  1:2247446‑M1/1‑68 (MQ=255)

                                  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                    >  1:881872‑M1/1‑66 (MQ=255)

                                   gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                  >  1:2871643‑M1/1‑65 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 >  2:1745467‑M1/1‑64 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 >  1:2044618‑M1/1‑64 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 >  2:1090586‑M1/1‑64 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 <  2:741597‑M1/100‑37 (MQ=255)

                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                               >  1:1315371‑M1/1‑62 (MQ=255)

                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                >  1:1077257‑M1/1‑62 (MQ=255)

                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             <  1:1103862‑M1/100‑41 (MQ=255)

                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  2:1936453‑M1/1‑60 (MQ=255)

                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  1:331601‑M1/1‑60 (MQ=255)

                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  1:2581756‑M1/1‑60 (MQ=255)

                                          gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCTATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                            <  1:1507388‑M1/99‑42 (MQ=255)

                                                naataatgCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                     >  2:695618‑M1/2‑52 (MQ=255)

                                                cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                     >  2:2518130‑M1/1‑52 (MQ=255)

                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                    >  1:2341300‑M1/1‑51 (MQ=255)

                                                  ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                   <  1:1996391‑M1/100‑51 (MQ=255)

                                                   taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                  >  2:2034392‑M1/1‑49 (MQ=255)

                                                   taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                  >  2:162849‑M1/1‑49 (MQ=255)

                                                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                 <  2:2878914‑M1/100‑53 (MQ=255)

                                                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                 <  2:1370946‑M1/100‑53 (MQ=255)

                                                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                <  1:2015552‑M1/100‑54 (MQ=255)

                                                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                            >  1:1366737‑M1/1‑43 (MQ=255)

                                                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                            >  2:2440755‑M1/1‑43 (MQ=255)

                                                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                            >  2:1128224‑M1/1‑43 (MQ=255)

                                                            aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         >  2:839885‑M1/1‑40 (MQ=255)

                                                            aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         <  2:248060‑M1/100‑61 (MQ=255)

                                                            aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         <  2:1055468‑M1/100‑61 (MQ=255)

                                                            aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         <  1:419128‑M1/100‑61 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgcta                                        >  1:2427908‑M1/1‑39 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                        <  1:328731‑M1/100‑62 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                        >  2:192988‑M1/1‑39 (MQ=255)

                                                               nagATCACACAAAACACTCAATTACTTAACATAAATGtataagtgatacttgtgttttggctgtatttcaatctgttaaaaagtttgtcagtacgcta                                        >  1:1134007‑M1/2‑37 (MQ=255)

                                                               cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                      >  2:1410765‑M1/1‑37 (MQ=255)

                                                               cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                       >  1:1516661‑M1/1‑37 (MQ=255)

                                                                aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                     <  1:2506405‑M1/100‑65 (MQ=255)

                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                    >  1:1448181‑M1/1‑35 (MQ=255)

                                                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                     <  1:1219698‑M1/97‑65 (MQ=255)

                                                                      cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                               >  1:45461‑M1/1‑30 (MQ=255)

                                                                      cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                               >  2:1913769‑M1/1‑30 (MQ=255)

                                                                        caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                             <  2:2564051‑M1/100‑73 (MQ=255)

                                                                        caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                             >  2:2102842‑M1/1‑28 (MQ=255)

                                                                           aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                          >  2:185576‑M1/1‑25 (MQ=255)

                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                      >  1:1192825‑M1/1‑21 (MQ=255)

                                                                                tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                     <  2:376394‑M1/100‑81 (MQ=255)

                                                                                tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                     <  2:1149470‑M1/100‑81 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    >  2:314253‑M1/1‑19 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    >  2:779633‑M1/1‑19 (MQ=255)

                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                   <  2:2078811‑M1/100‑83 (MQ=255)

                                                                                   aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                  >  2:1971902‑M1/1‑17 (MQ=255)

                                                                                     tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                <  2:1840550‑M1/100‑86 (MQ=255)

                                                                                     tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                <  2:1775245‑M1/100‑86 (MQ=255)

                                                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa               >  2:2811500‑M1/1‑14 (MQ=255)

                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat              >  2:2298585‑M1/1‑13 (MQ=255)

                                                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc             <  2:2024657‑M1/100‑89 (MQ=255)

                                                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaaaatgaaacaggaataatcg            >  1:650258‑M1/1‑11 (MQ=255)

                                                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg            >  1:1906003‑M1/1‑11 (MQ=255)

                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa          >  2:2676930‑M1/1‑9 (MQ=255)

                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa          >  2:921309‑M1/1‑9 (MQ=255)

                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa          >  1:789962‑M1/1‑9 (MQ=255)

                                                                                            cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa         >  1:1945796‑M1/1‑8 (MQ=255)

                                                                                            cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa         >  1:2726254‑M1/1‑8 (MQ=255)

                                                                                            cATAAATGtataagtcatacttttgatttgggtgtatttcaatctgctaaaaagttgttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa         >  2:2703107‑M1/1‑8 (MQ=255)

                                                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat        >  1:2187306‑M1/1‑7 (MQ=255)

                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       >  1:948349‑M1/1‑6 (MQ=255)

                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcaacgctagcaagcaaaaatgaaacaggaataatcgaaatg       >  1:2024281‑M1/1‑6 (MQ=255)

                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     >  2:1518960‑M1/1‑4 (MQ=255)

                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     >  1:1299264‑M1/1‑4 (MQ=255)

                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     <  2:45461‑M1/100‑97 (MQ=255)

                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     >  2:662223‑M1/1‑4 (MQ=255)

                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga    <  1:1418464‑M1/100‑98 (MQ=255)

                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga    >  1:31869‑M1/1‑3 (MQ=255)

                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga    <  1:2805761‑M1/100‑98 (MQ=255)

                                                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat   <  1:2844328‑M1/100‑99 (MQ=255)

                                                                                                   gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  <  1:705523‑M1/100‑100 (MQ=255)

                                                                                                   gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  <  2:397277‑M1/100‑100 (MQ=255)

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601

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