BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,677,840	C→A	5.4%	G171W (GGG→TGG)	yhjG ←	AsmA family protein YhjG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,677,840	0	C	A	5.4%	544.4 / 10.2	159	G171W (GGG→TGG)	yhjG	AsmA family protein YhjG
Reads supporting (aligned to +/- strand): ref base C (74/76); new base A (6/3); total (80/79)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 162 total aligned reads displayed.

TCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCTTGA  >  NC_000913/3677693‑3677989

tCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCtt                                                                                                                                                     <  1:1539529/150‑1 (MQ=255)

   ggCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGCTCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGAt                                                                                                                                                  >  1:2001856/1‑150 (MQ=255)

   ggCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGAt                                                                                                                                                  >  2:116360/1‑150 (MQ=255)

   ggCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAAGGGCGGATCGCCATTGAt                                                                                                                                                  <  1:1161031/150‑1 (MQ=255)

    gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCATTAACTTGGCTGAACGGAAGCGGCTTGCCTAACGGAGCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATc                                                                                                                                                 <  2:1880978/150‑1 (MQ=255)

    gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAAGGGCGGATCGCCATTGAtg                                                                                                                                                 <  2:2104648/150‑2 (MQ=255)

     ccAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCg                                                                                                                                                >  1:1120677/1‑150 (MQ=255)

     ccAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCg                                                                                                                                                <  2:1026846/150‑1 (MQ=255)

     ccAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGAACAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCg                                                                                                                                                >  2:1986639/1‑150 (MQ=255)

      cAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGa                                                                                                                                               >  1:183850/1‑150 (MQ=255)

      cAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAAGGGCGGATCGCCATTGAtgac                                                                                                                                               <  2:2423426/150‑4 (MQ=255)

       aaaaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaa                                                                                                                                              >  1:1723889/1‑150 (MQ=255)

         aaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                                            <  1:1067216/150‑1 (MQ=255)

         aaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                                            <  2:388376/150‑1 (MQ=255)

         aaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                                            <  2:345304/150‑1 (MQ=255)

         aaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                                            >  1:1752721/1‑150 (MQ=255)

          aaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGCCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAg                                                                                                                                           >  2:1994692/1‑150 (MQ=255)

             ttAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGTACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAat                                                                                                                                        <  2:1285990/149‑1 (MQ=255)

             gTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGAACCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAat                                                                                                                                        >  1:1668004/1‑150 (MQ=255)

              tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                                       <  1:2398591/150‑1 (MQ=255)

               aaTCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAatat                                                                                                                                      >  1:1915190/1‑150 (MQ=255)

                aTCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAAtatt                                                                                                                                     <  1:296233/150‑1 (MQ=255)

                 tCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAAtatta                                                                                                                                    <  1:796910/150‑1 (MQ=255)

                      cACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAg                                                                                                                               >  2:306420/1‑150 (MQ=255)

                      cACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATACTATCCAg                                                                                                                               >  2:2092994/1‑150 (MQ=255)

                         cTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCg                                                                                                                            >  1:756940/1‑150 (MQ=255)

                         cTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCg                                                                                                                            >  1:2436514/1‑150 (MQ=255)

                          ttttCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGa                                                                                                                           <  2:1642144/150‑1 (MQ=255)

                             tCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAAc                                                                                                                        >  2:1142980/1‑150 (MQ=255)

                               tttATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGa                                                                                                                      >  2:2311273/1‑150 (MQ=255)

                                  aTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAACTCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCa                                                                                                                   <  2:1561678/150‑1 (MQ=255)

                                  aTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCa                                                                                                                   >  1:924361/1‑150 (MQ=255)

                                     cGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGc                                                                                                                >  2:159025/1‑150 (MQ=255)

                                      gCTTTACCTTTCGCTCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCTAAAAGAATATTATCCAGCCGAAACGACCATGcc                                                                                                               <  2:682868/150‑1 (MQ=255)

                                      gCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGcc                                                                                                               <  2:1395063/150‑1 (MQ=255)

                                              tttCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCtt                                                                                                       <  2:923501/150‑1 (MQ=255)

                                              tttCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCtt                                                                                                       >  1:369930/1‑150 (MQ=255)

                                                 cGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGc                                                                                                    >  2:2457595/1‑150 (MQ=255)

                                                 cGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGc                                                                                                    <  2:266688/150‑1 (MQ=255)

                                                 cGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGc                                                                                                    <  1:2403478/150‑1 (MQ=255)

                                                  gATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCGTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGct                                                                                                   >  2:471097/1‑149 (MQ=255)

                                                      cAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcg                                                                                               <  1:1681534/150‑1 (MQ=255)

                                                       aGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcgc                                                                                              >  2:1031622/1‑150 (MQ=255)

                                                              tCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTg                                                                                       <  2:1367790/150‑1 (MQ=255)

                                                              tCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTg                                                                                       >  2:1629446/1‑150 (MQ=255)

                                                               cGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGt                                                                                      <  2:1752721/150‑1 (MQ=255)

                                                                gCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGtt                                                                                     <  1:1233279/150‑1 (MQ=255)

                                                                  tGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTat                                                                                   >  1:1191194/1‑150 (MQ=255)

                                                                   gAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCACATCCGCTTTGCTAACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTatc                                                                                  >  2:2168395/1‑150 (MQ=255)

                                                                     aCGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTatcat                                                                                >  1:484541/1‑150 (MQ=255)

                                                                       ggCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCATACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCg                                                                              >  2:159727/1‑150 (MQ=255)

                                                                          aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTg                                                                           <  1:2270272/150‑1 (MQ=255)

                                                                             ggCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCa                                                                        >  2:2033535/1‑150 (MQ=255)

                                                                               cTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAg                                                                      <  2:753818/150‑1 (MQ=255)

                                                                               cTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAg                                                                      <  2:1353441/150‑1 (MQ=255)

                                                                                  gCCTAAGGGAACAACAAAAATCTCCAGATCCGCATTGCTTACATTGTCATCAATGGCGATCCGCCCATGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTCTTAACATCGATGCCAAGAtt                                                                   >  1:1280525/1‑150 (MQ=255)

                                                                                    cTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTaa                                                                 >  1:1687008/1‑150 (MQ=255)

                                                                                     tAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTaaa                                                                <  1:1172600/150‑1 (MQ=255)

                                                                                      aaGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAAc                                                               <  1:1994692/150‑1 (MQ=255)

                                                                                         ggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTc                                                            >  2:93947/1‑150 (MQ=255)

                                                                                            tCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAg                                                         <  2:645510/150‑1 (MQ=255)

                                                                                               acaaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTa                                                      >  1:2057089/1‑150 (MQ=255)

                                                                                                caaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGATATCATAGTTGGCAAGATTAAACGTCCAGatat                                                     >  1:2228709/1‑150 (MQ=255)

                                                                                                  aaaaTCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTg                                                   <  2:369930/150‑1 (MQ=255)

                                                                                                        tcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCttt                                             >  2:2252425/1‑150 (MQ=255)

                                                                                                             aTCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCaga                                        <  1:306420/150‑1 (MQ=255)

                                                                                                                     tGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGaa                                <  1:1500715/150‑1 (MQ=255)

                                                                                                                       cTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATc                              >  2:1203305/1‑150 (MQ=255)

                                                                                                                         tACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAg                            >  2:2337129/1‑150 (MQ=255)

                                                                                                                         tACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAg                            >  2:327802/1‑150 (MQ=255)

                                                                                                                           cttttcgatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAACCGACCATGCCGCCGGCTTTGGATTAGCGTCTTTGTTATCAGCGTGGGCATGCTTAAACGGCCAGTTATTGTTCTTTTCGGCGAGGCGAATCAGgc                          >  2:849176/8‑150 (MQ=255)

                                                                                                                            tttGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGgcg                         <  1:2199028/150‑1 (MQ=255)

                                                                                                                              tGTCATCTATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGgcgcg                       <  2:519059/150‑1 (MQ=255)

                                                                                                                               tcgatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGgcgcgc                      >  2:1910587/4‑150 (MQ=255)

                                                                                                                                 gatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt                    >  1:189612/2‑150 (MQ=255)

                                                                                                                                 gatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt                    >  1:552835/2‑150 (MQ=255)

                                                                                                                                 catcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt                    >  1:2339728/1‑150 (MQ=255)

                                                                                                                                 catcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt                    >  1:632929/1‑150 (MQ=255)

                                                                                                                                  atcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTc                   >  2:1812659/1‑150 (MQ=255)

                                                                                                                                   tcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCg                  <  1:327430/150‑1 (MQ=255)

                                                                                                                                    caATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCgg                 >  1:1998539/1‑150 (MQ=255)

                                                                                                                                    caATGGCGATCCGCCCTTGATCG‑AAAGAATATTATCCAGCCGAAACGACCATGCCGCCGGCTTTGCATTAGCGTCGTTGTTTTCATCCTTGGCAAGATTAAACGTCCAATTTTTGTTCTTTTCAGAGAGGCCAATCAGGCGCGCGTCggg                >  2:2190798/1‑150 (MQ=255)

                                                                                                                                     aaTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTc                                           >  1:1049137/1‑123 (MQ=255)

                                                                                                                                     aaTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCggg                <  2:183850/150‑1 (MQ=255)

                                                                                                                                     aaTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCggg                <  2:1324070/150‑1 (MQ=255)

                                                                                                                                      aTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGc               <  2:1823757/150‑1 (MQ=255)

                                                                                                                                      aTGGCGATCCGCCCGTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGGTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGc               <  2:1315078/150‑1 (MQ=255)

                                                                                                                                       tGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGGCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCt              <  1:2033535/150‑1 (MQ=255)

                                                                                                                                         gCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCttt            >  1:1793798/1‑150 (MQ=255)

                                                                                                                                           gATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTc          >  1:1111007/1‑150 (MQ=255)

                                                                                                                                           gATCCGCCCTTGATCGAAAAGAATATAATCCAGCCGAAACCACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTc          >  2:904941/1‑150 (MQ=255)

                                                                                                                                               cGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGc      <  1:2311273/150‑1 (MQ=255)

                                                                                                                                               cGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGc      >  1:769205/1‑150 (MQ=255)

                                                                                                                                               cGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGc      <  1:159025/150‑1 (MQ=255)

                                                                                                                                                gCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCt     <  1:342880/150‑1 (MQ=255)

                                                                                                                                                gCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCt     <  1:159727/150‑1 (MQ=255)

                                                                                                                                                gCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCt     <  1:688837/150‑1 (MQ=255)

                                                                                                                                                gCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCt     <  1:720159/150‑1 (MQ=255)

                                                                                                                                                gCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCt     >  1:1648048/1‑150 (MQ=255)

                                                                                                                                                   cTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCTTGa  <  2:1067910/150‑1 (MQ=255)

TCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGCTTGA  >  NC_000913/3677693‑3677989

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: `‑`