BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,573,189	T→C	17.8%	pseudogene (776/1503 nt)	yjiT →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,189	0	T	C	17.8%	255.1 / 48.7	120	pseudogene (776/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand): ref base T (46/52); new base C (12/9); total (58/62)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

Only 100 of 122 total aligned reads displayed.

TTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTT  >  NC_000913/4573090‑4573286

ttGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                   <  1:2095727/100‑1 (MQ=255)

 tGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTa                                                                                                  <  1:2384821/100‑1 (MQ=255)

    gaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                                                               >  2:2192416/1‑100 (MQ=255)

      aGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGc                                                                                             <  1:987804/100‑1 (MQ=255)

       ggTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCt                                                                                            >  1:2403130/1‑100 (MQ=255)

       ggTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCt                                                                                            >  2:424743/1‑100 (MQ=255)

         ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                          <  1:743156/100‑1 (MQ=255)

         ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                          >  1:960691/1‑100 (MQ=255)

         ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgc                                                                                          <  1:1955463/100‑3 (MQ=255)

         tgcACCCTACCTGACGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                          >  1:1727459/3‑100 (MQ=255)

          tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                         <  2:769631/100‑1 (MQ=255)

          tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:602995/100‑4 (MQ=255)

          tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2317231/100‑4 (MQ=255)

          tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2316714/100‑4 (MQ=255)

          tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2231539/100‑4 (MQ=255)

          tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2214680/100‑4 (MQ=255)

           cACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                        >  1:1841131/1‑100 (MQ=255)

           cACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                        <  2:2540445/100‑1 (MQ=255)

            aCCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGcc                                                                                       >  1:2489445/1‑100 (MQ=255)

             cccTACCTTGCGACGCGAGATTCAAAAACCGGTAGAATATCAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGACTAGGcc                                                                                       >  2:456915/1‑99 (MQ=255)

              ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATGTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgcata                                                                                     <  1:1474808/100‑8 (MQ=255)

               cTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                    <  1:1508372/100‑1 (MQ=255)

               cTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                    >  2:1835159/1‑100 (MQ=255)

                tACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgcatagc                                                                                   <  1:673469/100‑10 (MQ=255)

                  ccTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACt                                                                                 >  2:1191932/1‑100 (MQ=255)

                    tGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTAc                                                                               <  2:1177803/100‑1 (MQ=255)

                         aaGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaa                                                                          <  2:2275467/100‑1 (MQ=255)

                            cGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAAt                                                                       >  2:2348631/1‑100 (MQ=255)

                                 ttCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATCagag                                                                  >  2:238160/1‑100 (MQ=255)

                                   ccAATACCGTTAGATGATGAAACAGGAACGCATTTTTTTAATGGACTGTTATGCGCAGCCACTGTAGATGCTAAGGCTAGACTACAAAAAAATAAGAGt                                                                 >  2:555069/1‑99 (MQ=255)

                                    cAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTAGACTACAAAAAAATAAGAGTAc                                                               <  2:1489424/100‑1 (MQ=255)

                                     aaTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcg                                                              >  1:109175/1‑100 (MQ=255)

                                      aTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgc                                                             >  2:1834482/1‑100 (MQ=255)

                                          cGTTAGATGATGAAACAGGAACTCATGTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                         <  1:853260/100‑1 (MQ=255)

                                          cGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                         <  1:860868/100‑1 (MQ=255)

                                          cGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                         >  2:1561641/1‑100 (MQ=255)

                                           gTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                        >  1:1381395/1‑100 (MQ=255)

                                             tAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCa                                                      >  1:1133274/1‑100 (MQ=255)

                                              aGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                     >  2:2185352/1‑100 (MQ=255)

                                               gatgatGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAt                                                    >  2:2043126/1‑100 (MQ=255)

                                                atgatgAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                   >  2:1426653/1‑100 (MQ=255)

                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  <  1:1165942/100‑1 (MQ=255)

                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  >  1:1254679/1‑100 (MQ=255)

                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  >  1:1933007/1‑100 (MQ=255)

                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  <  1:324349/100‑1 (MQ=255)

                                                     gaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                              <  1:592886/100‑1 (MQ=255)

                                                     gaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                              <  2:769639/100‑1 (MQ=255)

                                                      aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg                                             <  1:2211135/100‑1 (MQ=255)

                                                      aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg                                             <  1:2329601/100‑1 (MQ=255)

                                                      aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg                                             >  2:386063/1‑100 (MQ=255)

                                                        aCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGt                                           >  2:1790264/1‑100 (MQ=255)

                                                         cAGTAACGCATTTCATGAATTGACTGCTATACGCTGCCTCTGTAGAGGCTAAGCATCAACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGtt                                          <  2:129968/100‑2 (MQ=255)

                                                         cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                          <  2:110699/100‑1 (MQ=255)

                                                         cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                          >  1:2034501/1‑100 (MQ=255)

                                                         cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                          <  1:946334/100‑1 (MQ=255)

                                                          aGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCa                                         <  1:283427/100‑1 (MQ=255)

                                                           ggAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAg                                        <  1:1881014/100‑1 (MQ=255)

                                                           ggAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAg                                        >  2:1031027/1‑100 (MQ=255)

                                                             aaCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaa                                      >  1:95732/1‑100 (MQ=255)

                                                              aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaa                                     >  2:840121/1‑100 (MQ=255)

                                                              aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaa                                     <  1:134161/100‑1 (MQ=255)

                                                               cGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                    <  2:727451/100‑1 (MQ=255)

                                                                 cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGc                                  >  1:1585397/1‑100 (MQ=255)

                                                                   ttttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAt                                >  1:800168/1‑100 (MQ=255)

                                                                     ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAACAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATcc                              >  2:1226465/1‑100 (MQ=255)

                                                                     ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCCGAAAAGCATcc                              <  2:1967115/100‑1 (MQ=255)

                                                                      tttGAATGGACTGCTATGCGCAGCCTCTTTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCa                             <  2:324874/100‑1 (MQ=255)

                                                                      tttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCa                             >  1:1444316/1‑100 (MQ=255)

                                                                       ttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAg                            >  2:2025751/1‑100 (MQ=255)

                                                                         gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                          >  1:571591/1‑100 (MQ=255)

                                                                         gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                          >  1:1233303/1‑100 (MQ=255)

                                                                         gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGCCAGAAAAGCATCCAGAt                          >  1:2005827/1‑100 (MQ=255)

                                                                          aaTGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGTGTACGCGTTGCCAATTCTATTGGTCAGACAAGCATCCAGacg                         <  1:1737074/100‑3 (MQ=255)

                                                                            tGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAg                       <  1:1043333/100‑1 (MQ=255)

                                                                             ggACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGc                      <  1:61065/100‑1 (MQ=255)

                                                                               aCTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCtt                    <  1:424747/100‑1 (MQ=255)

                                                                                cTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTa                   <  2:1257616/100‑1 (MQ=255)

                                                                                 gtcgaggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAg                  >  2:306769/11‑100 (MQ=38)

                                                                                   cTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAggg                <  2:437852/100‑1 (MQ=255)

                                                                                    tATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt               >  1:2303735/1‑100 (MQ=255)

                                                                                     aggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt               >  1:1451532/7‑99 (MQ=255)

                                                                                     aggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTg              >  1:2358617/7‑100 (MQ=255)

                                                                                     aTGCGCAGCCTCTGAAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt               <  1:1015553/99‑1 (MQ=255)

                                                                                      tGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTg              >  1:1122795/1‑99 (MQ=255)

                                                                                      tGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa             >  1:373783/1‑100 (MQ=255)

                                                                                      tGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa             <  1:2134458/100‑1 (MQ=255)

                                                                                       gcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGAt            >  2:1954447/5‑100 (MQ=255)

                                                                                        cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATCAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATCAGCTTAGGGTGATa           >  2:16231/4‑100 (MQ=255)

                                                                                        cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  1:1757760/4‑100 (MQ=255)

                                                                                        cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  2:5043/4‑100 (MQ=255)

                                                                                        cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  1:1409087/4‑100 (MQ=255)

                                                                                        cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  1:1616619/4‑100 (MQ=255)

                                                                                        cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAATACGCGTTGCCAATTATATTGATCAGAAAAGCATCCAAATAAGCTTAggaaaata           >  2:12718/4‑94 (MQ=255)

                                                                                         gcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           <  2:2034487/99‑1 (MQ=255)

                                                                                         gcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAg          <  1:1526316/100‑1 (MQ=255)

                                                                                           tctctctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTa        <  1:1562020/97‑1 (MQ=255)

                                                                                              ctctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTaaga     >  2:1562009/1‑97 (MQ=255)

                                                                                                ctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATctct   <  2:1333093/100‑1 (MQ=255)

                                                                                                 tGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCtt  >  1:2382758/1‑100 (MQ=255)

                                                                                                   tAGAGGCTAAGCCACGACTACAAAAAACTCAGAGTACGCGGTGCCAATTCTATTGGACAGAAACGCATCCAGATGAGCTTAGGGTGATAGTATATCtt  >  2:1077859/1‑98 (MQ=255)

TTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTT  >  NC_000913/4573090‑4573286

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