BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,880,375	G→T	5.6%	L283M (CTG→ATG) ‡	recF ←	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,880,375	0	G	T	5.6%	436.8 / 10.5	144	L283M (CTG→ATG) ‡	recF	gap repair protein
Reads supporting (aligned to +/- strand): ref base G (76/60); new base T (2/6); total (78/66)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 146 total aligned reads displayed.

ATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCA  >  NC_000913/3880227‑3880524

aTCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGa                                                                                                                                                      >  1:449650/1‑150 (MQ=255)

 tCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTACGCGGCTGGCAAGCAGCCCGCGACGCTCTTCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACTCCGCGGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGAc                                                                                                                                                     <  1:889828/150‑1 (MQ=255)

 tCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGAc                                                                                                                                                     <  2:185537/150‑1 (MQ=255)

  cgcgCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACg                                                                                                                                                    >  2:145549/1‑150 (MQ=255)

    cgcTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTa                                                                                                                                                  >  1:461446/1‑150 (MQ=255)

       tGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGTTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAgg                                                                                                                                               <  1:798335/150‑1 (MQ=255)

       tGACAAAGACCTGTGAGTGCGGCGCTTTTAAGCGGCGGGCAAGCAGCCCGCGACGCGCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAgg                                                                                                                                               <  2:1570044/150‑1 (MQ=255)

        gACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGgc                                                                                                                                              >  2:1712927/1‑150 (MQ=255)

          cAAAGCCCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGCGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCGTTCACGGGTGCGGACCTCTCCTTGCGCCAGACGTACGgcgc                                                                                                                                            >  1:71991/1‑150 (MQ=255)

           aaaGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGca                                                                                                                                           <  2:43405/150‑1 (MQ=255)

            aaGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGcac                                                                                                                                          >  2:826548/1‑150 (MQ=255)

            aaGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGcac                                                                                                                                          <  1:1773740/150‑1 (MQ=255)

             aGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCTTACGTCTGGCGcaag                                                                                                                                         <  1:1770720/150‑3 (MQ=255)

             aGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCTTACGTCTGGCGcaag                                                                                                                                         <  1:1700518/150‑3 (MQ=255)

              gACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCTTACGTCTGGCGcaagg                                                                                                                                        <  1:1316477/150‑4 (MQ=255)

              gACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACAt                                                                                                                                        <  2:1222504/150‑1 (MQ=255)

              gACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACAt                                                                                                                                        <  1:49910/150‑1 (MQ=255)

                 cTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCTTACGTCTGGCGcaaggaga                                                                                                                                     <  1:662673/150‑7 (MQ=255)

                  tgtgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAAc                                                                                                                                    >  1:244992/1‑150 (MQ=255)

                  tgtgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAAc                                                                                                                                    <  1:370618/150‑1 (MQ=255)

                  tgggATTGCGTCGCCTATAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCCTCTCTCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAAc                                                                                                                                    <  2:292646/147‑1 (MQ=255)

                    tgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAg                                                                                                                                  <  2:303298/150‑1 (MQ=255)

                       ttGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCtt                                                                                                                               <  2:1070135/150‑1 (MQ=255)

                       ttGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCtt                                                                                                                               >  2:632133/1‑150 (MQ=255)

                           gTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCCCATCAACAGCTTAAGc                                                                                                                           >  1:547068/1‑150 (MQ=255)

                            tCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCt                                                                                                                          >  1:603028/1‑150 (MQ=255)

                             cGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTg                                                                                                                         <  1:937814/150‑1 (MQ=255)

                              gCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGc                                                                                                                        >  1:812204/1‑150 (MQ=255)

                               cTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATATATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGAAAGGCGCACATAAACAGCTTAAGCTGcc                                                                                                                       >  2:336537/1‑150 (MQ=255)

                                   tAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGCCGTAAGGCGCACATCAACAGCTTAAGCTGCCCAcg                                                                                                                   >  1:1566194/1‑150 (MQ=255)

                                     aGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCAcgcg                                                                                                                 >  1:159619/1‑150 (MQ=255)

                                      gCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGa                                                                                                                <  2:1347847/150‑1 (MQ=255)

                                           tGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAgg                                                                                                           <  2:1273637/150‑1 (MQ=255)

                                           tGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAgg                                                                                                           <  2:960298/150‑1 (MQ=255)

                                             gCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTa                                                                                                         <  2:72925/150‑1 (MQ=255)

                                             gCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTa                                                                                                         >  2:1472697/1‑150 (MQ=255)

                                              cAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTAt                                                                                                        >  1:976343/1‑150 (MQ=255)

                                                 gcagcCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCtt                                                                                                     <  2:1120350/150‑1 (MQ=255)

                                                   agcCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTcc                                                                                                   >  1:1205475/1‑150 (MQ=255)

                                                   agcCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTcc                                                                                                   <  2:1720319/150‑1 (MQ=255)

                                                    gcCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAACAGGTAGAGACACCGCCGCCCGCTTTCCCCGGTGAGGAACTCTCCTGGCGCCAGACGTACGGCTCACATCAACAGCTTAAGCTGACCACGCGAAAAGGAATCTACCa                                                                                                  >  1:1780273/1‑150 (MQ=255)

                                                       cgcgACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCg                                                                                               >  1:55850/1‑150 (MQ=255)

                                                        gcgACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCgg                                                                                              >  2:958477/1‑150 (MQ=255)

                                                                catcatcaAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGt                                                                                      >  2:1426904/1‑150 (MQ=255)

                                                                    atcaAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgc                                                                                  >  2:1476823/1‑150 (MQ=255)

                                                                       aaGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGGGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGcgcg                                                                               <  1:1076373/150‑1 (MQ=255)

                                                                             agagGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCGCATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGGCGGCGCGAATGCg                                                                         <  1:259104/150‑1 (MQ=255)

                                                                               agGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                       <  1:1462673/150‑1 (MQ=255)

                                                                               agGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                       <  2:159619/150‑1 (MQ=255)

                                                                                ggCAAAATCATCTATCAGGTAGAGACACCTCCGCCCGCTTTCAAGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCATAAGCTGCCAACTCGATAAGGTATCTTCCACCGGCACACCGTCGGCGCGAATGCGaaa                                                                      >  1:583918/1‑150 (MQ=255)

                                                                                ggCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACCGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTat                                                                      >  1:1689809/1‑149 (MQ=255)

                                                                                  cAAAATCATCTAACAGGTAGAGACACCGCCGCAGGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCTCAACAACAGCTTAAGCTGCCCACGGGATAAGGTATCTTCCACCGGCGCACCGTAGGCGGGACAGCGAAAGt                                                                    >  2:903770/1‑150 (MQ=255)

                                                                                          tcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtg                                                            >  1:1560124/1‑150 (MQ=255)

                                                                                           cTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtgt                                                           >  1:1512755/1‑150 (MQ=255)

                                                                                           cTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtgt                                                           >  1:402550/1‑150 (MQ=255)

                                                                                            tATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtgtg                                                          <  1:1480828/150‑1 (MQ=255)

                                                                                             aTCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGAAAGTCTCCTTTGTGc                                                         >  2:259084/1‑150 (MQ=255)

                                                                                                aGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCggg                                                      <  1:10190/150‑1 (MQ=255)

                                                                                                aGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCggg                                                      <  1:980604/150‑1 (MQ=255)

                                                                                                 ggTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGc                                                     >  1:861126/1‑150 (MQ=255)

                                                                                                 ggTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGc                                                     >  1:176086/1‑150 (MQ=255)

                                                                                                   tAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCg                                                   >  2:1079964/1‑150 (MQ=255)

                                                                                                       gaCACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTgcg                                               >  1:1561769/1‑150 (MQ=255)

                                                                                                       gaCACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTgcg                                               >  1:1057410/1‑150 (MQ=255)

                                                                                                         cacCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACTGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGcgcg                                             >  1:717045/1‑150 (MQ=255)

                                                                                                         cacCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGcgcg                                             >  1:1563720/1‑150 (MQ=255)

                                                                                                         cacCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGAACATAAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCAACCGGCCCACCGTCGGCGCGAATGAGTAAGTCCGCTTTGTGCGGGCCGTCcgcg                                             >  1:1490569/1‑150 (MQ=255)

                                                                                                           ccgccgCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGgt                                           >  1:611025/1‑150 (MQ=255)

                                                                                                            cgccgcCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGgtg                                          <  2:128514/150‑1 (MQ=255)

                                                                                                            cgccgcCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGgtg                                          <  1:1438355/150‑1 (MQ=255)

                                                                                                            cgccgcCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGgtg                                          <  1:790672/150‑1 (MQ=255)

                                                                                                              ccgcccgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTa                                        <  1:632133/150‑1 (MQ=255)

                                                                                                                gcccgcTTGCACGGGTGACGAACCCTCCTTGCGCCAGCCGTACGGCGCACATCAACAGCTAAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAgg                                      <  1:303123/150‑1 (MQ=255)

                                                                                                                     cTTTCCCGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGc                                 <  1:1814470/150‑1 (MQ=255)

                                                                                                                      tttCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCt                                <  2:601284/150‑1 (MQ=255)

                                                                                                                       ttCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTg                               <  1:1131265/150‑1 (MQ=255)

                                                                                                                        tCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTgg                              <  2:1753216/150‑1 (MQ=255)

                                                                                                                        tCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTgg                              <  2:228171/150‑1 (MQ=255)

                                                                                                                           cGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGa                           >  1:10866/1‑150 (MQ=255)

                                                                                                                             ggTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATc                         <  2:257140/150‑1 (MQ=255)

                                                                                                                               tGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATcgc                       <  2:726966/150‑1 (MQ=255)

                                                                                                                                gcggAACTCTCCATGCGCCAGACGGACGGCGCACAACAACAGCATAAGCAGCCCCCGCGATAAGGAATCTTCCACCGGCGCACCGTCGGCGCGCATGCGTAAGTCCGCTTCGTGCGGGCCGGGCGCGGTGTAGGTTAGCTGGCGATcgcg                      <  2:861940/148‑1 (MQ=255)

                                                                                                                                gAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATagcg                      >  2:823565/1‑150 (MQ=255)

                                                                                                                                gAGGAACTCTCCTTGCGCCAGACGGAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGCTAAGGTATCTTCCACCGGCGCACCGTCGGCGCTAATGCGTAAGTCCGCGTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATcgcg                      <  2:878807/150‑1 (MQ=255)

                                                                                                                                 aGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGt                     >  2:435432/1‑150 (MQ=255)

                                                                                                                                  ggAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATTGCGtt                    >  2:347877/1‑150 (MQ=255)

                                                                                                                                  ggAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGtt                    <  2:923903/150‑1 (MQ=255)

                                                                                                                                   gAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTc                   <  1:787033/150‑1 (MQ=255)

                                                                                                                                     aCTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTCAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCaa                 <  1:291834/150‑1 (MQ=255)

                                                                                                                                       tctcCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCaaaa               >  1:576567/1‑150 (MQ=255)

                                                                                                                                       tctcCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCaaaa               >  1:1529851/1‑150 (MQ=255)

                                                                                                                                        ctcCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGCTCGCGGTCAAAAt              >  2:1394708/1‑150 (MQ=255)

                                                                                                                                          ccTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTa            >  2:702109/1‑150 (MQ=255)

                                                                                                                                          ccTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTa            <  1:1585627/150‑1 (MQ=255)

                                                                                                                                           cTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTAc           >  1:1530300/1‑150 (MQ=255)

                                                                                                                                             tGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATATGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGt         <  1:260264/150‑1 (MQ=255)

                                                                                                                                                  cAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAg    >  2:1765909/1‑150 (MQ=255)

                                                                                                                                                  cAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTAGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCTCGTTCAAAATTCCGTTCCAg    >  2:1765908/1‑150 (MQ=255)

                                                                                                                                                    gACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCa  >  2:1440801/1‑150 (MQ=255)

                                                                                                                                                    gACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCa  <  1:1712927/150‑1 (MQ=255)

ATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCA  >  NC_000913/3880227‑3880524

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