BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,564,047	T→C	Y373C (TAC→TGC)	b1490 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,564,047	0	T	C	87.9%	692.8 / 77.5	273	Y373C (TAC→TGC)	b1490	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (13/20); new base C (119/121); total (132/141)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 273 total aligned reads displayed.

TCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCAC  >  NC_000913/1563958‑1564136

tctcCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGc                                                                                        <  2:5820871/93‑1 (MQ=255)

tctcCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGc                                                                                        >  2:170864/1‑93 (MQ=255)

 ctcCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCg                                                                                       >  1:4203023/1‑93 (MQ=255)

  tcCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCgg                                                                                      <  1:3267686/93‑1 (MQ=255)

   ccACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGa                                                                                     >  1:2513833/1‑93 (MQ=255)

    cACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGaa                                                                                    >  1:6361408/1‑93 (MQ=255)

     aCCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGaaa                                                                                   >  2:4999651/1‑93 (MQ=255)

     aCCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGaaa                                                                                   <  1:5351395/93‑1 (MQ=255)

     aCCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGaaa                                                                                   >  1:2254390/1‑93 (MQ=255)

      cccGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGaaaa                                                                                  >  2:6654468/1‑93 (MQ=255)

       ccGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGaaaaa                                                                                 >  2:1107250/1‑92 (MQ=255)

       ccGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAAc                                                                                 <  2:891237/93‑1 (MQ=255)

       ccGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAAc                                                                                 <  1:5635703/93‑1 (MQ=255)

        cGACTGCGAATACGTCCTGCGGTACGTACCGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACa                                                                                <  1:3656887/93‑1 (MQ=255)

          aCTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATa                                                                              <  2:6746939/93‑1 (MQ=255)

          aCTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATa                                                                              <  1:4941091/93‑1 (MQ=255)

          aCTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATa                                                                              <  1:6789025/93‑1 (MQ=255)

           cTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATaa                                                                             <  2:7028825/93‑1 (MQ=255)

           cTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATaa                                                                             <  1:1186484/93‑1 (MQ=255)

            tGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAAt                                                                            >  2:5370659/1‑93 (MQ=255)

            tGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAAt                                                                            <  2:6285208/93‑1 (MQ=255)

            tGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAAt                                                                            >  1:3909261/1‑93 (MQ=255)

             gCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATc                                                                           >  1:7097804/1‑93 (MQ=255)

              cGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCa                                                                          <  2:216856/93‑1 (MQ=255)

              cGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCa                                                                          <  1:2294050/93‑1 (MQ=255)

               gAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCac                                                                         <  1:1996505/93‑1 (MQ=255)

                aaTACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCTCCCCCGCAGCGGAAAACATAATCact                                                                        >  2:4148122/1‑93 (MQ=255)

                aaTACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCact                                                                        <  2:3693184/93‑1 (MQ=255)

                aaTACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCact                                                                        >  2:6557746/1‑93 (MQ=255)

                  tACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCactac                                                                      <  1:7295040/93‑1 (MQ=255)

                  tACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCactac                                                                      >  1:5082653/1‑93 (MQ=255)

                  tACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCactac                                                                      >  2:5691754/1‑93 (MQ=255)

                      ttCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCg                                                                  >  1:2929546/1‑93 (MQ=255)

                           cGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGt                                                             <  1:4381566/93‑1 (MQ=255)

                             gTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTg                                                           <  2:2201042/93‑1 (MQ=255)

                             gTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTg                                                           >  2:594637/1‑93 (MQ=255)

                             gTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTg                                                           >  2:1509970/1‑93 (MQ=255)

                             gTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTg                                                           <  2:4563732/93‑1 (MQ=255)

                              tACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGtttt                                                          >  1:2419853/1‑91 (MQ=255)

                               aCGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTc                                                         <  2:3279160/93‑1 (MQ=255)

                               aCGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTc                                                         >  2:1872984/1‑93 (MQ=255)

                               aCGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTc                                                         >  1:2763714/1‑93 (MQ=255)

                                cGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCa                                                        >  2:7471587/1‑93 (MQ=255)

                                 gTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCAt                                                       >  1:6492670/1‑93 (MQ=255)

                                 gTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCAt                                                       >  2:2461750/1‑93 (MQ=255)

                                  tAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATa                                                      >  1:2138198/1‑93 (MQ=255)

                                  tAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATa                                                      >  2:7351412/1‑93 (MQ=255)

                                   aaCGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATaa                                                     <  1:14456/93‑1 (MQ=255)

                                    aCGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATaaa                                                    <  1:7471587/93‑1 (MQ=255)

                                     cgtttcgtttTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATaaaa                                                   <  1:2988957/93‑1 (MQ=255)

                                       tttcgtttTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAgg                                                 <  2:2618125/93‑1 (MQ=255)

                                       tttcgtttTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAgg                                                 >  2:7162442/1‑93 (MQ=255)

                                         tcgtttTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGcc                                               >  1:6883005/1‑93 (MQ=255)

                                          cgtttTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCt                                              <  1:5761161/93‑1 (MQ=255)

                                          cgtttTCAGAAGCTTCAATCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCt                                              >  2:4903436/1‑93 (MQ=255)

                                            ttttCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTga                                            <  2:3804480/93‑1 (MQ=255)

                                            ttttCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTga                                            <  2:3970752/93‑1 (MQ=255)

                                            ttttCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTga                                            <  1:6717672/93‑1 (MQ=255)

                                            ttttCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTga                                            <  1:1437062/93‑1 (MQ=255)

                                              ttCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTgaga                                          <  1:3195308/93‑1 (MQ=255)

                                                 aGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGAc                                       >  1:6073228/1‑93 (MQ=255)

                                                 aGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGAc                                       >  2:4893681/1‑93 (MQ=255)

                                                 aGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGAc                                       >  1:6988648/1‑93 (MQ=255)

                                                  gAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACt                                      >  1:6822832/1‑93 (MQ=255)

                                                  gAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACt                                      <  1:1288771/93‑1 (MQ=255)

                                                   aaGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACtt                                     <  2:5165960/93‑1 (MQ=255)

                                                     gCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTa                                   <  2:1636/93‑1 (MQ=255)

                                                     gCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTa                                   <  1:1889102/93‑1 (MQ=255)

                                                      cTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTAc                                  <  2:4843523/93‑1 (MQ=255)

                                                      cTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTAc                                  >  1:4105037/1‑93 (MQ=255)

                                                      cTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTAc                                  <  1:1733090/93‑1 (MQ=255)

                                                      cTTCAGTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTAc                                  >  1:5569561/1‑93 (MQ=255)

                                                           gTCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGa                             >  2:3589030/1‑93 (MQ=255)

                                                            tCAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGaa                            <  2:5048263/93‑1 (MQ=255)

                                                             cAAAACAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAAt                           >  1:2539214/1‑93 (MQ=255)

                                                              aaaaCAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGATACTTTACGCAGAAtt                          >  2:591294/1‑93 (MQ=255)

                                                               aaaCAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAAttt                         >  2:3056041/1‑93 (MQ=255)

                                                               aaaCAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAAttn                         >  2:2431730/1‑92 (MQ=255)

                                                                 aCAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTtca                       >  1:4009225/1‑93 (MQ=255)

                                                                 aCAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTtca                       >  1:1546780/1‑93 (MQ=255)

                                                                  cAATGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTtcat                      <  1:5264597/93‑1 (MQ=255)

                                                                   aaTGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTtcatc                     >  2:658688/1‑93 (MQ=255)

                                                                   aaTGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTtcatc                     <  1:3283425/93‑1 (MQ=255)

                                                                    aTGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTtcatca                    >  1:2740258/1‑93 (MQ=255)

                                                                     tGATAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCAc                   >  2:811324/1‑93 (MQ=255)

                                                                       aTAAATTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCa                 <  2:2950701/93‑1 (MQ=255)

                                                                         aaaTTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAgt               <  2:1308586/93‑1 (MQ=255)

                                                                         aaaTTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAgt               <  1:106751/93‑1 (MQ=255)

                                                                         aaaTTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAgt               <  1:5517959/93‑1 (MQ=255)

                                                                          aaTTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAgtg              >  2:3433328/1‑93 (MQ=255)

                                                                          aaTTCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAgtg              <  2:3129983/93‑1 (MQ=255)

                                                                            ttCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGtt            >  2:179384/1‑93 (MQ=255)

                                                                             tCATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTa           <  2:2744936/93‑1 (MQ=255)

                                                                              cATCGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTAt          <  1:3363964/93‑1 (MQ=255)

                                                                               aTCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATg         <  1:1832393/93‑1 (MQ=255)

                                                                                 cGCCCCCGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGc       <  1:4610138/93‑1 (MQ=255)

                                                                                    ccccGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGcccc    >  1:6574697/1‑93 (MQ=255)

                                                                                     cccGCAGCGGAAAACATAATCACTACGGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCa   >  2:2297876/1‑93 (MQ=255)

                                                                                      ccGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCAc  <  2:6352474/93‑1 (MQ=255)

                                                                                      ccGCAGCGGAAAACATAATCACTACTGCGGACGTTGTCACAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCAc  >  1:6214150/1‑93 (MQ=255)

TCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCAC  >  NC_000913/1563958‑1564136

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