BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,802	A→G	100%	V222A (GTA→GCA)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	953,802	0	A	G	98.7%	487.2 / ‑3.1	149	V222A (GTA→GCA)	focA	formate channel
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (73/73); total (74/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

Only 100 of 149 total aligned reads displayed.

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTG  >  NC_000913/953705‑953897

  tgaccctgatgaATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCCAAATTCCCGGGATGCGAAGTCGCGGATTGCAAt                                                                                             >  2:195689/7‑100 (MQ=255)

  tATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAt                                                                                             <  1:2673679/100‑1 (MQ=255)

   gaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATa                                                                                            >  1:652712/5‑100 (MQ=255)

   gaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCAGTCCAAAATTCCGAAGATGCGAAGTCGCGAATTGCAATa                                                                                            >  1:786888/5‑100 (MQ=255)

    tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGTAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                           <  1:3241808/100‑1 (MQ=255)

    tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                           >  2:711741/1‑100 (MQ=255)

    tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                           >  2:2772549/1‑100 (MQ=255)

     ccgTGATGAATTTCATCACGGTCAGGTGAGAACAATTTTCCGGTTCAGAACCGACTGCGTTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          >  1:1753548/3‑100 (MQ=255)

     ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          >  1:1656526/3‑100 (MQ=255)

     ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          >  2:431498/3‑100 (MQ=255)

     ccgTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          >  1:344501/3‑100 (MQ=255)

     cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                           <  2:3853634/99‑1 (MQ=255)

      aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                         >  1:1158904/1‑100 (MQ=255)

      aGTGATGAAATTCATCACGGTCAGGAGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                         >  2:1645479/1‑100 (MQ=255)

         gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATc                                                                                      >  2:3092874/1‑100 (MQ=255)

           tgaAATTCCTCACGGTCAGGTGAAAAAAATTTTCCGGTGCAGAACCGACTTCGGTCCAAAATTCCGGAGATGCGAAGTCGCGGATTGCAATACCCAACgg                                                                                    >  2:3785243/1‑100 (MQ=255)

            gaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCggg                                                                                   >  2:849833/1‑100 (MQ=255)

             aatttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTACGGTCCAAAATTCCGGGAATGCGAAGTCGCGGATTACAATACCCATCGGGa                                                                                  >  1:550028/4‑100 (MQ=255)

             aaaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGa                                                                                  <  2:2211757/100‑1 (MQ=255)

              aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                                 >  1:1288786/1‑100 (MQ=255)

               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                                <  1:2357140/100‑1 (MQ=255)

               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATc                                                                                <  1:649917/100‑1 (MQ=255)

                  catcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATa                                                                             >  1:1150427/1‑100 (MQ=255)

                     caCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                          <  2:609445/100‑1 (MQ=255)

                     caCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACAGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                          <  1:960452/100‑1 (MQ=255)

                          tCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGtt                                                                     <  1:3631456/100‑1 (MQ=255)

                          tCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGtt                                                                     <  2:1998338/100‑1 (MQ=255)

                           cAGGTGAGAAAAATTATCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGttt                                                                    <  1:2899225/100‑1 (MQ=255)

                             ggtttgaAAAATGTTTCGGTGCATAACCTACTGCGGTCCGATATTGCGGGTATGCTAAGTCGCAGATTGCAATACCCGTCGGGATCATAAACATGTTTGc                                                                  <  2:3289661/95‑1 (MQ=255)

                             ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGACGATGCAAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                  >  1:368163/1‑100 (MQ=255)

                              gTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCg                                                                 <  2:2955694/100‑1 (MQ=255)

                              gTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCg                                                                 <  2:2865080/100‑1 (MQ=255)

                               tGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTACGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGa                                                                <  2:1362532/100‑1 (MQ=255)

                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             <  1:2984080/100‑1 (MQ=255)

                                   naaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACt                                                            >  1:607144/2‑100 (MQ=255)

                                   aaaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCTAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             <  2:919874/99‑1 (MQ=255)

                                   aaaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGAAATACACATCGGGATCATAAACATGTTTGCGATACt                                                            >  2:2465699/1‑100 (MQ=255)

                                   aaaaaTTTTCCGGTGCAGAACCGACTGCGGTACAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATagt                                                            >  2:2073049/1‑98 (MQ=255)

                                   aaaaaTTTTCCGGTGCAGAACCGACTGCGGTACAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATAGGGATCATAAACATGTTTGCGATACt                                                            >  2:1834959/1‑100 (MQ=255)

                                    aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGTGGATACGAAGAAGAAAATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                           >  1:193090/1‑100 (MQ=255)

                                    aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                           <  2:3151464/100‑1 (MQ=255)

                                    aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                           <  2:1175664/100‑1 (MQ=255)

                                     aaaTTTTCCTTTGCATAACCTACTGCGGTCCAAAATTCCGGTTATTCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtgt                                                          <  1:2703866/100‑1 (MQ=255)

                                       aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtgt                                                          >  2:3769589/1‑98 (MQ=255)

                                        ttttCCGGTGCAGAACCGACTGCTGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGTGATCATAAACATGTTTGCGATACTGTGc                                                        <  2:3801053/99‑1 (MQ=255)

                                          ttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCa                                                     <  1:3898207/100‑1 (MQ=255)

                                            ccGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                                                   <  1:2391992/100‑1 (MQ=255)

                                                tGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg                                               >  2:2655774/1‑100 (MQ=255)

                                                 gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                              >  2:3523322/1‑100 (MQ=255)

                                                 gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                              <  1:711750/100‑1 (MQ=255)

                                                   aGAACCGACTGCAGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                            <  1:280356/100‑1 (MQ=255)

                                                    gAACCGACTGCTGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                            <  1:2115633/99‑1 (MQ=255)

                                                    gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           <  2:3378937/100‑1 (MQ=255)

                                                    gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           <  2:3172109/100‑1 (MQ=255)

                                                    gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           <  1:3677499/100‑1 (MQ=255)

                                                     aaCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGTATTTCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                          <  1:1042131/100‑1 (MQ=255)

                                                      aCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCtcg                                           >  2:1197807/1‑96 (MQ=255)

                                                      aCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                          <  1:2517786/99‑1 (MQ=255)

                                                         gACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaaca                                      >  2:3203899/1‑100 (MQ=255)

                                                          aCTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaacaa                                     <  1:135168/100‑1 (MQ=255)

                                                           cTGCGGTCCAAAATTCCGGGGATACGAAGTCGCGGATTGCAATACCCATAGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa                                    >  2:862211/1‑100 (MQ=255)

                                                            tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGTATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                   <  1:2442466/100‑1 (MQ=255)

                                                            tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                   <  2:1162285/100‑1 (MQ=255)

                                                            tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                   <  2:2102237/100‑1 (MQ=255)

                                                             gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                                  >  2:346196/1‑100 (MQ=255)

                                                              cGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACAt                                 >  1:1111102/1‑100 (MQ=255)

                                                                gTCCAAAATTCCGGGGATGCGAAGTCGCTGATTGCAATACCCATCGGGATCATAAACATGTTTGCTATACTGTGCTCAAAACCGCAGACAACAAACATc                                >  1:705341/1‑99 (MQ=255)

                                                                gTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcg                               >  2:1648236/1‑100 (MQ=255)

                                                                  ccAAAATTCCGGGGATGCGAAGTCGTGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                             <  1:1865745/100‑1 (MQ=255)

                                                                  ccAAAATTCCGGGGATGCGAAGTCAAGGATTGCAATACCCATAGGGATCATAAACATGTTTGCGATACTGTGATAAAAACCGCTGACAACAAACATcgcg                             >  2:949659/1‑100 (MQ=255)

                                                                   cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                            >  1:3345581/1‑100 (MQ=255)

                                                                      aaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCg                         <  1:1497581/100‑1 (MQ=255)

                                                                      aaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCg                         <  1:3127907/100‑1 (MQ=255)

                                                                        ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                       <  2:28896/100‑1 (MQ=255)

                                                                          ccGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcag                     <  1:1623757/100‑1 (MQ=255)

                                                                          ccGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcag                     >  1:2374961/1‑100 (MQ=255)

                                                                           cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                    >  2:235163/1‑100 (MQ=255)

                                                                           cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                    >  2:245624/1‑100 (MQ=255)

                                                                           cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                    <  1:547194/100‑1 (MQ=255)

                                                                           cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                    >  1:428730/1‑100 (MQ=255)

                                                                            ggggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagca                   <  2:3490711/100‑1 (MQ=255)

                                                                                 tGCGAAGTCGCTGATTGCAATTCTCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCAc                  <  1:3191850/96‑1 (MQ=255)

                                                                                 tGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATa              >  2:1395889/1‑100 (MQ=255)

                                                                                  gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa             <  1:1235975/100‑1 (MQ=255)

                                                                                  gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa             >  2:2645523/1‑100 (MQ=255)

                                                                                  gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa             >  2:3653078/1‑100 (MQ=255)

                                                                                  gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGACAACAAACATCGCAACCGGCAACACCATaa             >  1:1095518/1‑100 (MQ=255)

                                                                                   cGAAGTCGCGGATTGCAATACCCATCGGGATCATAATCATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt            <  1:1979484/100‑1 (MQ=255)

                                                                                   cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt            >  1:96238/1‑100 (MQ=255)

                                                                                    gAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg           <  1:2284519/100‑1 (MQ=255)

                                                                                    gAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg           >  1:615832/1‑100 (MQ=255)

                                                                                       gTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAAc        >  1:3630825/1‑100 (MQ=255)

                                                                                       gTCGCGGATTGCAATACCCATCGGAATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACAACGCGAACGGCAGCACCAAAATGAAc        >  2:416502/1‑100 (MQ=255)

                                                                                        tCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACg       <  2:1753533/100‑1 (MQ=255)

                                                                                        tCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACg       <  1:3563864/100‑1 (MQ=255)

                                                                                         cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGaaag       >  1:2833543/1‑97 (MQ=255)

                                                                                         cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCACGACCAGCAGCACCATAATGAACGc      >  1:1551818/1‑100 (MQ=255)

                                                                                          gcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGAGACCAGCAGCACCATAATGAACGCt     >  1:1404089/1‑100 (MQ=255)

                                                                                           cgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCtt    >  1:1805505/1‑100 (MQ=255)

                                                                                             gATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTg  <  2:1327106/100‑1 (MQ=255)

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTG  >  NC_000913/953705‑953897

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