Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,368,810 A→T 5.2% S258C (AGC→TGC)  arnA → fused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,368,8100AT 5.2% 248.1 / 3.5 77S258C (AGC→TGC) arnAfused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base T (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

AATGGCATAAACCGGCATCCGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTGAATAG  >  NC_000913/2368664‑2368955
                                                                                                                                                  |                                                                                                                                                 
aaTGGCATAAACCGGCATCCGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGa                                                                                                                                                   <  6:351603/147‑1 (MQ=255)
 aTGGCATAAACCGGCATCCGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGc                                                                                                                                                 >  1:263511/1‑148 (MQ=255)
  tGGCATAATCCGGCATCCGTACTGCACAACATGGTACGGGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCg                                                                                                                                                <  3:98737/148‑1 (MQ=255)
  tGGCATAAACCGGCATCCGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCg                                                                                                                                                <  3:104767/148‑1 (MQ=255)
    gCATAAACCGGCATCCGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTg                                                                                                                                              <  4:430767/148‑1 (MQ=255)
         caCCCGCACCCGTACTGCACAACATGGTACGAGCCGCTGCCGATCCGTGGCCGGGTTCCTTCAGCTCTGTTGGCACTCAGAAATTCACCGTCTGGTCGGCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTc                                                                                                                                         <  2:384832/147‑1 (MQ=255)
              gCATCCGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTg                                                                                                                                    <  8:195304/148‑1 (MQ=255)
                   cGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCa                                                                                                                               <  1:212643/148‑1 (MQ=255)
                   cGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCa                                                                                                                               >  2:6612/1‑148 (MQ=255)
                         gCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTg                                                                                                                                       <  8:308552/134‑1 (MQ=255)
                         gCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTg                                                                                                                                       >  7:308552/1‑134 (MQ=255)
                         gCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCAc                                                                                                                              >  3:276568/1‑143 (MQ=255)
                         gCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCAc                                                                                                                              <  4:276568/143‑1 (MQ=255)
                           acaacaTGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGa                                                                                                                        <  1:171422/147‑1 (MQ=255)
                           acaacaTGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGa                                                                                                                        >  2:171422/1‑147 (MQ=255)
                           acaacaTGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGAt                                                                                                                       >  3:410632/1‑148 (MQ=255)
                           acaacaTGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGAt                                                                                                                       >  7:41684/1‑148 (MQ=255)
                            caacaTGGTACGTGCCGTTGCCGATCCGTGGCCTGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTg                                                                                                                                       <  8:327749/131‑1 (MQ=255)
                            caacaTGGTACGTGCCGTTGCCGATCCGTGGCCTGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTg                                                                                                                                       >  7:327749/1‑131 (MQ=255)
                                  ggTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTc                                                                                                                                         >  2:438217/1‑123 (MQ=255)
                                  ggTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTc                                                                                                                                         <  1:438217/123‑1 (MQ=255)
                                   gTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCtgtg                                                                                                               <  4:160147/148‑1 (MQ=255)
                                    tACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGTGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTgg                                                                                                              >  8:429809/1‑148 (MQ=255)
                                         gCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATg                                                                                                         >  4:93053/1‑148 (MQ=255)
                                              tGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCATCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGgcgcg                                                                                                    >  2:333180/1‑148 (MQ=255)
                                              tGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGgcgcg                                                                                                    >  1:447677/1‑148 (MQ=255)
                                                   aTCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTgg                                                                                                >  6:181017/1‑147 (MQ=255)
                                                    tCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCACAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGaa                                                                                              <  3:229780/148‑1 (MQ=255)
                                                    tCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGaa                                                                                              <  4:410632/148‑1 (MQ=255)
                                                    tCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGa                                                                                               >  6:354299/1‑147 (MQ=255)
                                                     ccGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGaaa                                                                                             >  8:476071/1‑148 (MQ=255)
                                                          gCCGGGCGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGCGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTTGAAATCGTc                                                                                        <  8:365254/148‑1 (MQ=255)
                                                           ccGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCa                                                                                       <  3:93053/148‑1 (MQ=255)
                                                               gTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTAATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCgg                                                                                   >  6:420295/1‑148 (MQ=255)
                                                                tGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGa                                                                                  <  2:263511/148‑1 (MQ=255)
                                                                     tCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAggc                                                                             >  2:469894/1‑148 (MQ=255)
                                                                      cAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAggcg                                                                            >  5:192479/1‑148 (MQ=255)
                                                                      cAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAggc                                                                             <  3:443131/147‑1 (MQ=255)
                                                                           aTGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGa                                                                       >  5:121651/1‑148 (MQ=255)
                                                                            tGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGTGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGAc                                                                      <  7:429809/148‑1 (MQ=255)
                                                                            tGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGa                                                                       <  1:202993/147‑1 (MQ=255)
                                                                                 gCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCt                                                                                                                                        <  8:379472/77‑1 (MQ=255)
                                                                                 gCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCt                                                                                                                                        >  7:379472/1‑77 (MQ=255)
                                                                                 gCAATCAGAAATTCACCGTCTGGTCGTCGCGTGGTCATCCTCATGCCATCAAAGCACTGCCGGGGAGCGTGATTTCTGTTGCGCCAGTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCAGCTGACAGGCGGGCGACGGCa                                                                  <  6:19509/147‑1 (MQ=255)
                                                                                 gCAATCAGAAATTCACCGTCTGGTCGTCGCGGGGTCATCCTCATGCCAGCAAAGCACACCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCAt                                                                 <  3:113772/148‑1 (MQ=255)
                                                                                      cAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGcc                                                                                                                                >  3:184183/1‑80 (MQ=255)
                                                                                      cAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGcc                                                                                                                                <  4:184183/80‑1 (MQ=255)
                                                                                      cAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTa                                                            >  1:87525/1‑148 (MQ=255)
                                                                                         aaaTTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGTGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCg                                                                                                            <  6:357113/97‑1 (MQ=255)
                                                                                         aaaTTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGTGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCg                                                                                                            >  5:357113/1‑97 (MQ=255)
                                                                                         aaaTTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGc                                                         <  1:469894/148‑1 (MQ=255)
                                                                                          aaTTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATg                                                          >  3:32992/1‑146 (MQ=255)
                                                                                           aTTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAg                                                       <  6:192479/148‑1 (MQ=255)
                                                                                                 cGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCg                                                 <  8:41684/148‑1 (MQ=255)
                                                                                                   tCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCa                                               <  2:87525/148‑1 (MQ=255)
                                                                                                                 gtTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAgcgc                                       <  7:181768/142‑1 (MQ=255)
                                                                                                                 gtTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAgcgc                                       >  8:181768/1‑142 (MQ=255)
                                                                                                                 gtTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGc                                 >  2:437839/1‑148 (MQ=255)
                                                                                                                   tCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAAtt                                            >  2:72608/1‑135 (MQ=255)
                                                                                                                   tCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAAtt                                            <  1:72608/135‑1 (MQ=255)
                                                                                                                     aTCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTggg                             <  5:420295/148‑1 (MQ=255)
                                                                                                                        cTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCt                          <  4:168017/148‑1 (MQ=255)
                                                                                                                          cATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTgg                        >  1:46665/1‑148 (MQ=255)
                                                                                                                          cATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTgg                        >  5:423447/1‑148 (MQ=255)
                                                                                                                            tGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTg                      >  5:248957/1‑148 (MQ=255)
                                                                                                                             gCCAGCAAAGCACAGCCGGGGTGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCgg                                                                                   >  7:188814/1‑86 (MQ=255)
                                                                                                                             gCCAGCAAAGCACAGCCGGGGTGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCgg                                                                                   <  8:188814/86‑1 (MQ=255)
                                                                                                                             gCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAgg                                                                              >  5:240989/1‑91 (MQ=255)
                                                                                                                             gCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAgg                                                                              <  6:240989/91‑1 (MQ=255)
                                                                                                                              ccAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCgg                                                                                   >  2:301935/1‑85 (MQ=255)
                                                                                                                              ccAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCgg                                                                                   <  1:301935/85‑1 (MQ=255)
                                                                                                                                  cAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGt                >  5:211385/1‑148 (MQ=255)
                                                                                                                                           gCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTg       <  2:46665/148‑1 (MQ=255)
                                                                                                                                             cGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTGaa     <  6:423447/148‑1 (MQ=255)
                                                                                                                                              ggggAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTGAAt    <  2:447677/148‑1 (MQ=255)
                                                                                                                                               gggAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTGaa     >  8:234479/1‑146 (MQ=255)
                                                                                                                                                ggAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTGAATag  <  6:248957/148‑1 (MQ=255)
                                                                                                                                                ggAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTGAATag  <  5:181017/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AATGGCATAAACCGGCATCCGTACTGCACAACATGGTACGTGCCGTTGCCGATCCGTGGCCGGGTGCCTTCAGCTATGTTGGCAATCAGAAATTCACCGTCTGGTCGTCGCGTGTTCATCCTCATGCCAGCAAAGCACAGCCGGGGAGCGTGATTTCTGTTGCGCCACTGCTGATTGCCTGTGGCGATGGCGCGCTGGAAATCGTCACCGGACAGGCGGGCGACGGCATTACTATGCAGGGCTCGCAATTAGCGCAGACGCTGGGCCTGGTGCAAGGTTCACGCTTGAATAG  >  NC_000913/2368664‑2368955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: