BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,721,938	A→G	V598V (GTA→GTG)	b1645 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,721,938	0	A	G	99.2%	410.3 / ‑4.6	133	V598V (GTA→GTG)	b1645	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (62/70); total (62/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 135 total aligned reads displayed.

TTACGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACA  >  NC_000913/1721849‑1722027

ttaCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCgg                                                                                        >  1:2278095/1‑93 (MQ=255)

ttaCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCgg                                                                                        >  1:2850179/1‑93 (MQ=255)

  aCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGa                                                                                      <  2:4364314/93‑1 (MQ=255)

  aCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGa                                                                                      >  2:66887/1‑93 (MQ=255)

     cTGGGGTGTAGTGCTGCTTAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATaa                                                                                   <  1:924089/93‑1 (MQ=255)

     cTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATaa                                                                                   <  2:2278095/93‑1 (MQ=255)

     cTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATaa                                                                                   <  2:3046287/93‑1 (MQ=255)

     cTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATaa                                                                                   <  2:2997823/93‑1 (MQ=255)

      tGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAAc                                                                                  <  1:1409220/93‑1 (MQ=255)

      gggggTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAAc                                                                                  <  2:858939/92‑1 (MQ=255)

       ggggTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACt                                                                                 <  2:2937465/93‑1 (MQ=255)

        gggTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTg                                                                                <  2:2419072/93‑1 (MQ=255)

         ggggtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGATCCGATCCGTTATCGCGAGTGCGGGATAACTGt                                                                               >  2:3421200/4‑93 (MQ=255)

          gtgtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTa                                                                              <  1:1565594/93‑1 (MQ=255)

          gtgtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTa                                                                              <  2:3963937/93‑1 (MQ=255)

            gtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAtt                                                                            <  2:1002613/93‑1 (MQ=255)

            gtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAtt                                                                            <  2:2365313/93‑1 (MQ=255)

             tAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAttt                                                                           >  2:1085700/1‑93 (MQ=255)

             tAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAttt                                                                           <  1:4556713/93‑1 (MQ=255)

              aGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTc                                                                          <  2:2567860/93‑1 (MQ=255)

              aGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTc                                                                          <  1:2820323/93‑1 (MQ=255)

              aGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTc                                                                          <  2:554317/93‑1 (MQ=255)

               gTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCa                                                                         <  2:2434163/93‑1 (MQ=255)

               gTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCa                                                                         <  1:38756/93‑1 (MQ=255)

                tgctgcTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCAc                                                                        >  1:4485985/1‑93 (MQ=255)

                tgctgcTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCAc                                                                        >  1:4485984/1‑93 (MQ=255)

                  ctgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTg                                                                      <  1:1661811/93‑1 (MQ=255)

                  ctgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTg                                                                      <  2:1852004/93‑1 (MQ=255)

                  ctgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTg                                                                      <  2:3361021/93‑1 (MQ=255)

                   tgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGt                                                                     >  1:243625/1‑93 (MQ=255)

                   tgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGt                                                                     >  1:743057/1‑93 (MQ=255)

                      tgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGc                                                                  >  2:2244452/1‑93 (MQ=255)

                      tgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGc                                                                  <  2:434927/93‑1 (MQ=255)

                       gAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCg                                                                 <  2:1642975/93‑1 (MQ=255)

                       gAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCg                                                                 <  2:1273428/93‑1 (MQ=255)

                       gAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCg                                                                 >  1:2138807/1‑93 (MQ=255)

                         aCTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCggg                                                               <  2:2710110/93‑1 (MQ=255)

                          cTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCgggg                                                              >  2:2587147/1‑93 (MQ=255)

                          cTGTCCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCgggg                                                              <  1:4308483/93‑1 (MQ=255)

                           tGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGa                                                             <  1:2922254/93‑1 (MQ=255)

                            gTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAg                                                            >  1:804889/1‑93 (MQ=255)

                            gTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAg                                                            <  1:595399/93‑1 (MQ=255)

                               ctcATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtga                                                         <  2:54338/93‑1 (MQ=255)

                               ctcATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtga                                                         <  1:3133982/93‑1 (MQ=255)

                                  aTGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtgatga                                                      >  1:3094177/1‑93 (MQ=255)

                                  aTGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtgatga                                                      >  1:4259947/1‑93 (MQ=255)

                                      tgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtga                                                  >  1:2347194/1‑93 (MQ=255)

                                      tgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCAAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtga                                                  <  1:3043150/93‑1 (MQ=255)

                                       gtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtgag                                                 <  2:2876828/93‑1 (MQ=255)

                                        tCTGGCAATTGCGCGGCTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCCCTGTTGCGGGGAGTGATGAGTGAGc                                                <  1:3005143/93‑1 (MQ=255)

                                          tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              >  1:99784/1‑93 (MQ=255)

                                          tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              <  2:3960193/93‑1 (MQ=255)

                                          tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              <  1:3834552/93‑1 (MQ=255)

                                          tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              <  2:906982/93‑1 (MQ=255)

                                           ggCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTg                                             <  1:1726935/93‑1 (MQ=255)

                                             cAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGc                                           <  2:2142233/93‑1 (MQ=255)

                                              aaTTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCg                                          >  2:1322058/1‑93 (MQ=255)

                                               aTTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGt                                         >  1:2254164/1‑93 (MQ=255)

                                                ttGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGtt                                        >  2:4241496/1‑93 (MQ=255)

                                                ttGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGtt                                        >  1:2246737/1‑93 (MQ=255)

                                                    gcgACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGc                                    <  2:3747624/93‑1 (MQ=255)

                                                      gACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaa                                  >  1:2281888/1‑93 (MQ=255)

                                                       aCTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaa                                 >  1:1851787/1‑93 (MQ=255)

                                                         tGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaaaa                               <  1:2596570/93‑1 (MQ=255)

                                                          gggAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAAt                              <  2:3468947/93‑1 (MQ=255)

                                                            gAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGATTGAGCGTGGCGTTCAGCAAAAATCa                            <  2:919533/93‑1 (MQ=255)

                                                            gAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCa                            <  2:4428188/93‑1 (MQ=255)

                                                             aaTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCAc                           <  2:487359/93‑1 (MQ=255)

                                                              aTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                          >  2:376891/1‑93 (MQ=255)

                                                              aTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                          <  2:358071/93‑1 (MQ=255)

                                                              aTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                          >  2:3239443/1‑93 (MQ=255)

                                                              aTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                          >  1:4127836/1‑93 (MQ=255)

                                                               tCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTg                         >  1:1388882/1‑93 (MQ=255)

                                                               tCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTg                         >  1:1591526/1‑93 (MQ=255)

                                                               tCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTg                         >  2:4476939/1‑93 (MQ=255)

                                                               tCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATAACTg                         >  2:4476940/1‑93 (MQ=255)

                                                               tCGCGTTCCGATCCGATATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTg                         <  1:3737638/93‑1 (MQ=255)

                                                                   gTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcgg                     >  1:562079/1‑93 (MQ=255)

                                                                    ttCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcggc                    >  2:1682737/1‑93 (MQ=255)

                                                                    ttCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcggc                    >  2:1100946/1‑93 (MQ=255)

                                                                    ctCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTagccgc                    >  1:1497621/2‑93 (MQ=255)

                                                                       cGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCcac                 <  1:2602452/93‑1 (MQ=255)

                                                                       cGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCcac                 >  1:2024733/1‑93 (MQ=255)

                                                                         aTCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCacac               >  2:826349/1‑93 (MQ=255)

                                                                          tCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACt              >  1:1580081/1‑93 (MQ=255)

                                                                          tCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACt              >  1:2158504/1‑93 (MQ=255)

                                                                             gTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTga           >  1:3466977/1‑93 (MQ=255)

                                                                             gTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTga           <  2:1659543/93‑1 (MQ=255)

                                                                              ttATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaa          <  2:143865/93‑1 (MQ=255)

                                                                               tATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaag         <  2:1090173/93‑1 (MQ=255)

                                                                               tATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaag         <  2:3941820/93‑1 (MQ=255)

                                                                               tATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaag         <  2:1521945/93‑1 (MQ=255)

                                                                                 tCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaagaa       >  1:2271831/1‑93 (MQ=255)

                                                                                  cgcgAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAAt      <  1:2726262/93‑1 (MQ=255)

                                                                                  cgcgAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAAt      <  2:3198320/93‑1 (MQ=255)

                                                                                   gcgAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAAtt     >  2:126514/1‑93 (MQ=255)

                                                                                     gAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTAc   >  2:2036541/1‑93 (MQ=255)

                                                                                     gAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTAc   <  1:84074/93‑1 (MQ=255)

                                                                                      aGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACa  <  1:2486910/93‑1 (MQ=255)

                                                                                      aGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACa  <  1:3212922/93‑1 (MQ=255)

TTACGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACA  >  NC_000913/1721849‑1722027

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