BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,989,886	T→A	34.7%	D5V (GAT→GTT)	yqeL ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,989,886	0	T	A	34.7%	162.8 / 157.4	150	D5V (GAT→GTT)	yqeL	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (63/35); new base A (25/27); total (88/62)
Fisher's exact test for biased strand distribution p-value = 5.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

Only 100 of 153 total aligned reads displayed.

ACTAAAATAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATCTTAG  >  NC_000913/2989787‑2989984

aCTAAAATAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAt                                                                                                    <  2:1766835/100‑1 (MQ=255)

aCTAAAATAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAt                                                                                                    <  2:80651/100‑1 (MQ=255)

 cTAAAATAGCTAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTAATc                                                                                                   >  2:2137011/1‑100 (MQ=255)

    aaaTAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCaan                                                                                                <  2:394054/100‑2 (MQ=255)

    aaaTAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAAc                                                                                                >  2:2368962/1‑100 (MQ=255)

    aaaTAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAAc                                                                                                <  2:19806/100‑1 (MQ=255)

     aaTAGATAACACAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACa                                                                                               >  2:2034351/1‑100 (MQ=255)

      aTAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACAt                                                                                              >  1:2255544/1‑100 (MQ=255)

       tAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATc                                                                                             >  1:49488/1‑100 (MQ=255)

         gATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCtt                                                                                           >  1:2169597/1‑100 (MQ=255)

          aTAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCttt                                                                                          >  1:1774798/1‑100 (MQ=255)

          aTAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCttt                                                                                          <  1:1281331/100‑1 (MQ=255)

           tAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTc                                                                                         >  2:2049830/1‑100 (MQ=255)

              cAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCAtt                                                                                      >  1:2048500/1‑100 (MQ=255)

                aaaaGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCAtttt                                                                                    >  2:726860/1‑100 (MQ=255)

                aaaaGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCAtttt                                                                                    <  1:257513/100‑1 (MQ=255)

                 aaaGATAGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTa                                                                                   <  1:794339/100‑1 (MQ=255)

                  aaGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAg                                                                                  <  2:836781/100‑1 (MQ=255)

                   aGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAg                                                                                  <  1:246392/99‑1 (MQ=255)

                    gATGGAACCTCGTTAAGCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGAc                                                                                <  2:1398658/100‑1 (MQ=255)

                    atctgAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGa                                                                                 <  2:1189156/95‑1 (MQ=255)

                       ggAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGACAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACAtc                                                                             >  2:2341351/1‑100 (MQ=255)

                        gAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACAtct                                                                            >  2:1314149/1‑100 (MQ=255)

                         aaCCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACAtctc                                                                           <  2:277747/100‑1 (MQ=255)

                          aCCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTcc                                                                          >  2:152632/1‑100 (MQ=255)

                           ccTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCt                                                                         <  1:2324665/100‑1 (MQ=255)

                           ccTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCt                                                                         >  1:1108542/1‑100 (MQ=255)

                             tccgatctCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTg                                                                       <  2:362261/93‑1 (MQ=255)

                             tCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTg                                                                       <  1:1629898/100‑1 (MQ=255)

                             tCGTTAATCATAAAACAATAAAATTAAATATTCTCGAAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTg                                                                       >  1:1645641/1‑100 (MQ=255)

                              cGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGc                                                                      <  1:44007/100‑1 (MQ=255)

                                ttAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCaa                                                                    <  2:476183/100‑1 (MQ=255)

                                    tCATAAAACAATAAAATGAAATATTCTCGCAGTATATGACAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCGTTTTAGACATATCATTGCAAttaa                                                                >  2:1475608/1‑97 (MQ=255)

                                     cATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTgagatcgga                                                               >  1:362264/1‑92 (MQ=255)

                                     cATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTaa                                                               <  1:2064827/100‑1 (MQ=255)

                                     cATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTaa                                                               >  2:6521/1‑100 (MQ=255)

                                     cATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTaa                                                               >  1:111789/1‑100 (MQ=255)

                                      ntAAAACAATAAAATTAAATATTCTCGCAGTATATAAAAGTCTAAAGCATAAAAAATTTGATCAACATCTTTAATTTTAGACATCTACTTGCAATATaaa                                                              >  2:2178456/2‑100 (MQ=255)

                                      aTAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTaaa                                                              >  2:902458/1‑100 (MQ=255)

                                      aTAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTaaa                                                              >  1:546392/1‑100 (MQ=255)

                                      aTAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTaaa                                                              >  1:155134/1‑100 (MQ=255)

                                      aTAAAACAATAAAATTAAATATTCACGCAGTATATGACAATCTAAAGAATAAAAGATTTGATAAACATCTTTCATTTTAGACATCTCCTTGCAATATaaa                                                              >  2:2178128/1‑100 (MQ=255)

                                       tAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTaaaa                                                             >  1:657559/1‑100 (MQ=255)

                                        aaaaCAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAat                                                            >  1:1981161/1‑100 (MQ=255)

                                         aaaCAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAata                                                           >  1:1589244/1‑100 (MQ=255)

                                         aaaCAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAata                                                           >  1:2326725/1‑100 (MQ=255)

                                         aaaCAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCATTGCAATGTataata                                                           >  2:1915586/1‑100 (MQ=255)

                                          aaCAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAatat                                                          <  1:623806/100‑1 (MQ=255)

                                            cAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCa                                                        <  1:385377/100‑1 (MQ=255)

                                              aTAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATg                                                      <  2:1101694/100‑1 (MQ=255)

                                              aTAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATg                                                      >  1:688322/1‑100 (MQ=255)

                                              aTAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATg                                                      >  2:2207119/1‑100 (MQ=255)

                                                 aaaTTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTAt                                                   >  1:1198414/1‑100 (MQ=255)

                                                     tAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATcaac                                               <  1:1726052/100‑1 (MQ=255)

                                                      aaaTATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATcaaca                                              >  2:2333663/1‑100 (MQ=255)

                                                      aaaTATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATcaaca                                              <  2:2201555/100‑1 (MQ=255)

                                                      aaaTATTCTCGCAGTATATGACAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATcaaca                                              <  1:2341379/100‑1 (MQ=255)

                                                             ctcGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACAtat                                        >  1:617216/1‑97 (MQ=255)

                                                                gCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTgg                                      >  2:15612/1‑98 (MQ=255)

                                                                gCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTg                                    <  2:1044145/100‑1 (MQ=255)

                                                                gCAGTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTg                                    >  1:2395090/1‑100 (MQ=255)

                                                                 cAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGc                                   >  2:1004932/1‑100 (MQ=255)

                                                                   gTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCaa                                 >  1:1717168/1‑100 (MQ=255)

                                                                   gTATATGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCaa                                 >  1:354886/1‑100 (MQ=255)

                                                                     atatGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAAt                               <  2:2005925/100‑1 (MQ=255)

                                                                      tatGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTTGTGCAAATc                              <  1:1362495/100‑1 (MQ=255)

                                                                      tatGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATc                              >  1:1020077/1‑100 (MQ=255)

                                                                      tatGGCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATc                              >  1:418832/1‑100 (MQ=255)

                                                                        tGGCGGTCTAAAGCATCCAAAATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAAAGTAAAATATCATGTATCAAACACATCTGGTGCAAAACCa                            >  2:442699/1‑100 (MQ=255)

                                                                        tGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCa                            >  2:1692447/1‑100 (MQ=255)

                                                                        tGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCa                            >  2:286880/1‑100 (MQ=255)

                                                                         ngCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCAt                           >  2:1214827/2‑100 (MQ=255)

                                                                         ggCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCAt                           <  1:1912610/100‑1 (MQ=255)

                                                                          gCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATTTCATGTATCAACAACATCTGGTGCAAATCCATg                          <  2:1579032/100‑1 (MQ=255)

                                                                          gCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATg                          <  1:796559/100‑1 (MQ=255)

                                                                          gCAGTCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATg                          <  1:1193792/100‑1 (MQ=255)

                                                                              tCTAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTc                      <  2:1897994/100‑1 (MQ=255)

                                                                               cTAAAGAATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTTGTGCAAATCCATGAGTCg                     <  2:737134/100‑1 (MQ=255)

                                                                                tAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGa                    >  1:1520612/1‑100 (MQ=255)

                                                                                tAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGa                    <  1:815954/100‑1 (MQ=255)

                                                                                tAAAGCATCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGa                    >  1:1752201/1‑100 (MQ=255)

                                                                                 aaaGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGaa                   >  2:536551/1‑100 (MQ=255)

                                                                                   aGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACt                 >  1:1491408/1‑100 (MQ=255)

                                                                                    gCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTc                >  2:1568031/1‑100 (MQ=255)

                                                                                     cATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCg               >  1:2256561/1‑100 (MQ=255)

                                                                                     cATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCg               >  2:846122/1‑100 (MQ=255)

                                                                                      aTCAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGa              >  1:1397010/1‑100 (MQ=255)

                                                                                       tCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGAc             >  2:1301183/1‑100 (MQ=255)

                                                                                        cAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACa            >  1:2140795/1‑100 (MQ=255)

                                                                                        cAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACa            >  2:1010452/1‑100 (MQ=255)

                                                                                        cAAAGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACa            >  2:2100924/1‑100 (MQ=255)

                                                                                         aaaGATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACAt           <  2:804148/100‑1 (MQ=255)

                                                                                            gATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCa        <  2:1607280/100‑1 (MQ=255)

                                                                                            gATTTGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCa        >  1:1791859/1‑100 (MQ=255)

                                                                                              tttGAACAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATc      >  1:1388435/1‑100 (MQ=255)

                                                                                                tGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATCtt    <  1:1010463/100‑1 (MQ=255)

                                                                                                tGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATCtt    <  1:2160052/100‑1 (MQ=255)

                                                                                                tGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATCtt    >  2:1366221/1‑100 (MQ=255)

                                                                                                tGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATCtt    >  2:1746430/1‑100 (MQ=255)

                                                                                                  aTCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATCTTAg  >  1:1944113/1‑100 (MQ=255)

ACTAAAATAGATAACAAAAAGATGGAACCTCGTTAATCATAAAACAATAAAATTAAATATTCTCGCAGTATATGGCAGTCTAAAGCATCAAAGATTTGATCAACATCTTTCATTTTAGACATCTCCTTGCAATGTAAAATATCATGTATCAACAACATCTGGTGCAAATCCATGAGTCGAACTCGACATTCATCTTAG  >  NC_000913/2989787‑2989984

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