BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,473,032	A→G	K19R (AAG→AGG)	b2361 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,473,032	0	A	G	91.4%	1055.6 / 69.7	395	K19R (AAG→AGG)	b2361	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (15/19); new base G (188/173); total (203/192)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 398 total aligned reads displayed.

ATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTC  >  NC_000913/2472942‑2473121

 tGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGct                                                                                        >  2:2922661/1‑92 (MQ=255)

 tGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCa                                                                                        >  1:6249449/1‑93 (MQ=255)

 tGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCa                                                                                        >  1:8160169/1‑93 (MQ=255)

 tGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCa                                                                                        <  2:7723531/93‑1 (MQ=255)

  ggTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCAt                                                                                       <  2:839611/93‑1 (MQ=255)

  ggTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGAGAGGCAt                                                                                       >  1:9317358/1‑93 (MQ=255)

   gTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCAtt                                                                                      <  1:5036052/93‑1 (MQ=255)

    tGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCAttt                                                                                     <  1:1334606/93‑1 (MQ=255)

     gAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCAtttt                                                                                    <  1:5232354/93‑1 (MQ=255)

     gAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCAtttt                                                                                    >  1:7775947/1‑93 (MQ=255)

       gggCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTtta                                                                                  >  2:118891/1‑93 (MQ=255)

         gCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTttatt                                                                                >  1:2991153/1‑93 (MQ=255)

          cTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTttatta                                                                               <  2:3480485/93‑1 (MQ=255)

          cTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTttatta                                                                               <  2:191797/93‑1 (MQ=255)

          cTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTttatta                                                                               <  2:6165509/93‑1 (MQ=255)

           ttCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACATAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTtattat                                                                              >  2:4467818/1‑93 (MQ=255)

           ttCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTtattat                                                                              <  2:7018168/93‑1 (MQ=255)

           ttCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTtattat                                                                              <  1:3564395/93‑1 (MQ=255)

           ttCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTtattat                                                                              <  1:1483361/93‑1 (MQ=255)

            tCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATg                                                                             >  2:2453935/1‑93 (MQ=255)

            tCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATg                                                                             >  2:7602887/1‑93 (MQ=255)

            tCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATg                                                                             <  1:6887602/93‑1 (MQ=255)

            tCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATg                                                                             >  2:7387844/1‑93 (MQ=255)

            tCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATg                                                                             >  2:6840243/1‑93 (MQ=255)

             ccTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGa                                                                            <  1:6651599/93‑1 (MQ=255)

             ccTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGa                                                                            >  1:978885/1‑93 (MQ=255)

             ccTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGa                                                                            >  1:3724045/1‑93 (MQ=255)

              cTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGAc                                                                           >  2:6416439/1‑93 (MQ=255)

              cTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGAc                                                                           >  2:8591440/1‑93 (MQ=255)

              cTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGAc                                                                           >  2:3150169/1‑93 (MQ=255)

               ttttACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACa                                                                          >  1:3158280/1‑93 (MQ=255)

               ttttACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACa                                                                          >  1:8328784/1‑93 (MQ=255)

                tttACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACaa                                                                         <  1:792583/93‑1 (MQ=255)

                tttACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACaa                                                                         <  1:4066682/93‑1 (MQ=255)

                 ttACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAg                                                                        <  1:9051290/93‑1 (MQ=255)

                  tACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGa                                                                       <  2:3015025/93‑1 (MQ=255)

                     ccAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGataa                                                                    >  1:7578064/1‑93 (MQ=255)

                      cAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGataaa                                                                   <  1:4428029/93‑1 (MQ=255)

                      cAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGataaa                                                                   <  2:9224487/93‑1 (MQ=255)

                      cAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGataaa                                                                   <  2:345520/93‑1 (MQ=255)

                       aaaCTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATAACAAGataaat                                                                  <  1:7393860/93‑1 (MQ=255)

                        aaCTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGataaata                                                                 <  2:2073170/93‑1 (MQ=255)

                        aaCTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGataaata                                                                 <  1:437067/93‑1 (MQ=255)

                        aaCTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGataaata                                                                 <  2:4266800/93‑1 (MQ=255)

                        aaCTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGataaata                                                                 <  2:5039994/93‑1 (MQ=255)

                         aCTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGataaataa                                                                <  2:4073577/93‑1 (MQ=255)

                           tCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAg                                                              <  2:324903/93‑1 (MQ=255)

                            cAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAga                                                             <  2:2677348/93‑1 (MQ=255)

                            cAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAga                                                             <  2:5429321/93‑1 (MQ=255)

                               cgcgGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAgacga                                                          <  1:9638119/93‑1 (MQ=255)

                               cgcgGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAgacga                                                          >  2:5644185/1‑93 (MQ=255)

                                gcgGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAgacgac                                                         <  1:2804279/93‑1 (MQ=255)

                                 cgGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACa                                                        >  1:5896948/1‑93 (MQ=255)

                                  gtggCAGGGANNATNCGGGGAGAACAAAATGTNATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGANAt                                                       <  2:4706729/90‑1 (MQ=255)

                                   gTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATc                                                      >  2:62065/1‑93 (MQ=255)

                                   gTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATc                                                      >  2:1752998/1‑93 (MQ=255)

                                     gCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGt                                                    <  1:8642040/93‑1 (MQ=255)

                                       aGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGa                                                  >  1:989162/1‑93 (MQ=255)

                                       aGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGa                                                  <  1:2025316/93‑1 (MQ=255)

                                          gcATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTa                                               >  2:9520057/1‑93 (MQ=255)

                                           cATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGCCAAGATAAATAAAGACGACATCGTGATTaa                                              >  1:5548680/1‑93 (MQ=255)

                                           cATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTaa                                              >  1:3201756/1‑93 (MQ=255)

                                           cATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTaa                                              >  1:1194487/1‑93 (MQ=255)

                                           cATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTaa                                              >  1:7131282/1‑93 (MQ=255)

                                            atatACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAAc                                             >  1:2543653/1‑93 (MQ=255)

                                             tataCGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACg                                            >  1:4250259/1‑93 (MQ=255)

                                             tataCGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACg                                            >  1:6013271/1‑93 (MQ=255)

                                             tataCGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACg                                            >  2:4900047/1‑93 (MQ=255)

                                                aCGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGata                                         <  1:8553537/93‑1 (MQ=255)

                                                aCGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGata                                         <  1:3464099/93‑1 (MQ=255)

                                                  gTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATc                                       <  2:6528148/93‑1 (MQ=255)

                                                  gTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATc                                       >  2:5178580/1‑93 (MQ=255)

                                                  gTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATc                                       <  2:8721198/93‑1 (MQ=255)

                                                  gTGGAGAACAAAATGTCATGTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATc                                       <  2:5295565/93‑1 (MQ=255)

                                                   tGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATcg                                      <  2:9703142/93‑1 (MQ=255)

                                                      agaACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCgg                                   >  2:8828057/1‑93 (MQ=255)

                                                      agaACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCgg                                   >  2:103709/1‑93 (MQ=255)

                                                      agaACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCgg                                   >  1:4095365/1‑93 (MQ=255)

                                                       gaACAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGCGATTAACGATATCGCGGt                                  >  1:951146/1‑93 (MQ=255)

                                                        aaCAAAATGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGtt                                 <  2:8062816/93‑1 (MQ=255)

                                                            aaaTGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAACGACGACATCGTGATTACCGATATCGCGGTTTccc                             >  2:1076661/1‑93 (MQ=255)

                                                              aTGTCATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCtt                           <  2:7849302/93‑1 (MQ=255)

                                                                  cATTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCaa                       <  1:2601250/93‑1 (MQ=255)

                                                                   atttatttAAACTTTTTCCAGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATAGCGCGGTTTCCCTTTCaaa                      >  1:6432877/1‑93 (MQ=255)

                                                                   aTTTATTAAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCaaa                      <  1:2400018/93‑1 (MQ=255)

                                                                     ttattaAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAata                    <  2:1365928/93‑1 (MQ=255)

                                                                       attaAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATc                  <  2:161928/93‑1 (MQ=255)

                                                                         taAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTg                >  2:5752881/1‑93 (MQ=255)

                                                                         taAAACTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTg                <  2:9415868/93‑1 (MQ=255)

                                                                          aaaaCTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGc               >  2:9089346/1‑93 (MQ=255)

                                                                           aaaCTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGcc              >  1:4706679/1‑93 (MQ=255)

                                                                           aaaCTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGcc              <  1:9169571/93‑1 (MQ=255)

                                                                            aaCTTTTTCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCg             >  2:3727029/1‑93 (MQ=255)

                                                                               tttttCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCtt          >  1:8864361/1‑93 (MQ=255)

                                                                               tttttCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGGGATTAAGGATATCGCGGTTTCCCTTTCAAATACCTGCCGCtt          >  1:1085699/1‑93 (MQ=255)

                                                                                  ttCCGGGAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGc       <  2:5900278/93‑1 (MQ=255)

                                                                                      gggAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGt   >  2:6073409/1‑93 (MQ=255)

                                                                                       ggAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCCAATATCTGCCGCTTTGCCGGTc  <  1:1825071/93‑1 (MQ=255)

                                                                                       ggAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTc  <  2:8902955/93‑1 (MQ=255)

                                                                                       ggAGGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTc  <  1:575495/93‑1 (MQ=255)

ATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTC  >  NC_000913/2472942‑2473121

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