BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,716,402	A→C	F274C (TTC→TGC)	tyrS ←	tyrosyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,716,402	0	A	C	91.6%	358.5 / 15.6	120	F274C (TTC→TGC)	tyrS	tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (5/5); new base C (61/48); total (67/53)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 121 total aligned reads displayed.

TGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCG  >  NC_000913/1716306‑1716497

tGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACt                                                                                              <  2:3613746/100‑1 (MQ=255)

 gCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACtt                                                                                             >  2:1195443/1‑100 (MQ=255)

  cTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACtag                                                                                            >  1:2590847/1‑98 (MQ=255)

    ccGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAg                                                                                          <  2:1866388/100‑1 (MQ=255)

     cGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAgg                                                                                         <  2:2753310/100‑1 (MQ=255)

     cGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAgg                                                                                         >  2:2528918/1‑100 (MQ=255)

     cGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAgg                                                                                         <  2:1971159/100‑1 (MQ=255)

       ccAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGaa                                                                                       >  1:2816198/1‑100 (MQ=255)

       ccAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGaa                                                                                       >  1:228613/1‑100 (MQ=255)

         aGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGc                                                                                     <  2:428644/100‑1 (MQ=255)

          gTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCg                                                                                    >  1:2006724/1‑100 (MQ=255)

           tACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCgg                                                                                   <  1:3611267/100‑1 (MQ=255)

            aCATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGt                                                                                  >  2:3528686/1‑100 (MQ=255)

            aCATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGt                                                                                  <  1:3968306/100‑1 (MQ=255)

            aCATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGt                                                                                  <  1:2765557/100‑1 (MQ=255)

            aCATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGt                                                                                  <  1:1087027/100‑1 (MQ=255)

               tACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTaaa                                                                               >  2:3041682/1‑100 (MQ=255)

                ncTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAAc                                                                              >  1:3243038/2‑100 (MQ=255)

                aCTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAAc                                                                              >  1:1773464/1‑100 (MQ=255)

                 cTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCCATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACg                                                                             >  1:897436/1‑100 (MQ=255)

                 cTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATCAAGGTGAAGCACTTCAGGAAGCGGTAAACg                                                                             >  1:267845/1‑100 (MQ=255)

                    ggCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCCCTTCAATGCCCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCg                                                                          <  1:3948305/100‑1 (MQ=255)

                            gTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGc                                                                  <  1:66003/100‑1 (MQ=255)

                            gTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGc                                                                  <  2:3192744/100‑1 (MQ=255)

                             tGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCa                                                                 >  1:3966188/1‑100 (MQ=255)

                              gCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAg                                                                >  1:2993841/1‑100 (MQ=255)

                              gCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAg                                                                >  1:1929035/1‑100 (MQ=255)

                               cTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAgt                                                               >  1:3158888/1‑100 (MQ=255)

                               cTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCAGCCGCAgt                                                               <  1:1922898/100‑1 (MQ=255)

                                tttACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCCgtg                                                              >  1:2072403/1‑100 (MQ=255)

                                 ttACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAgtgt                                                             >  1:2331016/1‑100 (MQ=255)

                                  tACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGtt                                                            >  2:845180/1‑100 (MQ=255)

                                    ccGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGa                                                          >  2:2230644/1‑100 (MQ=255)

                                    ccGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGa                                                          >  1:2132325/1‑100 (MQ=255)

                                    ccGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGa                                                          >  1:2243082/1‑100 (MQ=255)

                                     cGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGAt                                                         >  1:3146202/1‑100 (MQ=255)

                                      gCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATc                                                        >  1:726736/1‑100 (MQ=255)

                                      gCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATc                                                        >  2:2667768/1‑100 (MQ=255)

                                      gCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATc                                                        >  1:1760277/1‑100 (MQ=255)

                                       cTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATcc                                                       <  1:3963066/100‑1 (MQ=255)

                                       cTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGCTGCTcc                                                       >  2:1831864/1‑100 (MQ=255)

                                         gTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAg                                                     <  1:1151921/100‑1 (MQ=255)

                                          tttttATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATTCTTATAAAGGTGAAGCACTTCAGGATGCGGTAAACGTCGGCATCCGCAGTGTTGATCTAGa                                                    >  2:3944702/1‑100 (MQ=255)

                                          tttttATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGa                                                    <  2:1665486/100‑1 (MQ=255)

                                            tttATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCTAGAAc                                                  >  1:3832259/1‑100 (MQ=255)

                                               aTCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTgg                                               <  1:83388/100‑1 (MQ=255)

                                               aTCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTgg                                               <  1:114543/100‑1 (MQ=255)

                                                 cttcttcttCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTa                                             <  1:638827/100‑1 (MQ=255)

                                                 cttcttcttCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTa                                             <  2:1689427/100‑1 (MQ=255)

                                                  ttcttcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAg                                            <  2:3060496/100‑1 (MQ=255)

                                                  ttcttcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAg                                            <  2:1802859/100‑1 (MQ=255)

                                                  ttcttcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTAAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAg                                            >  1:1593103/1‑100 (MQ=255)

                                                   tcttcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCGTCCGCAGTGTTGATCCAGAACTGGTAGa                                           <  2:3701166/100‑1 (MQ=255)

                                                   tcttcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGa                                           >  2:138422/1‑100 (MQ=255)

                                                    cttcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGaa                                          >  2:2897044/1‑100 (MQ=255)

                                                      tcttcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGCGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAAtt                                        >  2:3906146/1‑100 (MQ=255)

                                                         tcCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGt                                     <  2:2658470/100‑1 (MQ=255)

                                                          ccAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGCCGGCAACCGCAGAGTTGATCCAGAACTGGTAGAATTTGTa                                    >  2:3818879/1‑100 (MQ=255)

                                                           cAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTAc                                   >  1:1796247/1‑100 (MQ=255)

                                                           cAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTAc                                   <  2:3404702/100‑1 (MQ=255)

                                                            aGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACg                                  <  1:693218/100‑1 (MQ=255)

                                                                     aTCTCTTCAATGCTCATAAAGGTGAAGCACTTTAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGttt                         <  2:732823/100‑1 (MQ=255)

                                                                     aTCTCTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGttt                         >  2:1530050/1‑100 (MQ=255)

                                                                      tctcTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGtttt                        >  1:3643132/1‑100 (MQ=255)

                                                                      tctcTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGtttt                        >  1:521985/1‑100 (MQ=255)

                                                                      tctcTTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGtttt                        >  1:2790867/1‑100 (MQ=255)

                                                                       ctctTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTttc                       <  2:3073477/100‑1 (MQ=255)

                                                                        tctTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTttct                      >  1:2390954/1‑100 (MQ=255)

                                                                        tctTAAATGCTCGTAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTttct                      <  2:2041742/100‑1 (MQ=255)

                                                                         ctTCAATGCTGATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTttctt                     >  1:3574832/1‑100 (MQ=255)

                                                                         ctTCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTttctt                     <  2:3311618/100‑1 (MQ=255)

                                                                           tCAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCg                   <  1:3339209/100‑1 (MQ=255)

                                                                           tCAATGCTCATAAAGGTGAAGCACGTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTAATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCg                   <  1:1366319/100‑1 (MQ=255)

                                                                            cAATGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCgn                  <  1:3254364/100‑2 (MQ=255)

                                                                             aaTGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGa                 <  2:2415608/100‑1 (MQ=255)

                                                                             aaTGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGa                 <  2:3400277/100‑1 (MQ=255)

                                                                             aaTGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGa                 >  2:847440/1‑100 (MQ=255)

                                                                              aTGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGAt                <  2:3689741/100‑1 (MQ=255)

                                                                              aTGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGAt                <  2:3630288/100‑1 (MQ=255)

                                                                               tGCTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATc               <  1:1273046/100‑1 (MQ=255)

                                                                                 cTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCa             >  2:1602871/1‑100 (MQ=255)

                                                                                 cTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCa             >  1:549367/1‑100 (MQ=255)

                                                                                 cTCATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTCTTGATCCAGAACGAGAAGAATTTGTACGGGCTGGTTTTCTTCGGATCCa             >  2:3936647/1‑100 (MQ=255)

                                                                                   cATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAAc           >  1:562862/1‑100 (MQ=255)

                                                                                   cATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAAc           <  1:2569431/100‑1 (MQ=255)

                                                                                   cATAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAAc           >  1:2286338/1‑100 (MQ=255)

                                                                                    aTAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAAcc          <  2:78496/100‑1 (MQ=255)

                                                                                    aTAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAAcc          <  1:1344166/100‑1 (MQ=255)

                                                                                    aTAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAAcc          >  1:976007/1‑100 (MQ=255)

                                                                                     tAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCa         >  2:2696623/1‑100 (MQ=255)

                                                                                     tAAAGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCa         >  1:2821224/1‑100 (MQ=255)

                                                                                     tAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTATGGGCTGGTTTTCTTCGGATCCAACCa         >  2:2822194/1‑100 (MQ=255)

                                                                                       aaGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGa       >  1:3642658/1‑100 (MQ=255)

                                                                                       aaGGTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGa       >  2:931554/1‑100 (MQ=255)

                                                                                         ggTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACt     >  1:2616100/1‑100 (MQ=255)

                                                                                         ggTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACt     >  1:3554660/1‑100 (MQ=255)

                                                                                         ggTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACt     >  1:4067852/1‑100 (MQ=255)

                                                                                         ggTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACt     >  1:58031/1‑100 (MQ=255)

                                                                                          gTGAAGCACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTg    <  2:2272122/100‑1 (MQ=255)

                                                                                            gaagcaCTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCTAGAACTGGGAGAATTTGTACGGGCTGGTATTCTTCGGATCCAACCAGACTgcg  >  2:2420494/1‑100 (MQ=255)

TGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCTTCTTCTTCCAGGGCGTTGATCTCTTCAATGCTCATAAAGGTGAAGAACTTCAGGAAGCGGTAAACGTCGGCATCCGCAGTGTTGATCCAGAACTGGTAGAATTTGTACGGGCTGGTTTTCTTCGGATCCAACCAGACTGCG  >  NC_000913/1716306‑1716497

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