BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,557,462	C→A	5.6%	E206D (GAG→GAT)	rtcB ←	RNA‑splicing ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,557,462	0	C	A	5.6%	529.5 / 15.3	144	E206D (GAG→GAT)	rtcB	RNA‑splicing ligase
Reads supporting (aligned to +/- strand): ref base C (64/72); new base A (4/4); total (68/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 150 total aligned reads displayed.

GTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTCCCAGGTGTTTATA  >  NC_000913/3557313‑3557611

gTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCAtc                                                                                                                                                       >  1:2317379/1‑150 (MQ=255)

   aCCTCGTTTTCCATCATCGCATCGCGGTTTAGGCTGGCATAAAGCTGCGCCCTGGCGACGGCTTTCAGGTAATCATCTAAGTATTCCGTACCTTCCATTTAGTACGCCAGATCACGCGACGGCATCGTCTCAAGCGTTTCCTGCTTctct                                                                                                                                                    <  2:1384502/150‑1 (MQ=255)

    ccACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCtt                                                                                                                                                   >  2:1829896/1‑150 (MQ=255)

    ccACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGAGTTTCCTGCATCTCtt                                                                                                                                                   >  1:1455168/1‑150 (MQ=255)

     cACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCttt                                                                                                                                                  >  2:2051834/1‑150 (MQ=255)

       cGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCttttt                                                                                                                                                <  2:336255/150‑1 (MQ=255)

         ttttCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTtgt                                                                                                                                              <  1:819633/150‑1 (MQ=255)

         ttttCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTtgt                                                                                                                                              <  1:961406/150‑1 (MQ=255)

           gtCCATCATCGCATCGCGGTGACGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGc                                                                                                                                            <  2:1394816/149‑1 (MQ=255)

              catcatCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAg                                                                                                                                         >  1:1416800/1‑150 (MQ=255)

                  atcgcatcgcGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCg                                                                                                                                     >  2:1738896/1‑150 (MQ=255)

                  atcgcatcgcGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCg                                                                                                                                     >  2:1738985/1‑150 (MQ=255)

                   tcgcatcgcGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGa                                                                                                                                    >  1:1022972/1‑150 (MQ=255)

                   tcgcatcgcGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGa                                                                                                                                    >  1:1002629/1‑150 (MQ=255)

                   tcgcatcgcGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGa                                                                                                                                    >  1:2216460/1‑150 (MQ=255)

                          gcgGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAaagta                                                                                                                             >  2:2032906/1‑150 (MQ=255)

                          gcgGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAaagta                                                                                                                             >  2:338369/1‑150 (MQ=255)

                          gcgGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCAAGATCGATAaagta                                                                                                                             >  2:2344438/1‑150 (MQ=255)

                          gcgGTTAAGGCTGGCAAAAAGCTCCGCCCAGGCCACGGATTTCAGGTAATCATCAAAGTATTCCGTACCATACATAAAGTACGCCAGATCACTCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTAATTGTGCCAGATCGATAaagta                                                                                                                             >  2:805238/1‑150 (MQ=255)

                           cgGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAaagtaa                                                                                                                            <  2:845498/150‑1 (MQ=255)

                           cgGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAaagtaa                                                                                                                            <  1:1349358/150‑1 (MQ=255)

                                aaGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTccc                                                                                                                       >  1:2098158/1‑150 (MQ=255)

                                aaGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTccc                                                                                                                       >  1:716640/1‑150 (MQ=255)

                                aaGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTccc                                                                                                                       >  1:2098108/1‑150 (MQ=255)

                                aaGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTccc                                                                                                                       >  1:1524156/1‑150 (MQ=255)

                                  ggCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGa                                                                                                                     <  2:1770210/150‑1 (MQ=255)

                                  ggCTGGCAAAAAGCTGCGCCCAGGCCAAGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGa                                                                                                                     <  2:1770220/150‑1 (MQ=255)

                                      ggCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGc                                                                                                                 <  2:217927/150‑1 (MQ=255)

                                        cAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGt                                                                                                               >  1:1310710/1‑150 (MQ=255)

                                         aaaaaGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGtt                                                                                                              >  2:2158669/1‑150 (MQ=255)

                                         aaaaaGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGtt                                                                                                              >  2:2158622/1‑150 (MQ=255)

                                         aaaaaGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGtt                                                                                                              >  2:1507654/1‑150 (MQ=255)

                                          aaaaGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACTCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTa                                                                                                             <  1:2174949/150‑1 (MQ=255)

                                           aaaGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTAc                                                                                                            >  2:248896/1‑150 (MQ=255)

                                             aGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCa                                                                                                          <  1:2073416/150‑1 (MQ=255)

                                              gCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCaa                                                                                                         <  1:2020740/150‑1 (MQ=255)

                                               cTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAAt                                                                                                        >  2:571000/1‑150 (MQ=255)

                                                 gcgcCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTc                                                                                                      >  2:1681603/1‑150 (MQ=255)

                                                 gcgcCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTc                                                                                                      <  2:1442425/150‑1 (MQ=255)

                                                    cccAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCcgc                                                                                                   >  1:1377543/1‑150 (MQ=255)

                                                    cccAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCcgc                                                                                                   <  1:2372528/150‑1 (MQ=255)

                                                       aGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTg                                                                                                >  2:74994/1‑150 (MQ=255)

                                                            aCGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCg                                                                                           <  2:910289/150‑1 (MQ=255)

                                                            aCGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCg                                                                                           <  1:1106896/150‑1 (MQ=255)

                                                            aCGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCg                                                                                           <  1:391131/150‑1 (MQ=255)

                                                               gCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAg                                                                                        >  1:1786810/1‑150 (MQ=255)

                                                                        tAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATa                                                                               >  2:822784/1‑150 (MQ=255)

                                                                         aaTCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATaa                                                                              <  1:64316/150‑1 (MQ=255)

                                                                         aaTCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATaa                                                                              <  1:1430298/150‑1 (MQ=255)

                                                                           tcatcaAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATc                                                                            <  1:474872/150‑1 (MQ=255)

                                                                           tcatcaAACTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATc                                                                            <  2:2303920/150‑1 (MQ=255)

                                                                              tcaAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCcac                                                                         <  1:1829896/150‑1 (MQ=255)

                                                                              tcaAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCcac                                                                         >  1:999443/1‑150 (MQ=255)

                                                                                aaaGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCacac                                                                       >  1:2082572/1‑150 (MQ=255)

                                                                                aaaGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCacac                                                                       >  2:735953/1‑150 (MQ=255)

                                                                                aaaGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCacac                                                                       >  2:1128850/1‑150 (MQ=255)

                                                                                aaaGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCacac                                                                       >  1:1700020/1‑150 (MQ=255)

                                                                                 aaGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATATCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACAcc                                                                      >  1:1781734/1‑150 (MQ=255)

                                                                                 aaGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATATCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACAcc                                                                      >  1:1932971/1‑150 (MQ=255)

                                                                                    tATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTgg                                                                   >  2:1215363/1‑150 (MQ=255)

                                                                                    tATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTgg                                                                   >  1:832421/1‑150 (MQ=255)

                                                                                    tATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTgg                                                                   >  1:801230/1‑150 (MQ=255)

                                                                                    tATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTgg                                                                   >  1:650266/1‑150 (MQ=255)

                                                                                       tCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTcc                                                                <  1:939000/150‑1 (MQ=255)

                                                                                       tCCGTACCTTCCATAAAATACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTcc                                                                <  1:1145954/150‑1 (MQ=255)

                                                                                          gTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACTTCTCAAGCGTTACCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCAGGTCCGAc                                                             >  2:1541247/1‑150 (MQ=255)

                                                                                          gTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTTCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGAc                                                             <  1:1955624/150‑1 (MQ=255)

                                                                                          gTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGACCCGGAGTGCAGCATAATCCACACCTGGTCCGAc                                                             >  2:470982/1‑150 (MQ=255)

                                                                                          gTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGAc                                                             >  2:78257/1‑150 (MQ=255)

                                                                                              cTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACtcag                                                         >  2:1518894/1‑149 (MQ=255)

                                                                                                 ccATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCaa                                                      <  1:252363/150‑1 (MQ=255)

                                                                                                  cATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAg                                                     <  1:403291/150‑1 (MQ=255)

                                                                                                  cATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAg                                                     <  1:569075/150‑1 (MQ=255)

                                                                                                   aTAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAgg                                                    <  1:1657292/150‑1 (MQ=255)

                                                                                                   aTAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAgg                                                    <  1:1868182/150‑1 (MQ=255)

                                                                                                       aGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGa                                                <  1:1681603/150‑1 (MQ=255)

                                                                                                         tACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGAtt                                              <  1:985850/150‑1 (MQ=255)

                                                                                                                 aTCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATaaa                                      <  1:1738896/150‑1 (MQ=255)

                                                                                                                 aTCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATaaa                                      <  1:1738985/150‑1 (MQ=255)

                                                                                                                    aCGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTg                                   <  1:564279/150‑1 (MQ=255)

                                                                                                                    aCGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTg                                   <  1:2051834/150‑1 (MQ=255)

                                                                                                                    aCGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCAACTCATCAAGGCAGATTTCAATAAAGTg                                   >  1:574744/1‑150 (MQ=255)

                                                                                                                    aCGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCAACTCATCAAGGCAGATTTCAATAAAGTg                                   >  1:649553/1‑150 (MQ=255)

                                                                                                                      gcgACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCTATAAAGTGGt                                 >  2:2079178/1‑150 (MQ=255)

                                                                                                                         aCGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTAc                              <  1:551378/150‑1 (MQ=255)

                                                                                                                            gCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCgg                           <  2:1786810/150‑1 (MQ=255)

                                                                                                                                  tctcAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTACGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCa                     <  1:444352/150‑1 (MQ=255)

                                                                                                                                  tctcAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCa                     <  1:822784/150‑1 (MQ=255)

                                                                                                                                    tcAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTCCCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATACAGTGGTTACCGGTTCCCCGc                   <  2:2215061/150‑1 (MQ=255)

                                                                                                                                    tcAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGc                   <  2:1966047/150‑1 (MQ=255)

                                                                                                                                    tcAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGc                   <  2:1809133/150‑1 (MQ=255)

                                                                                                                                    tcAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGc                   <  2:825155/150‑1 (MQ=255)

                                                                                                                                       aGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGtt                <  2:1310710/150‑1 (MQ=255)

                                                                                                                                        gCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTc               >  2:1501269/1‑150 (MQ=255)

                                                                                                                                        gCGTTTCCTGCATATCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTc               <  2:1932971/150‑1 (MQ=255)

                                                                                                                                        gCGTTTCCTGCATATCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTc               <  2:1482486/150‑1 (MQ=255)

                                                                                                                                                tGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTCCCAGgtgt       <  2:1828056/150‑1 (MQ=255)

                                                                                                                                                   aTCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTCCCAGGTGTTta    <  1:248896/150‑1 (MQ=255)

                                                                                                                                                    tctcTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTCCCAGGTGTTtat   >  1:1450680/1‑150 (MQ=255)

                                                                                                                                                     ctctTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTCCCAGGTGTTtata  >  1:1322531/1‑150 (MQ=255)

GTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCGACGGCAACGTCTCAAGCGTTTCCTGCATCTCTTTTTGTGCCAGATCGATAAAGTAAGTCCCGATGGCGTTACCAATTCCGCGTGAACCGGAGTGCAGCATAATCCACACCTGGTCCGACTCATCAAGGCAGATTTCAATAAAGTGGTTACCGGTTCCCAGCGTTCCCAGGTGTTTATA  >  NC_000913/3557313‑3557611

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