BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,361,329	G→A	G52S (GGC→AGC)	gltF →	periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,361,329	0	G	A	99.0%	350.6 / ‑4.7	104	G52S (GGC→AGC)	gltF	periplasmic protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (52/51); total (52/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 105 total aligned reads displayed.

GCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCGGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAGCACATTA  >  NC_000913/3361233‑3361424

gcaTTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCGGCa                                                                                              <  1:1436584/100‑1 (MQ=255)

     cAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCAttgtt                                                                                         <  1:431426/100‑1 (MQ=255)

     cAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCAttgtt                                                                                         <  2:343538/100‑1 (MQ=255)

      agggCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGATTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCAttgttg                                                                                        <  1:903192/97‑1 (MQ=255)

      aGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCAttgttg                                                                                        <  1:1244001/100‑1 (MQ=255)

      aGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCAttgttg                                                                                        <  2:371895/100‑1 (MQ=255)

      aGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCAttgttg                                                                                        <  1:2424211/100‑1 (MQ=255)

       gTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGa                                                                                       >  2:599611/1‑100 (MQ=255)

       gTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGa                                                                                       >  1:483223/1‑100 (MQ=255)

          ccATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTa                                                                                    >  1:116770/1‑100 (MQ=255)

          ccATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTa                                                                                    <  1:150268/100‑1 (MQ=255)

           cATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTAt                                                                                   >  2:2192513/1‑100 (MQ=255)

            aTGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATg                                                                                  <  1:1804865/100‑1 (MQ=255)

            aTGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATg                                                                                  >  1:1623185/1‑100 (MQ=255)

              ggCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATggg                                                                                <  2:1076314/100‑1 (MQ=255)

               gCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGa                                                                               >  1:2269309/1‑100 (MQ=255)

               gCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGa                                                                               <  1:499586/100‑1 (MQ=255)

               gCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGa                                                                               >  1:430271/1‑100 (MQ=255)

                cAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGaa                                                                              <  1:2618132/100‑1 (MQ=255)

                cAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGaa                                                                              <  2:1746241/100‑1 (MQ=255)

                  aCAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAAc                                                                            >  2:736468/1‑100 (MQ=255)

                    aGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAAcat                                                                          <  2:233651/100‑1 (MQ=255)

                         cAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATaa                                                                     >  1:1474854/1‑100 (MQ=255)

                            cTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAg                                                                  >  1:2006057/1‑100 (MQ=255)

                            cTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAg                                                                  >  1:952364/1‑100 (MQ=255)

                            cTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAg                                                                  >  2:267034/1‑100 (MQ=255)

                            cTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAg                                                                  >  1:894154/1‑100 (MQ=255)

                            cTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAg                                                                  >  1:1573155/1‑100 (MQ=255)

                            cTGATACACAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAg                                                                  <  1:1639738/100‑1 (MQ=255)

                             tGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATAATAACAGt                                                                 >  2:1411419/1‑100 (MQ=255)

                              gATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTg                                                                >  1:2616616/1‑100 (MQ=255)

                                  cAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCt                                                            <  1:1369332/100‑1 (MQ=255)

                                  cAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGAGCt                                                            >  2:1206223/1‑100 (MQ=255)

                                   aGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTg                                                           >  2:2008438/1‑100 (MQ=255)

                                   aGAATCAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTg                                                           <  1:207764/100‑1 (MQ=255)

                                    gaaacaaCGATCAACAGCAAGTATACCCCTGAGAAATGTAAACGAGTAGATATGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGa                                                          >  2:1930647/6‑100 (MQ=255)

                                    gAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGa                                                          >  1:470695/1‑100 (MQ=255)

                                     aaTTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGaa                                                         >  2:438154/1‑100 (MQ=255)

                                      aTTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAAt                                                        <  1:1938187/100‑1 (MQ=255)

                                      aTTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAAt                                                        <  2:1231740/100‑1 (MQ=255)

                                       ttAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATc                                                       <  1:708076/100‑1 (MQ=255)

                                       ttAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATc                                                       <  1:1590306/100‑1 (MQ=255)

                                       ttAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATc                                                       <  1:2160811/100‑1 (MQ=255)

                                        tAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATcc                                                      <  2:92858/100‑1 (MQ=255)

                                         aaCTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCg                                                     <  2:2104633/100‑1 (MQ=255)

                                          aCTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGa                                                    >  1:2507360/1‑100 (MQ=255)

                                            tATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGAcc                                                  <  1:1926810/100‑1 (MQ=255)

                                            tATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCAGAACAGTGCGCTGAATCCGAcc                                                  >  2:1936148/1‑100 (MQ=255)

                                            tATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGCTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGAcc                                                  <  2:2503228/100‑1 (MQ=255)

                                             atcatcGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCg                                                 <  2:114986/100‑1 (MQ=255)

                                              tcatcGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGGTGACTATGGGAAACATCATAGCTGTGCGCTGAATCCGACCgg                                                <  1:2251992/100‑1 (MQ=255)

                                                 tcGGCGGGTATCCCCCGGCTGCATGGACCCCGGTTGTTACGGGCAGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGCATCCGACCGGTaa                                             <  2:208306/100‑1 (MQ=255)

                                                 tcGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTaa                                             >  1:2389974/1‑100 (MQ=255)

                                                   ggCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTaaaa                                           >  1:256878/1‑100 (MQ=255)

                                                    gCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAg                                          <  1:1971901/100‑1 (MQ=255)

                                                    gCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAg                                          <  1:1417884/100‑1 (MQ=255)

                                                      gAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGta                                        <  2:881613/100‑1 (MQ=255)

                                                      gAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGta                                        <  2:392060/100‑1 (MQ=255)

                                                       aGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGtaa                                       >  1:357099/1‑100 (MQ=255)

                                                       aGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGtaa                                       >  1:1922164/1‑100 (MQ=255)

                                                         tataCCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGtaata                                     <  1:1453882/100‑1 (MQ=255)

                                                           taCCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATaaa                                   >  1:1582788/1‑100 (MQ=255)

                                                            aCCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAAc                                  >  1:1125307/1‑100 (MQ=255)

                                                             ccccTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACt                                 >  1:2289724/1‑100 (MQ=255)

                                                             ccccTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACt                                 >  1:370743/1‑100 (MQ=255)

                                                             ccccTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACt                                 >  2:1579369/1‑100 (MQ=255)

                                                             ccccTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACt                                 >  2:1578865/1‑100 (MQ=255)

                                                             ccccTGGTGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACCATGGGAGACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACt                                 >  2:2423545/1‑100 (MQ=255)

                                                                    tGCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCaa                          <  2:1719657/100‑1 (MQ=255)

                                                                     gCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAAc                         >  1:2043574/1‑100 (MQ=255)

                                                                     gCATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAAc                         <  2:1097388/100‑1 (MQ=255)

                                                                      cATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACt                        >  2:1055410/1‑100 (MQ=255)

                                                                      cATGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCAGACCGGTAAAAGTAATAAACTGGTTCAACt                        >  1:1631491/1‑100 (MQ=255)

                                                                       aTGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTg                       >  2:1913168/1‑100 (MQ=255)

                                                                        tGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTgg                      >  1:2628598/1‑100 (MQ=255)

                                                                        tGTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTgg                      <  1:1811149/100‑1 (MQ=255)

                                                                         gTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTggg                     >  2:858652/1‑100 (MQ=255)

                                                                         gTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTggg                     <  1:1236804/100‑1 (MQ=255)

                                                                         gTACCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTggg                     <  2:1334664/100‑1 (MQ=255)

                                                                           aCCCCAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGcc                   <  2:2466266/100‑1 (MQ=255)

                                                                            ccccAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGACTCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCg                  >  2:307446/1‑100 (MQ=255)

                                                                            ccccAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCg                  <  2:2245668/100‑1 (MQ=255)

                                                                             cccAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGc                 <  1:1431346/100‑1 (MQ=255)

                                                                             cccAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGc                 >  1:1186357/1‑100 (MQ=255)

                                                                             cccAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGc                 <  2:2293950/100‑1 (MQ=255)

                                                                              ccAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCa                >  1:1081902/1‑100 (MQ=255)

                                                                              ccAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCa                <  1:1842454/100‑1 (MQ=255)

                                                                               cAGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCaa               <  1:2461443/100‑1 (MQ=255)

                                                                                aGTTGTTACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGCGCCGCaaa              <  2:2485689/100‑1 (MQ=255)

                                                                                   tgttACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAAc           <  1:1457497/100‑1 (MQ=255)

                                                                                   tgttACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAAc           >  2:1586270/1‑100 (MQ=255)

                                                                                   tgttACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAAc           >  1:2194050/1‑100 (MQ=255)

                                                                                   tgttACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAAc           >  1:546253/1‑100 (MQ=255)

                                                                                     ttACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAg         >  1:2599718/1‑100 (MQ=255)

                                                                                      tACGGGCGGCAGCATTGTTGACTATGGGAAAGATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAGc        >  1:2115780/1‑100 (MQ=255)

                                                                                      tACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAGc        >  2:1647948/1‑100 (MQ=255)

                                                                                      tACGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAGc        >  1:515559/1‑100 (MQ=255)

                                                                                        cGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGTGCCGCAAAAACAGcac      <  1:666987/100‑1 (MQ=255)

                                                                                        cGGGCGGCAGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAGcac      <  1:667137/100‑1 (MQ=255)

                                                                                            cggcagcATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAGCACATTa  <  2:1497380/100‑1 (MQ=255)

GCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTATACCCCTGGTGCATGTACCCCAGTTGTTACGGGCGGCGGCATTGTTGACTATGGGAAACATCATAACAGTGCGCTGAATCCGACCGGTAAAAGTAATAAACTGGTTCAACTGGGCCGCAAAAACAGCACATTA  >  NC_000913/3361233‑3361424

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