BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,881,075	A→G	6.0%	G49G (GGT→GGC)	recF ←	DNA repair protein RecF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,881,075	0	A	G	6.0%	453.0 / 13.3	154	G49G (GGT→GGC)	recF	DNA repair protein RecF
Reads supporting (aligned to +/- strand): ref base A (66/77); new base G (7/2); total (75/79)
Fisher's exact test for biased strand distribution p-value = 8.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 155 total aligned reads displayed.

CGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTAC  >  NC_000913/3880927‑3881224

ggTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGATCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGAc                                                                                                                                                      <  2:112473/149‑1 (MQ=255)

  tCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCa                                                                                                                                                    <  1:235796/150‑1 (MQ=255)

    gATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATg                                                                                                                                                  <  1:1258333/150‑1 (MQ=255)

    gATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATg                                                                                                                                                  <  2:1990513/150‑1 (MQ=255)

     aTGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGCCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATgg                                                                                                                                                 >  2:36917/1‑150 (MQ=255)

      tGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGc                                                                                                                                                >  2:317365/1‑150 (MQ=255)

      tGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGc                                                                                                                                                <  1:1692618/150‑1 (MQ=255)

      tGCGGACTTTGCTGTCGCCCTGTTTCTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGc                                                                                                                                                >  1:2076561/1‑150 (MQ=255)

         ggACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGa                                                                                                                                             <  1:1039244/150‑1 (MQ=255)

         ggACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGa                                                                                                                                             >  2:1946943/1‑150 (MQ=255)

            cTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCg                                                                                                                                          <  1:2485562/150‑1 (MQ=255)

            cTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCg                                                                                                                                          <  1:2079958/150‑1 (MQ=255)

             tttGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGt                                                                                                                                         <  1:704144/150‑1 (MQ=255)

              ttGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGta                                                                                                                                        >  2:1942184/1‑150 (MQ=255)

              ttGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGta                                                                                                                                        >  1:59898/1‑150 (MQ=255)

               tGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGtat                                                                                                                                       <  2:807695/150‑1 (MQ=255)

                gCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGtata                                                                                                                                      <  1:365516/150‑1 (MQ=255)

                 cTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCTCGGCCATGGTCGGGCGtttcg                                                                                                                                     <  1:1112538/150‑5 (MQ=255)

                 cTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCTCGGCCATGGTCGGGCGtttcg                                                                                                                                     <  1:1121035/150‑5 (MQ=255)

                    tCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATg                                                                                                                                  <  1:2204978/150‑1 (MQ=255)

                      gCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGc                                                                                                                                <  1:1748522/150‑1 (MQ=255)

                          tGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTcc                                                                                                                            >  2:2092904/1‑150 (MQ=255)

                           gTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGATGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCa                                                                                                                           <  1:1864303/150‑1 (MQ=255)

                             ttGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGc                                                                                                                         >  1:1424036/1‑150 (MQ=255)

                                 cTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGta                                                                                                                                        <  1:2290071/131‑1 (MQ=255)

                                    tGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGCCCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTgg                                                                                                                  >  2:1726767/1‑150 (MQ=255)

                                     ggTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGt                                                                                                                 <  2:1021339/150‑1 (MQ=255)

                                      gTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGtt                                                                                                                >  1:50185/1‑150 (MQ=255)

                                      gTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGtt                                                                                                                >  1:1848456/1‑150 (MQ=255)

                                       ttAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGttt                                                                                                               >  2:1729619/1‑150 (MQ=255)

                                         aaGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTg                                                                                                             <  2:2044696/150‑1 (MQ=255)

                                          aGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGc                                                                                                            >  1:568244/1‑150 (MQ=255)

                                           gCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGcc                                                                                                           <  1:380711/150‑1 (MQ=255)

                                           gCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGACAGTGGAGAACAAACGCCTCCTGCCCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGcc                                                                                                           <  2:859117/150‑1 (MQ=255)

                                                tCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAAAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGc                                                                                                      >  2:1590207/1‑150 (MQ=255)

                                                     gTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTg                                                                                                 <  1:1802641/150‑1 (MQ=255)

                                                       ctcGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGc                                                                                               >  2:769143/1‑150 (MQ=255)

                                                       ctcGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGc                                                                                               <  1:999015/150‑1 (MQ=255)

                                                          gcgcTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCAcc                                                                                            <  1:2470468/150‑1 (MQ=255)

                                                           cgcTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCt                                                                                           >  2:2333861/1‑150 (MQ=255)

                                                             ctctTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAc                                                                                         >  2:2163489/1‑150 (MQ=255)

                                                              tctTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAcc                                                                                        <  2:1351311/150‑1 (MQ=255)

                                                              tctTCGCCCTGTAATCGCCCGTGGACAACAAACGCCTCCTGCTCATGGTGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAcc                                                                                        <  1:1651961/150‑1 (MQ=255)

                                                                 tCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGACAGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGCTTTTGCCACTGCCGTTGGCACCTACCAGa                                                                                     >  1:759680/1‑150 (MQ=255)

                                                                    cccTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAg                                                                                  <  1:2042864/150‑1 (MQ=255)

                                                                     ccTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGt                                                                                 <  2:324879/150‑1 (MQ=255)

                                                                       tGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTa                                                                               >  1:2409505/1‑150 (MQ=255)

                                                                       tGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTa                                                                               <  2:2358644/150‑1 (MQ=255)

                                                                       tGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTa                                                                               <  2:636418/150‑1 (MQ=255)

                                                                         tAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCGGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTaaa                                                                             <  2:2076561/150‑1 (MQ=255)

                                                                         tAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTaaa                                                                             <  2:397070/150‑1 (MQ=255)

                                                                                ccGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGga                                                                      <  2:1014042/150‑1 (MQ=255)

                                                                                 cGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGgag                                                                     <  2:50185/150‑1 (MQ=255)

                                                                                    ggAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATa                                                                  >  1:1956741/1‑150 (MQ=255)

                                                                                      agaACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAg                                                                >  1:5941/1‑150 (MQ=255)

                                                                                      agaACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAg                                                                >  1:603485/1‑150 (MQ=255)

                                                                                       gaACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAgg                                                               >  1:2235218/1‑150 (MQ=255)

                                                                                       gaACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAgg                                                               >  1:781153/1‑150 (MQ=255)

                                                                                            aaCGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCgaga                                                          >  1:1846392/1‑150 (MQ=255)

                                                                                            aaCGCCTCCTGCTCATGCCGAAAGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCCACCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAACAGAAAGTTAAAGCCGGGAGATAAGTCgaga                                                          >  2:1898052/1‑150 (MQ=255)

                                                                                             aCGCCTCATGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGAt                                                         <  1:2387708/150‑1 (MQ=255)

                                                                                                 ctcctGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCcgc                                                     <  1:1942184/150‑1 (MQ=255)

                                                                                                  tcctGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCcgcg                                                    <  1:1925885/150‑1 (MQ=255)

                                                                                                     ggCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCGGGGTTTGCCACTGCGGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGCCGAGATCCGCGGtt                                                 <  2:1852772/149‑1 (MQ=255)

                                                                                                       cTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTc                                               >  2:2123402/1‑150 (MQ=255)

                                                                                                         cATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCaa                                             >  2:2487348/1‑150 (MQ=255)

                                                                                                         cATGGCGAATGACGCGACCAATCTGCAAACTCCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCTTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCaa                                             >  2:620792/1‑150 (MQ=255)

                                                                                                              cGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTg                                        <  2:1765502/150‑1 (MQ=255)

                                                                                                                  tGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGa                                    <  2:1424036/150‑1 (MQ=255)

                                                                                                                       cgACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATcg                               <  2:568244/150‑1 (MQ=255)

                                                                                                                          ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                            <  2:1812588/150‑1 (MQ=255)

                                                                                                                          ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                            <  2:1894613/150‑1 (MQ=255)

                                                                                                                          ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                            <  2:1779968/150‑1 (MQ=255)

                                                                                                                          ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                            <  1:1533237/150‑1 (MQ=255)

                                                                                                                           cAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGa                           <  2:1689381/150‑1 (MQ=255)

                                                                                                                              tCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATca                        <  2:2485528/150‑1 (MQ=255)

                                                                                                                              tCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATca                        <  2:1036319/150‑1 (MQ=255)

                                                                                                                               cTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaa                       <  1:1876568/150‑1 (MQ=255)

                                                                                                                               cTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaa                       >  1:782740/1‑150 (MQ=255)

                                                                                                                               cTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAACGCGAGATCCGCGGGTTCAATGTTGCGGAAATCGCGGATcaa                       <  2:2155602/150‑1 (MQ=255)

                                                                                                                                tGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaac                      >  1:1045021/1‑150 (MQ=255)

                                                                                                                                   aaaCTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAg                   >  1:37008/1‑150 (MQ=255)

                                                                                                                                      catctataCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCgg                >  1:1998047/8‑150 (MQ=255)

                                                                                                                                      cTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCgg                >  1:1196638/1‑150 (MQ=255)

                                                                                                                                       tGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCggg               <  1:325197/150‑1 (MQ=255)

                                                                                                                                        gCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGt              >  1:242842/1‑150 (MQ=255)

                                                                                                                                         cGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTg             <  2:1050101/150‑1 (MQ=255)

                                                                                                                                         cGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTg             <  2:759680/150‑1 (MQ=255)

                                                                                                                                          tataCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGa            >  1:1413221/4‑150 (MQ=255)

                                                                                                                                          tataCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGAATTTGGCACATACAAGAAAGTTAAAGCCGGGAGATAATACGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGa            >  2:1831554/4‑150 (MQ=255)

                                                                                                                                           aaaCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAg           >  2:804139/1‑150 (MQ=255)

                                                                                                                                           aaaCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAg           >  1:332582/1‑150 (MQ=255)

                                                                                                                                             aCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAggg         >  1:948566/1‑150 (MQ=255)

                                                                                                                                             aCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAggg         >  1:290971/1‑150 (MQ=255)

                                                                                                                                             aCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAggg         >  1:1146024/1‑150 (MQ=255)

                                                                                                                                               gCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGTCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGAc       <  2:1655995/150‑1 (MQ=255)

                                                                                                                                               gCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGAc       <  2:1262344/150‑1 (MQ=255)

                                                                                                                                                cccGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACa      >  1:1359791/1‑150 (MQ=255)

                                                                                                                                                cccGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACa      >  2:153541/1‑150 (MQ=255)

                                                                                                                                                   gACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGAAACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGAGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTa   >  2:57607/1‑150 (MQ=255)

CGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTAC  >  NC_000913/3880927‑3881224

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