BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	645,034	A→G	V257A (GTG→GCG)	ybdS ←	YbdS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	645,034	0	A	G	99.6%	741.5 / ‑4.7	231	V257A (GTG→GCG)	ybdS	YbdS
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (116/114); total (116/115)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 239 total aligned reads displayed.

TCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATACCG  >  NC_000913/644944‑645123

tctctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGcgcg                                                                                         >  2:1267752/1‑93 (MQ=255)

tctctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGcgcg                                                                                         >  2:4254542/1‑93 (MQ=255)

 ctctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGt                                                                                        <  1:1382718/93‑1 (MQ=255)

  tctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTa                                                                                       <  2:5688389/93‑1 (MQ=255)

  tctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTa                                                                                       >  1:1238700/1‑93 (MQ=255)

   ctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAg                                                                                      >  2:1763246/1‑93 (MQ=255)

    tGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAgg                                                                                     <  2:5854289/93‑1 (MQ=255)

    tGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGGAgg                                                                                     >  2:2834376/1‑93 (MQ=255)

      cGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaa                                                                                   <  2:204152/93‑1 (MQ=255)

       gtgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaaa                                                                                  >  1:3793366/1‑93 (MQ=255)

       gtgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaaa                                                                                  >  2:3022177/1‑93 (MQ=255)

       gtgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaaa                                                                                  >  1:2487054/1‑93 (MQ=255)

         gtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGc                                                                                <  2:1109456/93‑1 (MQ=255)

          tCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGcc                                                                               <  2:1375773/93‑1 (MQ=255)

           cAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCa                                                                              <  2:6335744/93‑1 (MQ=255)

           cAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCa                                                                              <  2:866318/93‑1 (MQ=255)

            aGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCAc                                                                             >  1:758203/1‑93 (MQ=255)

            aGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCAc                                                                             <  2:290049/93‑1 (MQ=255)

            aGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTACGAAAGCCAc                                                                             >  2:1602208/1‑93 (MQ=255)

             gcgcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACg                                                                            >  2:694684/1‑93 (MQ=255)

             gcgcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACg                                                                            <  1:5250489/93‑1 (MQ=255)

              cgcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACgg                                                                           <  1:2577817/93‑1 (MQ=255)

              cgcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACgg                                                                           <  2:4103555/93‑1 (MQ=255)

                cACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGgcg                                                                         >  1:1941221/1‑93 (MQ=255)

                 aCCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGgcgc                                                                        >  1:5767947/1‑93 (MQ=255)

                   ccATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCaa                                                                      >  2:6478228/1‑93 (MQ=255)

                   ccATGGTTTTAAGTTCATCACCCGCCAAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCaa                                                                      >  1:3373782/1‑93 (MQ=255)

                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                     >  2:4324319/1‑93 (MQ=255)

                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                     >  1:5111434/1‑93 (MQ=255)

                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGCTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                     >  1:3265672/1‑93 (MQ=255)

                      tGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGa                                                                   >  2:361243/1‑93 (MQ=255)

                       ggTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGAt                                                                  <  1:5394892/93‑1 (MQ=255)

                       ggTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGAt                                                                  <  2:466347/93‑1 (MQ=255)

                       ggTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGAt                                                                  <  2:4968443/93‑1 (MQ=255)

                        gTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATa                                                                 >  1:1561119/1‑93 (MQ=255)

                         ttttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATaa                                                                >  1:1862322/1‑93 (MQ=255)

                          tttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAg                                                               >  2:4433309/1‑93 (MQ=255)

                           ttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGc                                                              <  1:4634316/93‑1 (MQ=255)

                              aGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAg                                                           <  2:1510110/93‑1 (MQ=255)

                              aGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAg                                                           >  2:932511/1‑93 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          >  2:5933572/1‑93 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          <  2:386981/93‑1 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          <  2:5780677/93‑1 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          <  1:2362404/93‑1 (MQ=255)

                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                         >  1:5625477/1‑93 (MQ=255)

                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                         <  2:416632/93‑1 (MQ=255)

                                   atcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAAc                                                      <  2:356157/93‑1 (MQ=255)

                                    tcccccGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAAcc                                                     <  2:2777251/90‑1 (MQ=255)

                                    tcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAAcc                                                     <  2:6566525/93‑1 (MQ=255)

                                      aCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAAcccc                                                   >  2:5122700/1‑93 (MQ=255)

                                       cccGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCa                                                  <  1:300935/93‑1 (MQ=255)

                                       cccGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCa                                                  >  1:2018609/1‑93 (MQ=255)

                                       cccGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCa                                                  >  2:5266695/1‑93 (MQ=255)

                                           cccAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCg                                              >  1:4108316/1‑93 (MQ=255)

                                            ccAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCgg                                             >  1:5016804/1‑93 (MQ=255)

                                             cAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGc                                            >  1:6099101/1‑93 (MQ=255)

                                              aGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCa                                           <  1:5455746/93‑1 (MQ=255)

                                              aGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCa                                           <  1:4539754/93‑1 (MQ=255)

                                              aGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCa                                           >  1:2983831/1‑93 (MQ=255)

                                                gTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAgg                                         <  2:5912457/93‑1 (MQ=255)

                                                gTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAgg                                         >  2:2350872/1‑93 (MQ=255)

                                                 ttGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGa                                        <  2:2430639/93‑1 (MQ=255)

                                                  tGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGaa                                       <  1:6480009/93‑1 (MQ=255)

                                                     cACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCa                                    <  1:2147841/93‑1 (MQ=255)

                                                      aCTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCag                                   >  1:3916572/1‑93 (MQ=255)

                                                      aCTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCag                                   >  1:6282749/1‑93 (MQ=255)

                                                        ttcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCagag                                 <  1:1015899/93‑1 (MQ=255)

                                                         tcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAgg                                <  2:3302096/93‑1 (MQ=255)

                                                         tcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAgg                                <  2:3714484/93‑1 (MQ=255)

                                                           ttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGaa                              >  2:6203041/1‑93 (MQ=255)

                                                             cACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAAcc                            >  1:2793038/1‑93 (MQ=255)

                                                             cACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAAcc                            >  1:2192869/1‑93 (MQ=255)

                                                              aCTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCa                           <  1:3622770/93‑1 (MQ=255)

                                                               cTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCAc                          <  1:342445/93‑1 (MQ=255)

                                                               cTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCAc                          <  2:4583230/93‑1 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                         >  1:1181008/1‑93 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                         <  1:1413133/93‑1 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                         <  1:483117/93‑1 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                         <  1:4876659/93‑1 (MQ=255)

                                                                 gtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGGACCACTg                        >  1:639783/1‑93 (MQ=255)

                                                                  tgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGCACCACTGc                       <  1:6570176/93‑1 (MQ=255)

                                                                  tgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGc                       >  2:2581953/1‑93 (MQ=255)

                                                                  tgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGc                       <  1:419200/93‑1 (MQ=255)

                                                                  tgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGc                       >  1:988896/1‑93 (MQ=255)

                                                                       gATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCa                  <  2:3723482/93‑1 (MQ=255)

                                                                        aTTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAg                 <  2:135900/93‑1 (MQ=255)

                                                                           tCCGGTATGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTa              <  1:6232764/93‑1 (MQ=255)

                                                                            ccGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATACCCCCAGCCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTaa             <  2:1410198/93‑1 (MQ=255)

                                                                             cGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAAt            >  2:107327/1‑93 (MQ=255)

                                                                              ggTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATc           <  1:3556655/93‑1 (MQ=255)

                                                                              ggTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAAACAGCTAATc           >  2:4730589/1‑93 (MQ=255)

                                                                                   gTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGAt      >  1:785664/1‑93 (MQ=255)

                                                                                   gTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGAt      >  2:6411724/1‑93 (MQ=255)

                                                                                   gTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGAt      <  1:3853624/93‑1 (MQ=255)

                                                                                    tACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATa     >  2:5523054/1‑93 (MQ=255)

                                                                                    tACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATa     >  1:5802583/1‑93 (MQ=255)

                                                                                     aCAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATAc    >  2:1535437/1‑93 (MQ=255)

                                                                                     aCAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATAc    >  1:1747115/1‑93 (MQ=255)

                                                                                      cAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATAcc   <  2:2103465/93‑1 (MQ=255)

                                                                                       aGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATACCg  <  1:467860/93‑1 (MQ=255)

TCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATACCG  >  NC_000913/644944‑645123

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