BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,694	A→G	9.5%	T92A (ACG→GCG)	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,180,694	0	A	G	9.5%	365.0 / 24.7	139	T92A (ACG→GCG)	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base A (56/67); new base G (5/8); total (64/75)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-02

Only 100 of 139 total aligned reads displayed.

GCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAA  >  NC_000913/3180547‑3180839

gCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACg                                                                                                                                                 <  1:1324488/150‑1 (MQ=255)

 cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGCACTGCACCAGATGTGCCTGAAAGTGGGGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATCCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGc                                                                                                                                                <  1:614351/150‑1 (MQ=255)

 cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGc                                                                                                                                                <  2:238359/150‑1 (MQ=255)

 cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGc                                                                                                                                                >  1:154675/1‑150 (MQ=255)

 cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTG‑CTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCa                                                                                                                                               >  1:835542/1‑150 (MQ=255)

    aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGTCAGCCATGACTGGCGCAca                                                                                                                                             <  2:1272637/150‑2 (MQ=255)

    aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGAAcc                                                                                                                                             >  1:1040050/1‑150 (MQ=255)

       ttGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGc                                                                                                                                          <  1:813721/150‑1 (MQ=255)

         gAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGCTGACCAAATTCCGCATTCCAAAACACACATGGAGTTTTGTGCGCCAGTCCTGGCTGACGCACCAGCCa                                                                                                                                        >  1:232476/1‑150 (MQ=255)

          aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAt                                                                                                                                       >  2:1708303/1‑150 (MQ=255)

          aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAt                                                                                                                                       >  1:458060/1‑150 (MQ=255)

          aaaaGCTGGAAGAAGGCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAt                                                                                                                                       <  1:723966/150‑1 (MQ=255)

           aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtg                                                                                                                                      <  2:1583995/150‑2 (MQ=255)

           aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATc                                                                                                                                      <  2:1198700/150‑1 (MQ=255)

             aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtgac                                                                                                                                    <  2:12029/150‑4 (MQ=255)

             aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGc                                                                                                                                    <  1:1688785/150‑1 (MQ=255)

             aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGc                                                                                                                                    <  1:1695381/150‑1 (MQ=255)

               cTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCtt                                                                                                                                  <  1:342144/150‑1 (MQ=255)

                tGGAAGAAGTCACCGCCGAACTGCACCAGATTTGCCTGAACGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTACGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCttt                                                                                                                                 <  1:941836/150‑1 (MQ=255)

                tGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCttt                                                                                                                                 <  2:1445972/150‑1 (MQ=255)

                 ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTTATGGCTGACGCACCAGCCATCGCTTTa                                                                                                                                >  1:850586/1‑150 (MQ=255)

                 ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTa                                                                                                                                >  1:683085/1‑150 (MQ=255)

                 ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGCAGTTTTGTGCGCCAGTCATGGATTACGCACCAGCCATCGCTTTa                                                                                                                                >  1:1125780/1‑150 (MQ=255)

                   aagaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                                                              <  2:1131972/150‑1 (MQ=255)

                   aagaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                                                              <  2:565879/150‑1 (MQ=255)

                   aagaagTCACCGCCGAACACCACCCGATGAGCCTTAAAGTGGTGGAAAAAGTGATCGCCAGCGACGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGAGTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                                                              <  1:864579/150‑1 (MQ=255)

                     gaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGCGGAAAAAGTGATCGCAAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcg                                                                                                                            <  2:1600789/150‑1 (MQ=255)

                     gaagGCACCGCCGAAATGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTTATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcg                                                                                                                            <  2:833072/150‑1 (MQ=255)

                      aagGCACCGCCGAACTGCACCTGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAACAACTCCTGGAGTTATGTGCGCCAGTCATGGCTGTCTCACCAGACATCGCTTTATTcgc                                                                                                                           >  2:887440/1‑150 (MQ=255)

                       agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                                                          <  1:793040/150‑1 (MQ=255)

                          cACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCt                                                                                                                       >  2:1151097/1‑150 (MQ=255)

                          cACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCt                                                                                                                       <  2:402851/150‑1 (MQ=255)

                            ccgccgAACTGCACCATATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAACCATCGCTTTATTCGCGTCTTg                                                                                                                     <  2:596239/150‑1 (MQ=255)

                            ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                     <  2:31202/150‑1 (MQ=255)

                            ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                     <  2:411449/150‑1 (MQ=255)

                            ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                     <  2:325101/150‑1 (MQ=255)

                             cgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGTTCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGa                                                                                                                    >  2:501654/1‑150 (MQ=255)

                             cgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGa                                                                                                                    >  2:1651773/1‑150 (MQ=255)

                              gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGACTGACGCACCAGCCATCGCtt                                                                                                                                  >  2:1698737/1‑135 (MQ=255)

                              gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGACTGACGCACCAGCCATCGCtt                                                                                                                                  <  1:1698737/135‑1 (MQ=255)

                               ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                  >  1:1010320/1‑150 (MQ=255)

                               ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                  <  1:679400/150‑1 (MQ=255)

                               ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                  >  2:1656406/1‑150 (MQ=255)

                                 gAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTg                                                                                                                <  1:901200/150‑1 (MQ=255)

                                    cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                             <  2:976714/150‑1 (MQ=255)

                                    cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                             <  1:179522/150‑1 (MQ=255)

                                     tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGt                                                                                                            >  2:421062/1‑150 (MQ=255)

                                      gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                           <  2:305259/150‑1 (MQ=255)

                                      gCACCAGATGAGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCAGATGCCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGCCGCACCAGCCATCGCTTTATTCGCGACTTGATCTGGCGTg                                                                                                           >  2:100551/1‑150 (MQ=255)

                                       cACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTTATCTGGCGTgg                                                                                                          <  1:981617/150‑1 (MQ=255)

                                           aGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGCCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATg                                                                                                      >  2:892061/1‑150 (MQ=255)

                                           aGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATg                                                                                                      <  1:547901/150‑1 (MQ=255)

                                            gATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                                                                     >  1:47726/1‑150 (MQ=255)

                                            gATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                                                                     <  1:224853/150‑1 (MQ=255)

                                                 gCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTg                                                                                                <  1:272265/150‑1 (MQ=255)

                                                   cTGATAGTGGTGGAAAATGTGATCGCCAGCGATGAGCTGATGACCAAATTGAGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATGGATTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                              <  1:264742/150‑1 (MQ=255)

                                                   cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCTAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                              >  2:897896/1‑150 (MQ=255)

                                                   cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                              >  2:987270/1‑150 (MQ=255)

                                                            gtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACCCCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGAACTGGTGAGCCTaaa                                                                                     >  1:1303483/1‑150 (MQ=255)

                                                             tggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGAGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAAc                                                                                    >  1:194634/1‑150 (MQ=255)

                                                             gggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAAc                                                                                    <  2:1578808/149‑1 (MQ=255)

                                                              ggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACt                                                                                   <  2:654383/150‑1 (MQ=255)

                                                                aaaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTc                                                                                 <  2:906081/150‑1 (MQ=255)

                                                                 aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTTATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                <  2:850586/150‑1 (MQ=255)

                                                                 aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                <  2:806871/150‑1 (MQ=255)

                                                                     gTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAACCTTCTGGaa                                                                            >  1:368947/1‑150 (MQ=255)

                                                                      tGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaa                                                                           >  1:923605/1‑150 (MQ=255)

                                                                       gATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaaa                                                                          <  2:773569/150‑1 (MQ=255)

                                                                         tCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATa                                                                        <  1:1637808/150‑1 (MQ=255)

                                                                            ccAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACg                                                                     <  1:1155798/150‑1 (MQ=255)

                                                                              aGCGATGAGCTTATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTTTGCGCCATTCATGGCTGACTCACCAGCCATCGCTTTATTCTCGTCTTGATCTGGCGTGGGATGTCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                   <  1:1613359/150‑1 (MQ=255)

                                                                                  aTGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGGGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATa                                                               >  1:143713/1‑150 (MQ=255)

                                                                                      gCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGcc                                                           <  2:187030/150‑1 (MQ=255)

                                                                                            gaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTc                                                     >  1:1599092/1‑150 (MQ=255)

                                                                                            gaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGAGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTc                                                     >  1:1605461/1‑150 (MQ=255)

                                                                                             aCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCa                                                    >  1:1192690/1‑150 (MQ=255)

                                                                                              ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                   <  1:501654/150‑1 (MQ=255)

                                                                                              ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                   >  2:1469220/1‑150 (MQ=255)

                                                                                                aaaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta                                                 >  1:629368/1‑150 (MQ=255)

                                                                                                 aaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtat                                                >  2:1442594/1‑150 (MQ=255)

                                                                                                  aTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtata                                               <  1:743090/150‑1 (MQ=255)

                                                                                                    tCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACg                                             <  1:531141/150‑1 (MQ=255)

                                                                                                        cATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggc                                         >  2:1353767/1‑150 (MQ=255)

                                                                                                              aaaaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGtt                                   >  2:1553664/1‑150 (MQ=255)

                                                                                                                   acacCTGGAGTTGTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTc                              <  2:908810/150‑1 (MQ=255)

                                                                                                                      ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGGCGTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCGTTCAGt                           >  2:1040383/1‑150 (MQ=255)

                                                                                                                            gTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATct                     >  1:1727090/1‑150 (MQ=255)

                                                                                                                            gTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATct                     >  1:1647728/1‑150 (MQ=255)

                                                                                                                               tggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggc                  >  2:413913/3‑150 (MQ=255)

                                                                                                                                ggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggct                 >  1:1279397/2‑150 (MQ=255)

                                                                                                                                 gtgCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctg                >  2:742032/1‑150 (MQ=255)

                                                                                                                                   gcgcCAGGCATGGCTGACGTACCAGCCATCGCTTTATTCGCGTCTTGATATGGCGTGGGATGGCACTGTTGATCCTAAACTTCTGGAATATATCGCGGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctggt              <  2:232476/150‑2 (MQ=255)

                                                                                                                                     taaagcgaTGGCTGGTGCGTCAGCCATCGCTTTTTTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAg            >  1:113565/8‑150 (MQ=255)

                                                                                                                                      ccAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGa           >  1:616668/1‑150 (MQ=255)

                                                                                                                                          cgaTGGCTGGTGCGTCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAg       >  1:757078/3‑150 (MQ=255)

                                                                                                                                           cATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGc      <  1:1614785/150‑1 (MQ=255)

                                                                                                                                            aTGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCt     >  1:845087/1‑150 (MQ=255)

                                                                                                                                             tGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCtt    >  1:970417/1‑150 (MQ=255)

                                                                                                                                              ggCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTa   >  1:399884/1‑150 (MQ=255)

                                                                                                                                               gCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTaa  >  1:146471/1‑150 (MQ=255)

GCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAA  >  NC_000913/3180547‑3180839

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