BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,637	T→C	7.7%	D467G (GAC→GGC) ‡	proS ←	proline‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,383,637	0	T	C	7.7%	437.0 / 17.5	143	D467G (GAC→GGC) ‡	proS	proline‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base T (68/64); new base C (5/6); total (73/70)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 152 total aligned reads displayed.

CCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTACTTG  >  NC_002947/1383490‑1383786

ccTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGcgtcg                                                                                                                                                     <  2:1737767/150‑1 (MQ=255)

 cTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGcgtcgt                                                                                                                                                    >  1:3741279/1‑150 (MQ=255)

   gggCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                                  <  2:1559243/150‑1 (MQ=255)

   gggCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                                  <  1:1502492/150‑1 (MQ=255)

     gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                                <  1:5318672/150‑1 (MQ=255)

     gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                                <  1:3314146/150‑1 (MQ=255)

     gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                                <  2:2947362/150‑1 (MQ=255)

     gCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGACGCCCTGGCTCCGTTCCa                                                                                                                                                <  2:1224908/150‑1 (MQ=255)

      cTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGACGCCCTGGCTCCGTTCCAg                                                                                                                                               <  1:827267/150‑1 (MQ=255)

       tGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAga                                                                                                                                              <  1:1715723/150‑1 (MQ=255)

       tGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACAAGGGCATCATCTGGAACGACGCCCTGGCTCCGTTCCa                                                                                                                                                <  2:4222271/150‑1 (MQ=37)

        ggTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAgat                                                                                                                                             >  2:2692468/1‑150 (MQ=255)

            ttCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGatgatg                                                                                                                                         <  2:1484833/150‑1 (MQ=255)

            ttCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGatgatg                                                                                                                                         <  1:2715087/150‑1 (MQ=255)

               ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                                      >  1:5255654/1‑150 (MQ=255)

               ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                                      >  1:4490727/1‑150 (MQ=255)

               ttGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                                      >  1:2508902/1‑150 (MQ=255)

                  tCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTg                                                                                                                                   >  2:1525087/1‑150 (MQ=255)

                  tCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTg                                                                                                                                   >  2:4697224/1‑150 (MQ=255)

                       gtcgtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTgtcgtc                                                                                                                              <  1:1306490/150‑1 (MQ=255)

                       gtcgtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTgtcgtc                                                                                                                              <  1:2692468/150‑1 (MQ=255)

                        tcgtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCa                                                                                                                             <  2:5441408/150‑1 (MQ=255)

                          gtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATa                                                                                                                           <  2:3444021/150‑1 (MQ=255)

                           tcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAg                                                                                                                          <  1:1018650/150‑1 (MQ=255)

                           tcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAg                                                                                                                          <  1:1522201/150‑1 (MQ=255)

                             cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                                        >  1:2039338/1‑150 (MQ=255)

                             cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                                        >  1:2039333/1‑150 (MQ=255)

                             cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                                        >  1:1638657/1‑150 (MQ=255)

                             cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                                        <  1:373775/150‑1 (MQ=255)

                             cagcagCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                                        <  1:4569968/150‑1 (MQ=255)

                              agcagcACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGctc                                                                                                                       >  1:1062933/1‑150 (MQ=255)

                              agcagcACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGctc                                                                                                                       >  1:3903709/1‑150 (MQ=255)

                                     cTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGa                                                                                                                <  1:1469114/150‑1 (MQ=255)

                                      ttCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGAt                                                                                                               <  2:4490727/150‑1 (MQ=255)

                                      ttCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGAt                                                                                                               >  2:3274276/1‑150 (MQ=255)

                                        cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATgg                                                                                                             <  1:2301507/150‑1 (MQ=255)

                                        cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATgg                                                                                                             >  1:3549344/1‑150 (MQ=255)

                                         aTAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGc                                                                                                            >  1:3929203/1‑150 (MQ=255)

                                          tAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTc                                                                                                                  >  1:2613758/1‑143 (MQ=255)

                                          tAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTc                                                                                                                  <  2:2613758/143‑1 (MQ=255)

                                          tAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCa                                                                                                           >  1:3114322/1‑150 (MQ=255)

                                          gaGCCGGCGGCAGTCAGGTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCa                                                                                                           <  2:5379809/149‑1 (MQ=255)

                                             ccGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCt                                                                                                        >  2:127315/1‑150 (MQ=255)

                                             ccGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCt                                                                                                        <  1:4697224/150‑1 (MQ=255)

                                                 cggcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCaa                                                                                                    <  2:2304494/150‑1 (MQ=255)

                                                     agtcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACcac                                                                                                >  1:3597791/1‑150 (MQ=255)

                                                       tcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGTCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCacac                                                                                              <  1:2864494/150‑1 (MQ=255)

                                                          gtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACggg                                                                                           >  2:1464488/1‑150 (MQ=255)

                                                           ttCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGa                                                                                          <  1:4729313/150‑1 (MQ=255)

                                                             cGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGaca                                                                                        <  1:1548886/150‑1 (MQ=255)

                                                               gCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACg                                                                                      >  2:3953545/1‑150 (MQ=255)

                                                                cATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                                     >  2:707630/1‑150 (MQ=255)

                                                                  tACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCa                                                                                   >  1:618993/1‑150 (MQ=255)

                                                                     aGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATg                                                                                <  2:3489641/150‑1 (MQ=255)

                                                                        ttGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                                  >  2:2966342/1‑81 (MQ=255)

                                                                        ttGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                                  <  1:2966342/81‑1 (MQ=255)

                                                                         tGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGt                                                                            >  1:5064523/1‑150 (MQ=255)

                                                                          gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                           <  2:1063863/150‑1 (MQ=255)

                                                                          gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                           <  2:1063938/150‑1 (MQ=255)

                                                                          gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                           <  1:1076557/150‑1 (MQ=255)

                                                                          gTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                           <  1:1542796/150‑1 (MQ=255)

                                                                             ggTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                                  >  2:3592779/1‑76 (MQ=255)

                                                                             ggTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTc                                                                                                                                                  <  1:3592779/76‑1 (MQ=255)

                                                                              gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCa                                                                                                           >  2:956429/1‑114 (MQ=255)

                                                                              gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCaa                                                                       >  1:85142/1‑150 (MQ=255)

                                                                              gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCaa                                                                       <  2:1898672/150‑1 (MQ=255)

                                                                              gTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCaa                                                                       >  2:3421668/1‑150 (MQ=255)

                                                                               ttGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGACGCCCTGGCTCCGTTCCAg                                                                                                                                               <  2:5085364/77‑1 (MQ=25)

                                                                                tGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAAcc                                                                     >  2:310400/1‑150 (MQ=255)

                                                                                  cTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCa                                                                   <  2:2478587/150‑1 (MQ=255)

                                                                                  cTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCa                                                                   >  2:4281184/1‑150 (MQ=255)

                                                                                  cTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCa                                                                   >  2:4281205/1‑150 (MQ=255)

                                                                                   ttCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCAt                                                                  >  2:3091746/1‑150 (MQ=255)

                                                                                   ttCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCAt                                                                  >  2:788517/1‑150 (MQ=255)

                                                                                   ttCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCAt                                                                  >  2:4431718/1‑150 (MQ=255)

                                                                                       cgGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGAc                                                              <  2:4795710/150‑1 (MQ=255)

                                                                                           ccACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTa                                                          <  2:4998049/150‑1 (MQ=255)

                                                                                            cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTAc                                                         >  1:4271612/1‑150 (MQ=255)

                                                                                                    tttCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCt                                                 >  1:1009997/1‑150 (MQ=255)

                                                                                                                   ggACCAGGGCA‑‑ATCTGGACCGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCa                                  <  1:1827634/150‑1 (MQ=255)

                                                                                                                   ggACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCa                                  <  1:5089824/150‑1 (MQ=255)

                                                                                                                    gACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGGCACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAg                                 >  2:3852251/1‑150 (MQ=255)

                                                                                                                     aCCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGc                                >  2:1906984/1‑150 (MQ=255)

                                                                                                                        aGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                                     <  2:3365777/94‑1 (MQ=255)

                                                                                                                        aGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                                     >  1:3365777/1‑94 (MQ=255)

                                                                                                                         gggCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATc                                                                 >  2:1051750/1‑113 (MQ=255)

                                                                                                                                aTCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCAc                       >  2:1732310/1‑150 (MQ=255)

                                                                                                                                 tCTGGAACGACGCCCTGGCTCCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATTGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACt                      >  1:4169430/1‑150 (MQ=255)

                                                                                                                                 tCTGGAACGACGCCCTGGCTCCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTAc                                                         <  1:3041694/115‑1 (MQ=38)

                                                                                                                                 tCTGGAACGACGCCCTGGCTCCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACt                      >  1:1944458/1‑150 (MQ=255)

                                                                                                                                  cTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACtt                     >  2:1998072/1‑150 (MQ=255)

                                                                                                                                       aCGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAgcg                >  1:5033493/1‑150 (MQ=255)

                                                                                                                                            gCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                           <  2:2043373/86‑1 (MQ=255)

                                                                                                                                            gCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCg           >  2:4005671/1‑150 (MQ=255)

                                                                                                                                               aGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTg        >  2:413041/1‑150 (MQ=255)

                                                                                                                                                  gCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTAc     <  2:3864625/150‑1 (MQ=255)

                                                                                                                                                    gtcgtTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTACtt   <  1:3274276/150‑1 (MQ=255)

                                                                                                                                                    gtcgtTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTACtt   <  2:3741279/150‑1 (MQ=255)

                                                                                                                                                     tcgtTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTACTTg  <  2:4175709/150‑1 (MQ=255)

                                                                                                                                                     tcgtTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATAGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTACTTg  <  2:4124003/150‑1 (MQ=255)

CCTGGGCTGGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCA‑‑ATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCTGTACTTG  >  NC_002947/1383490‑1383786

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