BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	403,334	C→T	T18T (ACC→ACT)	psiF →	PsiF family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	403,334	0	C	T	99.3%	543.9 / ‑4.4	135	T18T (ACC→ACT)	psiF	PsiF family protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (80/54); total (81/54)
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.52e-01

Only 100 of 139 total aligned reads displayed.

CGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGA  >  NC_000913/403186‑403480

cGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGt                                                                                                                                                  <  3:790255/151‑1 (MQ=255)

  cTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCg                                                                                                                                                >  3:350/1‑151 (MQ=255)

   tGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCgg                                                                                                                                               >  3:223084/1‑151 (MQ=255)

    gCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCgg                                                                                                                                               <  2:42144/150‑1 (MQ=255)

     ccGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGgcg                                                                                                                                             <  3:330172/151‑1 (MQ=255)

     ccGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGgcg                                                                                                                                             >  1:164488/1‑151 (MQ=255)

      cGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGgcgc                                                                                                                                            >  2:1099841/1‑151 (MQ=255)

      cGGGGGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGgcgc                                                                                                                                            >  1:35522/1‑151 (MQ=255)

           ggtTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCg                                                                                                                                       <  2:1257680/151‑1 (MQ=255)

               tttttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCgaga                                                                                                                                    <  3:1248668/150‑1 (MQ=255)

               tttttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCgag                                                                                                                                     >  4:32235/1‑149 (MQ=255)

                 tttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCgagaga                                                                                                                                  >  3:1012359/1‑150 (MQ=255)

                 tttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCgagaga                                                                                                                                  <  2:1272779/150‑1 (MQ=255)

                  ttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGaa                                                                                                                                 >  3:1407223/1‑150 (MQ=255)

                    gCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACtt                                                                                                                              <  1:186623/151‑1 (MQ=255)

                      tGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACttt                                                                                                                             >  1:150845/1‑150 (MQ=255)

                        ttAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTaa                                                                                                                           >  1:198348/1‑150 (MQ=255)

                         tAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAAc                                                                                                                          <  3:907203/150‑1 (MQ=255)

                          aGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAAccc                                                                                                                        >  2:609493/1‑151 (MQ=255)

                           gCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAAcccc                                                                                                                       >  3:994476/1‑151 (MQ=255)

                           gCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAAcccc                                                                                                                       <  2:600651/151‑1 (MQ=255)

                             aaCCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCcac                                                                                                                     >  2:771896/1‑151 (MQ=255)

                             aaCCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCcac                                                                                                                     >  4:433307/1‑151 (MQ=255)

                             aaCCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCcac                                                                                                                     <  4:198346/151‑1 (MQ=255)

                              aCCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCcac                                                                                                                     >  1:405175/1‑150 (MQ=255)

                               ccAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCcaca                                                                                                                    >  3:172899/1‑150 (MQ=255)

                               ccAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCcaca                                                                                                                    >  3:1030142/1‑150 (MQ=255)

                                     ttAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTGTAAACACTGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaac                                                                                                               >  2:755465/1‑149 (MQ=255)

                                      tAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaaca                                                                                                              <  3:546833/149‑1 (MQ=255)

                                      tAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGc                                                                                                            >  2:207417/1‑151 (MQ=255)

                                       aaTGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCg                                                                                                           >  3:646640/1‑151 (MQ=255)

                                       aaTGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCg                                                                                                           >  3:1196889/1‑151 (MQ=255)

                                       aaTGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGG‑TTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGt                                                                                                          >  3:918251/1‑151 (MQ=255)

                                          ggCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTa                                                                                                         <  3:999893/150‑1 (MQ=255)

                                            cAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATg                                                                                                       >  3:1390723/1‑150 (MQ=255)

                                             aGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAc                                                                                                     >  2:743431/1‑151 (MQ=255)

                                               aTCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcc                                                                                                    >  1:64005/1‑150 (MQ=255)

                                                 cACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctcc                                                                                                 >  4:165776/1‑151 (MQ=255)

                                                 cACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctc                                                                                                  <  3:70820/150‑1 (MQ=255)

                                                   cGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctcct                                                                                                <  2:1076387/150‑1 (MQ=255)

                                                   cGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTg                                                                                               <  2:834191/151‑1 (MQ=255)

                                                   cGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTg                                                                                               >  3:1165533/1‑151 (MQ=255)

                                                    gggCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGt                                                                                              >  2:931325/1‑151 (MQ=255)

                                                    gggCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGt                                                                                              >  4:114287/1‑151 (MQ=255)

                                                         cATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATca                                                                                         <  3:395495/151‑1 (MQ=255)

                                                         cATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATc                                                                                          >  2:1207853/1‑150 (MQ=255)

                                                             cGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATcagc                                                                                       >  2:121893/1‑149 (MQ=255)

                                                              gCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAgg                                                                                    <  3:677848/151‑1 (MQ=255)

                                                                 cATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCg                                                                                  <  3:1077426/150‑1 (MQ=255)

                                                                     gTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgc                                                                             <  3:1226863/151‑1 (MQ=255)

                                                                       tAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgcg                                                                            <  3:662991/150‑1 (MQ=255)

                                                                        aaaaCATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgcgc                                                                           <  2:978536/150‑1 (MQ=255)

                                                                         aaaCATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgcgc                                                                           >  2:57255/1‑149 (MQ=255)

                                                                            cATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCg                                                                      >  3:803120/1‑151 (MQ=255)

                                                                             aTGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGt                                                                     >  2:568752/1‑151 (MQ=255)

                                                                              tGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGtt                                                                    >  2:693659/1‑151 (MQ=255)

                                                                               gAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTg                                                                   >  4:445380/1‑151 (MQ=255)

                                                                                  gagGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGaaa                                                                <  3:511809/151‑1 (MQ=255)

                                                                                  gagGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGaaa                                                                <  4:251736/151‑1 (MQ=255)

                                                                                       aTGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAgggg                                                            >  3:739663/1‑150 (MQ=255)

                                                                                        tGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGAt                                                          >  3:1161498/1‑151 (MQ=255)

                                                                                          gTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGc                                                        <  1:218848/151‑1 (MQ=255)

                                                                                            gCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTc                                                      >  3:877130/1‑151 (MQ=255)

                                                                                              tATGAAAATAACATTACTGGTTAACTTGCTTTTCGGTATGGTTTTTTTAACCACTGTCGGCGCTGTCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGCAATCAGCAGGCGACGGCGCAGGCGTTGAACGGGGATGCTCg                                                     >  3:1406826/1‑150 (MQ=255)

                                                                                                 gAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAg                                                 <  1:395556/151‑1 (MQ=255)

                                                                                                  aaaaTAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTa                                                   >  2:37928/1‑148 (MQ=255)

                                                                                                   aaaTAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGAc                                               <  2:493415/151‑1 (MQ=255)

                                                                                                   aaaTAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGAc                                               <  3:121908/151‑1 (MQ=255)

                                                                                                    aaTAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGAcc                                              >  2:531936/1‑151 (MQ=255)

                                                                                                      tAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCt                                             >  3:1258450/1‑150 (MQ=255)

                                                                                                      tAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCt                                             >  3:706427/1‑150 (MQ=255)

                                                                                                       aaCATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTAc                                           >  2:429969/1‑151 (MQ=255)

                                                                                                        aCATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTAc                                           >  3:1248987/1‑150 (MQ=255)

                                                                                                        aCATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACa                                          <  2:223059/151‑1 (MQ=255)

                                                                                                         cATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACAt                                         >  3:196005/1‑151 (MQ=255)

                                                                                                           ttACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACAt                                         <  3:1032341/149‑1 (MQ=255)

                                                                                                           ttACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATg                                        >  1:201331/1‑150 (MQ=255)

                                                                                                                ggTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGa                                   <  2:1407152/150‑1 (MQ=255)

                                                                                                                ggTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGAt                                  >  2:829027/1‑151 (MQ=255)

                                                                                                                  ttACCTTGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTg                                >  2:732521/1‑151 (MQ=255)

                                                                                                                        tGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTga                           >  2:770954/1‑150 (MQ=255)

                                                                                                                        tGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTga                           <  2:350/150‑1 (MQ=255)

                                                                                                                        tGCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTg                            <  2:639128/149‑1 (MQ=255)

                                                                                                                         gCTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTga                           >  3:87231/1‑149 (MQ=255)

                                                                                                                          cTTTTCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTga                           >  2:646249/1‑148 (MQ=255)

                                                                                                                            tttCGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTgaaga                        <  2:542083/149‑1 (MQ=255)

                                                                                                                               cGGTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGc                   >  3:973227/1‑151 (MQ=255)

                                                                                                                                ggTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGc                   >  2:324822/1‑150 (MQ=255)

                                                                                                                                 gTCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCaa                 <  3:1088734/151‑1 (MQ=255)

                                                                                                                                  tCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCa                  >  2:1092380/1‑149 (MQ=255)

                                                                                                                                  tCTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCa                  >  3:409431/1‑149 (MQ=255)

                                                                                                                                   cTGGTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGt               >  1:66780/1‑151 (MQ=255)

                                                                                                                                     ggTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAg                <  3:612459/148‑1 (MQ=255)

                                                                                                                                     ggTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTc              <  3:531970/150‑1 (MQ=255)

                                                                                                                                      gTTTTTTTAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTg            >  2:278570/1‑151 (MQ=255)

                                                                                                                                          ttttAACCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGcc        >  2:799696/1‑151 (MQ=255)

                                                                                                                                              aaCCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTgg     <  1:14306/150‑1 (MQ=255)

                                                                                                                                              aaCCACTGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGc    <  1:433311/151‑1 (MQ=255)

                                                                                                                                               accactGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCg   <  3:519211/151‑1 (MQ=255)

                                                                                                                                                 cactGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGa  >  4:292285/1‑150 (MQ=255)

CGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGA  >  NC_000913/403186‑403480

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