BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,466,501	G→T	G75G (GGC→GGA)	bfr ←	bacterioferritin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,466,501	0	G	T	98.2%	382.1 / 0.2	109	G75G (GGC→GGA)	bfr	bacterioferritin
Reads supporting (aligned to +/- strand): ref base G (1/1); new base T (42/65); total (43/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 110 total aligned reads displayed.

TCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACA  >  NC_000913/3466365‑3466634

tcatcaTATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTccc                                                                                                                                    <  2:16055/140‑1 (MQ=255)

  atcatATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAg                                                                                                                                  <  2:1091606/140‑1 (MQ=255)

  atcatATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAg                                                                                                                                  <  1:1684296/140‑1 (MQ=255)

  atcatATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAg                                                                                                                                  <  2:322264/140‑1 (MQ=255)

  atcatATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAg                                                                                                                                  >  2:1657345/1‑140 (MQ=255)

  atcatATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAg                                                                                                                                  >  2:1629096/1‑140 (MQ=255)

   tcatATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAgg                                                                                                                                 <  1:1629096/140‑1 (MQ=255)

   tcatATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAgg                                                                                                                                 <  1:1552823/140‑1 (MQ=255)

    catATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACACAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGt                                                                                                                                >  2:31686/1‑140 (MQ=255)

      tatCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTcc                                                                                                                              >  1:674302/1‑140 (MQ=255)

       atCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCt                                                                                                                             <  1:1391870/140‑1 (MQ=255)

         cgcgGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGt                                                                                                                           <  2:1281921/140‑1 (MQ=255)

            ggCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAg                                                                                                                        >  2:1093948/1‑140 (MQ=255)

            ggCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAg                                                                                                                        >  2:1176919/1‑140 (MQ=255)

             gCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGt                                                                                                                       >  1:978064/1‑140 (MQ=255)

             gCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGt                                                                                                                       >  2:34072/1‑140 (MQ=255)

             gCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGt                                                                                                                       >  2:402809/1‑140 (MQ=255)

              cTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGtt                                                                                                                      <  2:1246584/140‑1 (MQ=255)

                 aCGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTgg                                                                                                                   >  2:289331/1‑140 (MQ=255)

                 aCGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTgg                                                                                                                   >  1:346849/1‑140 (MQ=255)

                  cGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGa                                                                                                                  <  2:1166822/140‑1 (MQ=255)

                  cGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGa                                                                                                                  <  1:1176919/140‑1 (MQ=255)

                  cGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACACAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGa                                                                                                                  <  1:31686/140‑1 (MQ=255)

                        cATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCt                                                                                                            >  1:262453/1‑140 (MQ=255)

                          tGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTc                                                                                                          <  1:823465/140‑1 (MQ=255)

                          tGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTc                                                                                                          <  2:60240/140‑1 (MQ=255)

                          tGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTc                                                                                                          <  1:531237/140‑1 (MQ=255)

                           gAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTcc                                                                                                         >  2:964676/1‑140 (MQ=255)

                            aaCGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCa                                                                                                        >  1:947113/1‑140 (MQ=255)

                            aaCGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCa                                                                                                        >  1:196661/1‑140 (MQ=255)

                                 tATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGaaaa                                                                                                   >  2:1429629/1‑140 (MQ=255)

                                   tCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAg                                                                                                 >  1:1703051/1‑140 (MQ=255)

                                    cGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGa                                                                                                <  1:94931/140‑1 (MQ=255)

                                      gCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACTTTCCAGAAAAAGAAt                                                                                              <  1:1119887/140‑1 (MQ=255)

                                      gCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAAt                                                                                              <  1:34072/140‑1 (MQ=255)

                                      gCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAAt                                                                                              <  2:1005622/140‑1 (MQ=255)

                                      gCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAAt                                                                                              <  1:289331/140‑1 (MQ=255)

                                         tAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgcg                                                                                           <  2:1098195/140‑1 (MQ=255)

                                         tAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgcg                                                                                           <  1:158970/140‑1 (MQ=255)

                                         tAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgcg                                                                                           <  2:843686/140‑1 (MQ=255)

                                         tAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgcg                                                                                           <  1:417541/140‑1 (MQ=255)

                                         tAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgc                                                                                            >  1:115866/1‑139 (MQ=255)

                                         tAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgc                                                                                            <  2:115866/139‑1 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  1:744191/140‑1 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  1:1639793/140‑1 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  2:923742/140‑1 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  1:807038/140‑1 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  2:346849/140‑1 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         >  2:1551922/1‑140 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  2:164260/140‑1 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         >  1:561170/1‑140 (MQ=255)

                                           aCCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  1:1201448/140‑1 (MQ=255)

                                                    cTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAAc                                                                                <  1:1346111/140‑1 (MQ=255)

                                                      cACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         >  1:1177894/1‑129 (MQ=255)

                                                      cACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                         <  2:1177894/129‑1 (MQ=255)

                                                      cACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGa                                                                              >  2:296867/1‑140 (MQ=255)

                                                          cAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCgg                                                                          <  2:912626/140‑1 (MQ=255)

                                                          cAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCgg                                                                          <  1:1429629/140‑1 (MQ=255)

                                                           aaaTTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGc                                                                         <  1:578047/140‑1 (MQ=255)

                                                           aaaTTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGc                                                                         <  2:1243767/140‑1 (MQ=255)

                                                                    gcgcCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcat                                                                <  1:599740/140‑1 (MQ=255)

                                                                     cgcCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatc                                                               >  1:454106/1‑140 (MQ=255)

                                                                      gcCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatct                                                              >  1:1269691/1‑140 (MQ=255)

                                                                        cATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatctca                                                            <  1:1018333/140‑1 (MQ=255)

                                                                         aTCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatctcat                                                           <  1:964676/140‑1 (MQ=255)

                                                                          tCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatctcatc                                                          <  1:50909/140‑1 (MQ=255)

                                                                            cAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCaa                                                        <  2:1165374/140‑1 (MQ=255)

                                                                            cAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCaa                                                        <  2:1714673/140‑1 (MQ=255)

                                                                             aGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAAt                                                       <  1:1093948/140‑1 (MQ=255)

                                                                             aGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAAt                                                       <  1:1551922/140‑1 (MQ=255)

                                                                               cTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATgg                                                     <  2:947113/140‑1 (MQ=255)

                                                                                    gTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCtcatca                                                         <  1:182884/131‑1 (MQ=255)

                                                                                    gTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCtcatca                                                         >  2:182884/1‑131 (MQ=255)

                                                                                    gTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCa                                                <  1:1180587/140‑1 (MQ=255)

                                                                                    gTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCa                                                <  1:402809/140‑1 (MQ=255)

                                                                                      gCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATg                                              <  2:482783/140‑1 (MQ=255)

                                                                                       ccAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGa                                             <  2:674302/140‑1 (MQ=255)

                                                                                       ccAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGa                                             <  2:1703051/140‑1 (MQ=255)

                                                                                         aGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATa                                           <  1:42925/140‑1 (MQ=255)

                                                                                           aTCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACt                                         <  1:321595/140‑1 (MQ=255)

                                                                                             cAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCGGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTcc                                       <  1:1241147/140‑1 (MQ=255)

                                                                                              aGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCca                                      <  2:1090467/140‑1 (MQ=255)

                                                                                                aaCGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGa                                             >  1:579635/1‑131 (MQ=255)

                                                                                                aaCGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGa                                             <  2:579635/131‑1 (MQ=255)

                                                                                                  cGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCAcatc                                  >  1:898279/1‑140 (MQ=255)

                                                                                                    cagcaTTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCATGTAAGTTTGGAAGACCTTCCAGAAAAAGAATACGCTCAATATATCGAGCGGAGTGTTAAGTGTCATCAAGGGATTCATGATACTCAAcatcat                                <  2:1651087/140‑1 (MQ=255)

                                                                                                     agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACt                                         >  2:1725817/1‑130 (MQ=255)

                                                                                                     agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACt                                         <  1:1725817/130‑1 (MQ=255)

                                                                                                     agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGCCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCTCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACCCCACATCAtt                               >  1:487702/1‑140 (MQ=255)

                                                                                                     agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGa                                             >  2:816571/1‑126 (MQ=255)

                                                                                                     agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGa                                             <  1:816571/126‑1 (MQ=255)

                                                                                                       caTTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATGATTga                             >  2:1062330/1‑140 (MQ=255)

                                                                                                        aTTTCCTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTgag                            <  2:1269691/140‑1 (MQ=255)

                                                                                                            ccTCAACATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGt                        >  2:403241/1‑140 (MQ=255)

                                                                                                                 aCATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTgag                   >  2:590776/1‑140 (MQ=255)

                                                                                                                  cATCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTgaga                  >  1:267394/1‑140 (MQ=255)

                                                                                                                   aTCTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGAc                 >  1:1394519/1‑140 (MQ=255)

                                                                                                                     cTTCACCAATGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGAccc               >  1:1113585/1‑140 (MQ=255)

                                                                                                                              tGTTCAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTa      >  2:456853/1‑140 (MQ=255)

                                                                                                                                  cAGTTTTCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACa  >  1:454966/1‑140 (MQ=255)

TCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACA  >  NC_000913/3466365‑3466634

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