BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	218,544	C→A	5.9%	D78Y (GAT→TAT)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	218,544	0	C	A	5.9%	425.0 / 2.5	145	D78Y (GAT→TAT)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base C (66/68); new base A (2/9); total (68/77)
Fisher's exact test for biased strand distribution p-value = 6.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

CTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGATCATCCC  >  NC_000913/218395‑218688

cTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATc                                                                                                                                                  >  2:882298/1‑150 (MQ=255)

atGTTTGTAAGAGGTTTTGTCTTTAAGAATATGGTCATTTATAAGTTCTTCATGAGTTGGTCCGAGTACGAACGGTCGCTCGCCACGGTCAACAAAACGCTGCAGTTCCGGACCGTACTGTTCCCAACGAGCACTCTCTTGCCACAAATc                                                                                                                                                  <  2:1656302/149‑1 (MQ=255)

  gTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCgg                                                                                                                                                <  1:661568/150‑1 (MQ=255)

  gTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCAAGAGTTGGGACGAGGACGAACGGACGCTCGCCACGGTCACCAAAACGCAGCAGTTCCGGACCGTACTGTTCACAACGCCCACTCTCTTGCCACAAATCgg                                                                                                                                                >  2:137635/1‑150 (MQ=255)

   tttGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACTCTCGCCACGGTCAACAAAACGAAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggc                                                                                                                                               <  1:729279/150‑1 (MQ=255)

   tttGTAAGAGCTAAGCTCGTTACGAATCAGGGCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggc                                                                                                                                               <  1:272558/150‑1 (MQ=255)

      gTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggctgg                                                                                                                                            >  2:277931/1‑150 (MQ=255)

      gTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggctgg                                                                                                                                            <  1:1998178/150‑1 (MQ=255)

       tAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggctggc                                                                                                                                           <  1:1989496/150‑1 (MQ=255)

       tAAGAGCTAAGCGCGTTACGAATCAGGTCAGTGATACCTTCTTCACGAGTTGGGCCGAGTACGACCGGACGCGCGCCCCGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggctggc                                                                                                                                           <  1:1277096/150‑1 (MQ=255)

        aaGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggctggct                                                                                                                                          >  2:2007938/1‑150 (MQ=255)

          gagCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGa                                                                                                                                        <  2:769523/150‑1 (MQ=255)

          gagCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGa                                                                                                                                        >  2:386242/1‑150 (MQ=255)

           agCTACGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGCCCGTACTGTTCCCACCGCCCACTGGCTTGCCCCAAATCGGATGGCTGaa                                                                                                                                       >  2:589318/1‑150 (MQ=255)

           agCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAACTCGGCTGGCTGaa                                                                                                                                       <  1:1201016/150‑1 (MQ=255)

           agCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGaa                                                                                                                                       <  2:954806/150‑1 (MQ=255)

           agCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGaa                                                                                                                                       >  2:143488/1‑150 (MQ=255)

            gCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCCCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCACCGACCACTCTCTTGCCACAACTCGGCTGGCTGAac                                                                                                                                      >  2:563708/1‑150 (MQ=255)

            gCTAAGCGCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAac                                                                                                                                      <  1:424315/150‑1 (MQ=255)

                aGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAaccacc                                                                                                                                  <  1:1550582/150‑1 (MQ=255)

                   tCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGc                                                                                                                               <  2:1499266/150‑1 (MQ=255)

                       tACGAATCAGGGCAGTGATAACTTCTTCATGAGTTGGGCCTAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCg                                                                                                                           <  2:1528037/150‑1 (MQ=255)

                         cGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGACCGGACGCTCGCCCCGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCAGGCTGAACCACCGGCCTCGac                                                                                                                         >  2:1156188/1‑150 (MQ=255)

                             tCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCt                                                                                                                     <  1:2088314/150‑1 (MQ=255)

                               aGGTCAGTGATAAGTTATTCATGAGTTGGGCCGAGTACGAACGGACGCTCTCCACGGGCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCg                                                                                                                   <  1:1263525/150‑1 (MQ=255)

                                  tCAGTGATTTCTTCTTTTTGAGTTGGGCCGAGTACGAACGGACGCTCGGCTCTGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGATAGAAATAGGCTGGCTGAACCACCGGGATCGACACCTCGATc                                                                                                                <  2:1704601/150‑1 (MQ=255)

                                  tCAGTGATAACTTCTTCATGAGGTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATc                                                                                                                <  2:480092/150‑1 (MQ=255)

                                     gTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCa                                                                                                             >  1:191561/1‑150 (MQ=255)

                                      tGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCAc                                                                                                            <  1:1554866/150‑1 (MQ=255)

                                           aCTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCg                                                                                                       <  2:1773533/150‑1 (MQ=255)

                                           aCTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCg                                                                                                       <  1:789082/150‑1 (MQ=255)

                                              tcttcATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCgttg                                                                                                    <  1:544173/150‑1 (MQ=255)

                                               tttcATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCgttgt                                                                                                   <  2:2049618/149‑1 (MQ=255)

                                                ttcCTGAGTTGTGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCATCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTTTTTTTTTTATGATTGGTGTACCACCGTGATCTACACCTCGATCGCACCGGCgggggg                                                                                                  <  2:976189/150‑6 (MQ=255)

                                                ttcATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCgttgtt                                                                                                  <  1:698477/150‑1 (MQ=255)

                                                 tcATGAGTTGGGCCGAGTACGAACTGACGCTAGCCACGGTCAACAAAACGCAGCAGTTCCGGACAGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTc                                                                                                 <  1:272025/150‑1 (MQ=255)

                                                  cATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTcc                                                                                                >  1:905640/1‑149 (MQ=255)

                                                  cATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCa                                                                                                >  1:1155429/1‑150 (MQ=255)

                                                   aTGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCAGGACCGTACTGTTCCCAACGACCACTCTCTTGCATAAAATAGGCTTGCTGTACCACCGGGATCTACACCTCGATCGCACCGGCGTTGTTCAt                                                                                               <  2:1297490/150‑1 (MQ=255)

                                                       gTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATctct                                                                                           <  1:864419/150‑1 (MQ=255)

                                                       gTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATctct                                                                                           >  1:1465557/1‑150 (MQ=255)

                                                          gggCCGAGTACGCACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGTACCACCGTCATCTACACCTCGATCGCACCGGCGTTGTTCATCTCTTca                                                                                        <  2:1390224/150‑1 (MQ=255)

                                                          gggCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAGTCGGCTGGCTGAACCACCGGCATCGCCCCCTCGATCGCACCGGCGTGGTGCATCTCTTcc                                                                                        >  2:1122837/1‑149 (MQ=255)

                                                          gggCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGCTCGCACCGGCGTTGGTCATCTCTTcc                                                                                        >  2:187418/1‑149 (MQ=255)

                                                          gggCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTca                                                                                        <  2:1972911/150‑1 (MQ=255)

                                                           ggCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTTTTTTTTTTATGATTGGGGTACCACCGTGATCTACACCTCGATCGCACCGGCGTTGTTCATCGCTGcac                                                                                       <  2:1424956/150‑1 (MQ=255)

                                                             ccGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTcacgc                                                                                     <  1:143488/150‑1 (MQ=255)

                                                              cGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCATcacgca                                                                                    >  1:2088183/1‑150 (MQ=255)

                                                                    cGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTGCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGt                                                                              <  1:1495926/150‑1 (MQ=255)

                                                                       aCGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGtttt                                                                           <  1:1469649/150‑1 (MQ=255)

                                                                        cGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTc                                                                          <  1:394330/150‑1 (MQ=255)

                                                                         ggACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCg                                                                         <  1:1846623/150‑1 (MQ=255)

                                                                          gACGCTCGCCTCGGTCAACAAAACGCAGGCGTGACGGACCGTACTGTTCCCAACGACTACTCTCTTGTTTTATAAATGATGGCTGTACCACCGGGATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGa                                                                        <  2:1303329/150‑1 (MQ=255)

                                                                          gACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGCCCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCGCTTCACGCACGATGTTTTCGa                                                                        >  2:1251640/1‑150 (MQ=255)

                                                                           aCGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGAc                                                                       <  2:1413183/150‑1 (MQ=255)

                                                                                cGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATAGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACttttt                                                                  >  2:2075370/1‑150 (MQ=255)

                                                                                  ccACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTc                                                                <  1:2087884/150‑1 (MQ=255)

                                                                                   aaCGGTCAACAAAACGCAGCAGGTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCa                                                               <  2:1131470/149‑1 (MQ=255)

                                                                                    aCGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAg                                                              <  1:1274945/150‑1 (MQ=255)

                                                                                     cGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCCCGCACGATGTTTTCGACTTTTTTCAGa                                                             >  1:846422/1‑150 (MQ=255)

                                                                                        tCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAAcg                                                          <  1:1943850/150‑1 (MQ=255)

                                                                                         caacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGCCCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAAcgc                                                         >  1:1500810/1‑150 (MQ=255)

                                                                                          aacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTTTTTTTTTTATGATTGGTGTACCACCGTGATCTACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAAcgcg                                                        <  2:10343/150‑1 (MQ=255)

                                                                                          aacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAAcgcg                                                        >  2:804372/1‑150 (MQ=255)

                                                                                          aacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTAATCTCTTCCCGCACGATGTTTTCGACTTTTTTCAGAAcgcg                                                        >  1:1030731/1‑150 (MQ=255)

                                                                                             aaaaCGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCAc                                                     >  2:807663/1‑150 (MQ=255)

                                                                                                     gcagTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCg                                             >  2:799265/1‑150 (MQ=255)

                                                                                                        gTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTGGCCTTTTTTCAGCTCGCGCACGCCGGTCGGcc                                          >  1:355574/1‑149 (MQ=255)

                                                                                                          tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGc                                        >  2:1788483/1‑150 (MQ=255)

                                                                                                          tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGc                                        <  1:934669/150‑1 (MQ=255)

                                                                                                          tCCGGACCGTAATCTTCCCAACGACCACACTCTTGCCACAAATAGGATGGCTGAACCACCCGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCTACTTTTTTCAGCACGCGCACGCCCGTCCGCAGc                                        >  2:618231/1‑150 (MQ=255)

                                                                                                           ccGGACCGTACTGTTCCCAACGCCCCCTCTCTTGCCACAAATCGGCTGGCTGAACCCCCGGCATCGACGCCGCGGTCGCCCCGGCGTTGTTCAGCTCTTGCCGCCCGAGGGTGTCGGCGTTGTTCCGGCCGCGCCCGCCGGTGGGCTGcc                                       >  2:728744/1‑150 (MQ=255)

                                                                                                           ccGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGCTCGCACCGGCGTTGTTCATCTCTTAACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGcc                                       >  2:164813/1‑150 (MQ=255)

                                                                                                           ccGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGAATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGcc                                       >  2:1682566/1‑150 (MQ=255)

                                                                                                             ggACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGAGCGCACCGGCGTAGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAg                                     <  1:91647/150‑1 (MQ=255)

                                                                                                              gACCGTACTGTTCCCAACGACCACTCTCTAGCCACAAATCGGCTGGCTGAACCACCGGCAGCGACACCTCGATCGCACCGGCGTTGATCATCTCTTCACGCACGATGTTTGCGACTTTTTTCAGAACGCGCACGCCGGTCTGCAGCCAgg                                    <  1:1362623/150‑1 (MQ=255)

                                                                                                                ccGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATGTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGta                                  <  2:1316298/150‑1 (MQ=255)

                                                                                                                ccGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCCCGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGta                                  >  2:755660/1‑150 (MQ=255)

                                                                                                                  gTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGTCTGAACCACCGTCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGtata                                <  2:1395660/150‑1 (MQ=255)

                                                                                                                   tACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGAAGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGTCAGCCAGGtatat                               >  1:314386/1‑150 (MQ=255)

                                                                                                                           ccAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGa                       >  2:492026/1‑150 (MQ=255)

                                                                                                                             aaCGACCACTCTCTTGCCACAAATCTGATGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAgg                     <  2:679535/150‑1 (MQ=255)

                                                                                                                                 aCCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAg                 <  2:1456843/150‑1 (MQ=255)

                                                                                                                                   cACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCt               <  1:1825258/150‑1 (MQ=255)

                                                                                                                                   cACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCt               <  1:28935/150‑1 (MQ=255)

                                                                                                                                   cACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCt               >  1:866813/1‑150 (MQ=255)

                                                                                                                                     ctctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTg             >  2:1207281/1‑150 (MQ=255)

                                                                                                                                     ctctctTGCCACAAATAGGCTCGCTGAACCACCGGCATCGACACATCGATCGCGCCGGCGTTGTTCATCTCTTCCCGCACGAAGTTTTGGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGAATATACCCCGGAGGCCAGCTTg             >  2:1743775/1‑150 (MQ=255)

                                                                                                                                        tctTGCCACTAATCGGCTGGCTGTACCACGGGCATCGACTCCTAGATCGCACGGGCGTTGTTCATCTCTTCACGCACGATGTTTTCTTCTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCgg          <  2:63356/150‑1 (MQ=255)

                                                                                                                                        tctTGCCACAAATCGGCTGGCTGAACCACCGGCAGCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCgg          <  2:598147/150‑1 (MQ=255)

                                                                                                                                          ttGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGat        <  1:2007938/150‑1 (MQ=255)

                                                                                                                                          ttGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGat        <  2:191561/150‑1 (MQ=255)

                                                                                                                                          ttGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGAGGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGat        <  1:1158870/150‑1 (MQ=255)

                                                                                                                                           tGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCTCGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGctc       >  1:430618/1‑150 (MQ=255)

                                                                                                                                           tGCCACAAATCGGCTGGCGGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGatc       <  1:623648/150‑1 (MQ=255)

                                                                                                                                            gCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGatca      >  2:504215/1‑150 (MQ=255)

                                                                                                                                             ccACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGACCGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGatcat     >  1:675105/1‑150 (MQ=255)

                                                                                                                                             acacaAAACGGCTGGCTTAACCACATGCATCGACAACTAGATCGCACAGGCCTTGTTAATATATTCACGCACGATGTTTTCGACTTTTTTCAGCACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGatcat     >  2:26099/2‑150 (MQ=255)

                                                                                                                                              cacaAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGAGGGTTTCGACGTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAAGCCAGCTTGCGGatcagc    >  1:1846145/1‑148 (MQ=255)

                                                                                                                                                caAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGATCATccc  >  1:523006/1‑150 (MQ=255)

CTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCACGCACGATGTTTTCGACTTTTTTCAGAACGCGCACGCCGGTCGGCAGCCAGGTATATAACCCGGAGGCCAGCTTGCGGATCATCCC  >  NC_000913/218395‑218688

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