BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,193,342	G→T	5.8%	R82L (CGT→CTT) ‡	PP_3684 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,193,342	0	G	T	5.8%	938.9 / 34.4	276	R82L (CGT→CTT) ‡	PP_3684	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (128/132); new base T (7/9); total (135/141)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 278 total aligned reads displayed.

CGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGAC  >  NC_002947/4193193‑4193491

 gTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGACGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                      >  1:100289/1‑150 (MQ=255)

 gTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                      <  1:535219/150‑1 (MQ=255)

   cTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGa                                                                                                                                                    <  2:3379749/150‑1 (MQ=255)

   cTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGa                                                                                                                                                    <  2:1825854/150‑1 (MQ=255)

    ttCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAg                                                                                                                                                   >  1:993075/1‑150 (MQ=255)

    gtCCCCGACCGCAGACCAGACGCAGAGATACGGCATGAGCGCCAATTCATCGAGCCGGTGAACAGCCTCAAACCCTCGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGCGATACGTGATGTGATCCTTGAGAGCCGTGAg                                                                                                                                                   <  1:6927260/149‑1 (MQ=255)

       ccAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGct                                                                                                                                                >  1:3360781/1‑150 (MQ=255)

        cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                               >  1:5208710/1‑150 (MQ=255)

        cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                               >  1:919731/1‑150 (MQ=255)

        cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                               >  1:1080271/1‑150 (MQ=255)

        cAGACCGAAGACCAGAAGCAGAGATAAGGCATCAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTCAGGGTCACAGACATAGGCGGGATGTCGATCGACCCTGAACAGATACGTGATGTGATCCAAGAGAGCCGTGAGGctc                                                                                                                                               >  2:5650545/1‑150 (MQ=255)

         aGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctct                                                                                                                                              >  1:6470084/1‑150 (MQ=255)

           aCCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCa                                                                                                                                            >  1:556363/1‑150 (MQ=255)

            ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCaa                                                                                                                                           <  2:5188820/150‑1 (MQ=255)

             cGAAGACCAGAGGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAg                                                                                                                                          <  1:542804/150‑1 (MQ=255)

             cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                          <  1:2912028/150‑1 (MQ=255)

              gAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAgg                                                                                                                                         <  1:2413197/150‑2 (MQ=255)

              gAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGc                                                                                                                                         >  2:6801613/1‑150 (MQ=255)

                aGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCAc                                                                                                                                       <  1:6990923/150‑1 (MQ=255)

                   ccAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTg                                                                                                                                    >  1:2596362/1‑150 (MQ=255)

                   ccAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTg                                                                                                                                    <  2:5347638/150‑1 (MQ=255)

                    cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                                                                                                                   >  2:3304291/1‑150 (MQ=255)

                    cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                                                                                                                   >  2:977819/1‑150 (MQ=255)

                     aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                                                                                                                  >  1:799832/1‑150 (MQ=255)

                     aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                                                                                                                  <  1:3687173/150‑1 (MQ=255)

                      gAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                                                                                                                 <  1:4442657/150‑1 (MQ=255)

                      gAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                                                                                                                 <  1:4028987/150‑1 (MQ=255)

                      gAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                                                                                                                 <  1:6141017/150‑1 (MQ=255)

                       aaGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCg                                                                                                                                >  1:155007/1‑150 (MQ=255)

                             agaTAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCgg                                                                                                                          >  2:4995645/1‑150 (MQ=255)

                               aTAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgcg                                                                                                                        <  2:1430718/150‑1 (MQ=255)

                                tAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgcgc                                                                                                                       >  2:5453246/1‑150 (MQ=255)

                                       tGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaa                                                                                                                >  1:2518971/1‑150 (MQ=255)

                                        gAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                               <  1:6701603/150‑1 (MQ=255)

                                        gAGCGCCAATTCAGCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                               <  2:2782189/150‑1 (MQ=255)

                                          gcgcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGc                                                                                                             >  2:240300/1‑150 (MQ=255)

                                           cgcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCg                                                                                                            <  1:2337786/150‑1 (MQ=255)

                                           cgcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCg                                                                                                            >  2:4034498/1‑150 (MQ=255)

                                            gcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                          >  1:4149143/1‑119 (MQ=255)

                                            gcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGa                                                                                                           >  1:1714586/1‑150 (MQ=255)

                                            gcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGAACTTGGTCGACACGGCGCATCGGAAAGCGa                                                                                                           >  2:331529/1‑150 (MQ=255)

                                             ccAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGaa                                                                                                          <  2:7325528/150‑1 (MQ=255)

                                             ccAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGaa                                                                                                          <  1:7104806/150‑1 (MQ=255)

                                              cAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAAc                                                                                                         <  1:6534410/150‑1 (MQ=255)

                                              cAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAAc                                                                                                         >  1:5454765/1‑150 (MQ=255)

                                                aTTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTAGGTCGACCCGGCGCATCGGAAAGCGAACaa                                                                                                       >  1:2212889/1‑150 (MQ=255)

                                                   cATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAAc                                                                                                    <  2:5758349/150‑1 (MQ=255)

                                                       gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                                <  2:2962416/150‑1 (MQ=255)

                                                       gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                                <  2:3088603/150‑1 (MQ=255)

                                                                aaCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCtt                                                                                       <  1:5031122/150‑1 (MQ=255)

                                                                aaCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCtt                                                                                       <  2:6873958/150‑1 (MQ=255)

                                                                 aCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACt                                                                                                                                      <  2:4552175/102‑1 (MQ=255)

                                                                     ccTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGt                                                                                  >  1:6636747/1‑150 (MQ=255)

                                                                       tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                          <  2:2493274/92‑1 (MQ=39)

                                                                       tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGttt                                                                                >  2:6004998/1‑150 (MQ=255)

                                                                       tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGttt                                                                                >  1:4297026/1‑150 (MQ=255)

                                                                       tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGttt                                                                                >  1:3295220/1‑150 (MQ=255)

                                                                       tCAAAACATTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTAGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTAGTGGtta                                                                                >  1:6875632/1‑149 (MQ=255)

                                                                        cAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTc                                                                               >  1:5809692/1‑150 (MQ=255)

                                                                         aaaaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCa                                                                              <  2:4000540/150‑1 (MQ=255)

                                                                              cTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGAtct                                                                         >  2:6242690/1‑150 (MQ=255)

                                                                                 gtgggTCACTGACCGAGTCGGGATGTCGATCGACCCTGAAGAGTTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCaa                                                                      <  2:1224592/148‑1 (MQ=255)

                                                                                  aGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAg                                                                     <  2:2041920/150‑1 (MQ=255)

                                                                                      tCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgc                                                                                                                         >  1:4796540/1‑94 (MQ=255)

                                                                                      tCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgc                                                                                                                         <  2:4796540/94‑1 (MQ=255)

                                                                                         cTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTa                                                              <  2:6542863/150‑1 (MQ=255)

                                                                                          tGACCGAGGCGGGATGTCGATCGACCCTGAAGAGAGACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTGGTGGTTTCAGATCTCAAGCCCTCTac                                                             >  2:4745552/1‑149 (MQ=255)

                                                                                            aCCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAAc                                                           >  1:7240866/1‑150 (MQ=255)

                                                                                            aCCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAAc                                                           >  1:364342/1‑150 (MQ=255)

                                                                                              cGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTg                                                         >  2:1174397/1‑150 (MQ=255)

                                                                                               gAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGc                                                        <  1:1778054/150‑1 (MQ=255)

                                                                                                aGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCa                                                       <  2:4685639/150‑1 (MQ=255)

                                                                                                    gggATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGa                                                   >  2:2728204/1‑150 (MQ=255)

                                                                                                     ggATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGAc                                                  <  1:3549452/150‑1 (MQ=255)

                                                                                                      gATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACt                                                 >  1:4777271/1‑150 (MQ=255)

                                                                                                          tcgatcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCaa                                             <  2:4111345/150‑1 (MQ=255)

                                                                                                            gatcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACt                                           <  1:1249201/150‑1 (MQ=255)

                                                                                                            gatcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACt                                           >  1:4592788/1‑150 (MQ=255)

                                                                                                            gatcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACt                                           >  1:290884/1‑150 (MQ=255)

                                                                                                              tcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGc                                         >  2:3545334/1‑150 (MQ=255)

                                                                                                                  cccTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTAt                                     >  1:6861684/1‑150 (MQ=255)

                                                                                                                   ccTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACAGCGTAtt                                    >  2:7178716/1‑150 (MQ=255)

                                                                                                                     agAAGACATACGTGATGTGATCCAAGAGAGCCGAGACGCACTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAACCCCTCTAAACTGCAGCGACTTCAACTGCGTAtttt                                  <  2:7415495/149‑1 (MQ=255)

                                                                                                                         gagaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGc                              <  1:6801613/150‑1 (MQ=255)

                                                                                                                           gaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGccc                            <  1:4880508/150‑1 (MQ=255)

                                                                                                                              aCGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCttt                         >  1:7184120/1‑150 (MQ=255)

                                                                                                                               cGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTg                        <  2:5574049/150‑1 (MQ=255)

                                                                                                                                   atgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCt                    <  1:5884213/150‑1 (MQ=255)

                                                                                                                                    tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                   <  1:410672/150‑1 (MQ=255)

                                                                                                                                    tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                   <  1:1926142/150‑1 (MQ=255)

                                                                                                                                       gatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGTGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGAc                >  1:6333206/1‑150 (MQ=255)

                                                                                                                                         tgcTTGAGAGCCTCACGGCTCTCAAGGACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACaa              >  1:3460032/3‑150 (MQ=255)

                                                                                                                                          ccTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAAc             >  2:542847/1‑150 (MQ=255)

                                                                                                                                           cTTGAGAGCCTCACGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACt            >  2:6720731/1‑150 (MQ=255)

                                                                                                                                             tGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTcc          <  2:3751582/150‑1 (MQ=255)

                                                                                                                                             tGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTcc          <  1:2053980/150‑1 (MQ=255)

                                                                                                                                             tGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTcc          <  2:6967909/150‑1 (MQ=255)

                                                                                                                                              gagaGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCg         >  2:1059738/1‑150 (MQ=255)

                                                                                                                                                gagCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGta       <  2:4199157/150‑1 (MQ=255)

                                                                                                                                                  gCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGtata     <  2:477283/150‑1 (MQ=255)

CGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGAC  >  NC_002947/4193193‑4193491

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