BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	691,390	G→C	16.7%	A132G (GCC→GGC)	ybeX ←	putative ion transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	691,390	0	G	C	16.7%	322.8 / 56.0	125	A132G (GCC→GGC)	ybeX	putative ion transport
Reads supporting (aligned to +/- strand): ref base G (52/52); new base C (10/11); total (62/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 128 total aligned reads displayed.

GTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGG  >  NC_000913/691293‑691488

gtgGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCt                                                                                                  >  1:1043849/1‑100 (MQ=255)

  ggTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTc                                                                                                >  2:1009782/1‑100 (MQ=255)

   gTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                               <  2:696943/100‑1 (MQ=255)

     aaCGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGc                                                                                             >  1:576914/1‑100 (MQ=255)

       cGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                                           >  1:2048743/1‑100 (MQ=255)

        gCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                          <  2:1799522/100‑2 (MQ=255)

         cTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCg                                                                                         >  2:959611/1‑100 (MQ=255)

         cTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCg                                                                                         <  1:522525/100‑1 (MQ=255)

         cTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgg                                                                                         <  2:1069961/100‑3 (MQ=255)

         cTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgg                                                                                         <  1:386937/100‑3 (MQ=255)

          tGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgga                                                                                        <  1:1016984/100‑4 (MQ=255)

          tGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgg                                                                                         <  2:2559138/99‑3 (MQ=255)

           gagaGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCCGAAGCCTTCAGCAtggac                                                                                       <  1:1560084/100‑5 (MQ=255)

            agagCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtggaca                                                                                      <  2:1466303/100‑6 (MQ=255)

            agagCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtggaca                                                                                      <  2:1466008/100‑6 (MQ=255)

             gagCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgc                                                                                     >  2:1862686/1‑100 (MQ=255)

              agCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcg                                                                                    <  2:885088/100‑1 (MQ=255)

              agCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcg                                                                                    <  2:2517677/100‑1 (MQ=255)

               gCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcgc                                                                                   >  2:1221408/1‑100 (MQ=255)

                cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCa                                                                                  >  1:18797/1‑100 (MQ=255)

                       ctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACgg                                                                           >  2:262727/1‑100 (MQ=255)

                        tttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgc                                                                          >  1:773208/1‑100 (MQ=255)

                        tttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgc                                                                          >  2:392349/1‑100 (MQ=255)

                         ttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgca                                                                         >  1:748101/1‑100 (MQ=255)

                         ttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgca                                                                         >  1:2021472/1‑100 (MQ=255)

                         ttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgca                                                                         >  1:1524443/1‑100 (MQ=255)

                         ttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgca                                                                         >  2:2360738/1‑100 (MQ=255)

                          tCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgcag                                                                        >  1:1736758/1‑100 (MQ=255)

                          tCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgcag                                                                        <  2:1068466/100‑1 (MQ=255)

                           cAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgcagc                                                                       <  1:534507/100‑1 (MQ=255)

                             gCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCATCATCGCTGCGCATAAACGGCAGCaa                                                                     >  2:1517504/1‑100 (MQ=255)

                              cATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgcagcacg                                                                    >  1:1558711/1‑98 (MQ=255)

                              cATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAg                                                                    <  2:1443646/100‑1 (MQ=255)

                                tCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTc                                                                  <  1:2444873/100‑1 (MQ=255)

                                tCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTc                                                                  >  1:604727/1‑100 (MQ=255)

                                     tCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                                             <  2:2158860/100‑1 (MQ=255)

                                     tCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                                             >  2:2535345/1‑100 (MQ=255)

                                      cTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                                             <  1:1221404/99‑1 (MQ=255)

                                      cTACGCGCTTACTTTCAGAAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGc                                                            <  1:357218/100‑1 (MQ=255)

                                       tACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                                             <  1:835987/98‑1 (MQ=255)

                                       tACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGcc                                                           >  1:840656/1‑100 (MQ=255)

                                        aCGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCa                                                          <  1:34126/100‑1 (MQ=255)

                                           cgcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCa                                                       >  1:1096964/1‑100 (MQ=255)

                                            gcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAg                                                      >  2:1252916/1‑100 (MQ=255)

                                             cTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGa                                                     >  2:2483264/1‑100 (MQ=255)

                                                aCTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAAta                                                  <  2:870302/100‑2 (MQ=255)

                                                 cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATcc                                                 >  1:697869/1‑100 (MQ=255)

                                                 cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATcc                                                 >  1:1481162/1‑100 (MQ=255)

                                                 cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATcc                                                 >  1:136751/1‑100 (MQ=255)

                                                   ttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCt                                               >  1:75802/1‑100 (MQ=255)

                                                     cAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTc                                             <  1:262724/100‑1 (MQ=255)

                                                     cAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTc                                             <  1:936329/100‑1 (MQ=255)

                                                      aGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCa                                            <  1:1219323/100‑1 (MQ=255)

                                                      aGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCa                                            >  1:1648268/1‑100 (MQ=255)

                                                      aGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCa                                            <  2:2214740/100‑1 (MQ=255)

                                                       ggAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCaa                                           >  1:411470/1‑100 (MQ=255)

                                                         aaCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtg                                         <  2:2275560/100‑1 (MQ=255)

                                                         aaCGACAACCGCCTGACGTAACACTTTGTACATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGTAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtg                                         <  1:1553951/100‑1 (MQ=255)

                                                          aCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtgt                                        >  1:1935104/1‑100 (MQ=255)

                                                           cGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtgtg                                       <  1:1303345/100‑1 (MQ=255)

                                                            gACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGa                                      <  2:356839/100‑1 (MQ=255)

                                                             aCAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGTAGCAAGTCTTTCGCCATCAGAATCCTTTCAATGTGAt                                     <  1:1933649/100‑1 (MQ=255)

                                                             aCAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGAt                                     >  2:782951/1‑100 (MQ=255)

                                                                  cGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                <  2:1935114/100‑1 (MQ=255)

                                                                  cGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                <  2:411472/100‑1 (MQ=255)

                                                                  cGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                <  2:1395714/100‑1 (MQ=255)

                                                                  cGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                <  2:1524452/100‑1 (MQ=255)

                                                                    ccTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                <  1:1983697/98‑1 (MQ=255)

                                                                    ccTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTc                              <  1:1599049/100‑1 (MQ=255)

                                                                    ccTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTc                              <  2:751868/100‑1 (MQ=255)

                                                                       gACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTc                           <  1:1178459/100‑1 (MQ=255)

                                                                       gACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTc                           >  1:588348/1‑100 (MQ=255)

                                                                        aCGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCg                          <  1:1431763/100‑1 (MQ=255)

                                                                         cGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                         >  1:2585145/1‑100 (MQ=255)

                                                                         cGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                         <  1:2297328/100‑1 (MQ=255)

                                                                              cacTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAACCGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATc                    <  1:248248/100‑1 (MQ=255)

                                                                                cTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCAc                  >  2:1279029/1‑100 (MQ=255)

                                                                                cTTTGTCCATGCTGAAGGCTTCAGCATCACTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCAc                  >  2:743461/1‑100 (MQ=255)

                                                                                   tGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg               <  2:588350/100‑1 (MQ=255)

                                                                                   tGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg               <  1:2570632/100‑1 (MQ=255)

                                                                                   tGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg               <  2:1321726/100‑1 (MQ=255)

                                                                                   tGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg               >  2:2081301/1‑100 (MQ=255)

                                                                                   cagcgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg               >  1:2578469/6‑100 (MQ=255)

                                                                                   cagcgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg               >  2:1737644/6‑100 (MQ=255)

                                                                                    agcgaTGCTAAAGCCTGCAGCATCGCCGCGCATAAACGGCTGCACGTCTTTCGCCATCAGAATCCCTTCAATGTGGTGTTTTCCTTCGCTAATCACCgg               >  1:403253/5‑99 (MQ=255)

                                                                                     gcgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCggg              >  2:1395904/4‑99 (MQ=255)

                                                                                     gcgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGa             >  2:337935/4‑100 (MQ=255)

                                                                                      ccATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGaa            >  1:862627/1‑100 (MQ=255)

                                                                                       gaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGaaa           >  2:1660553/2‑100 (MQ=255)

                                                                                       gaTGCTGAAGCCGTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGaaa           >  2:1704959/2‑100 (MQ=255)

                                                                                       cATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGaaa           >  2:407862/1‑100 (MQ=255)

                                                                                         tGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACg         <  2:2126852/100‑1 (MQ=255)

                                                                                          gCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGt        >  1:229901/1‑100 (MQ=255)

                                                                                           cTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg       <  1:1089694/100‑1 (MQ=255)

                                                                                           cTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg       <  2:774517/100‑1 (MQ=255)

                                                                                           cTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg       >  2:2220472/1‑100 (MQ=255)

                                                                                           cTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg       >  1:205522/1‑100 (MQ=255)

                                                                                                ggCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTgg  >  2:598201/1‑100 (MQ=255)

                                                                                                ggCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTgg  >  2:2575718/1‑100 (MQ=255)

                                                                                                ggCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTgg  >  2:1353145/1‑100 (MQ=255)

GTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGG  >  NC_000913/691293‑691488

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