Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,476,115 A→C N227H (AAT→CAT)  yjgL → SopA‑central‑domain‑like hexapeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,476,1150AC98.8% 569.9 / ‑5.9 162N227H (AAT→CAT) yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (89/71);  total (89/73)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

TCTACAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGAATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAGGACTTAA  >  NC_000913/4476018‑4476211
                                                                                                 |                                                                                                
 cTACAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCAtt                                                                                               <  1:1735842/100‑1 (MQ=255)
  tACAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCAttt                                                                                              >  1:3854779/1‑100 (MQ=255)
  tACAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCAttt                                                                                              >  1:630842/1‑100 (MQ=255)
   aCAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTa                                                                                             <  1:2799353/100‑1 (MQ=255)
   aCAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTa                                                                                             <  1:2083298/100‑1 (MQ=255)
    cAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTaa                                                                                            >  1:2450695/1‑100 (MQ=255)
     aGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTaaa                                                                                           >  1:837541/1‑100 (MQ=255)
      gAAATTAAATAATCTACAGCAGTTAAATAATCTACTGAAGTTAAATAATATACCGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAAt                                                                                          >  2:3687899/1‑100 (MQ=38)
      gAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAAt                                                                                          >  2:827478/1‑100 (MQ=255)
         aTTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGAt                                                                                       <  2:193454/100‑1 (MQ=255)
           tAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATtc                                                                                     >  1:1598737/1‑100 (MQ=255)
             aataatCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATtctc                                                                                   >  2:3445636/1‑100 (MQ=255)
             aataatCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATtctc                                                                                   >  2:1029479/1‑100 (MQ=255)
              ataatCGACAGAAGTTAAATAAACTACTGAATGTAAATAATATACTTGTGTTAAATAATCCGCAAGAGTTAAATAATCCGCAGCATTTAAATGATTCTCa                                                                                  <  2:819937/100‑1 (MQ=255)
               taatCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAg                                                                                 <  1:2568096/100‑1 (MQ=255)
                aatCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAgg                                                                                >  2:656278/1‑100 (MQ=255)
                  tCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAg                                                                              <  2:1847892/100‑1 (MQ=255)
                  gcTACAGAAGTTAAATAATCTACTGAAGTTAAATATTATCCAGGGGTGAAATAATCCTCAGGAGTTATAGGATCCGCAGAATTTAAATGATTCTCAGGAg                                                                              <  2:243209/99‑1 (MQ=255)
                   cTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGt                                                                             <  2:2709457/100‑1 (MQ=255)
                    tACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGtt                                                                            >  1:599532/1‑100 (MQ=255)
                    tACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGtt                                                                            >  2:930578/1‑100 (MQ=255)
                      cAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTaa                                                                          >  1:3514824/1‑100 (MQ=255)
                      cAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTaa                                                                          >  2:1448083/1‑100 (MQ=255)
                      cAGAAGTTAAATAATCTACCGAAGTTAAACAATATACAGGGGTGAAAAAATCCTATGGAGATAAATAATCCGCAGCATTTAAATGATTCTTAGGTGTTaa                                                                          <  2:948800/100‑1 (MQ=25)
                        gAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTaaan                                                                        <  1:2681371/100‑2 (MQ=255)
                           gTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAAc                                                                     >  1:3974108/1‑100 (MQ=255)
                            ttAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACt                                                                    >  1:901345/1‑100 (MQ=255)
                              aaaTAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCg                                                                  <  2:3307980/100‑1 (MQ=255)
                              aaaTAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCg                                                                  >  2:1229958/1‑100 (MQ=255)
                                ataatCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCCTTTAAATGATTCTCAGGAGTTAAATAACTCgaa                                                                >  2:2516285/1‑98 (MQ=255)
                                ataatCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCa                                                                <  1:1027146/100‑1 (MQ=255)
                                ataatCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCa                                                                >  1:2514164/1‑100 (MQ=255)
                                 taatCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAg                                                               >  2:2094956/1‑100 (MQ=255)
                                  aatCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAgg                                                              <  1:2463570/100‑1 (MQ=255)
                                  aatCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAgg                                                              <  1:833625/100‑1 (MQ=255)
                                      tACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAAtt                                                          >  2:44024/1‑100 (MQ=255)
                                      tACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAAtt                                                          >  2:3927883/1‑100 (MQ=255)
                                       aCTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTa                                                         >  1:3559518/1‑100 (MQ=255)
                                        cTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTaa                                                        >  2:336346/1‑100 (MQ=255)
                                         tGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTaaa                                                       >  2:3451802/1‑100 (MQ=255)
                                          gAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAAt                                                      <  2:3674107/100‑1 (MQ=255)
                                          gAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAAt                                                      >  2:1982111/1‑100 (MQ=255)
                                          gAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAAt                                                      >  1:988201/1‑100 (MQ=255)
                                          gAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCACAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAAt                                                      >  1:3350774/1‑100 (MQ=255)
                                            aGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAg                                                    <  1:744825/100‑1 (MQ=255)
                                            aGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAg                                                    <  1:1168143/100‑1 (MQ=255)
                                            aGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAg                                                    >  1:3285568/1‑100 (MQ=255)
                                             gTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGt                                                   >  1:2317301/1‑100 (MQ=255)
                                             gTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGt                                                   >  1:413301/1‑100 (MQ=255)
                                             gTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAAGTAAATAGt                                                   >  2:269416/1‑100 (MQ=255)
                                              ttAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTc                                                  <  2:2927618/100‑1 (MQ=255)
                                               tAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTcc                                                 <  1:1376565/100‑1 (MQ=255)
                                                 aataatATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCcac                                               >  1:2933377/1‑100 (MQ=255)
                                                   taataTACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAg                                             >  2:2389303/1‑100 (MQ=255)
                                                     atatACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGa                                           >  1:1028117/1‑100 (MQ=255)
                                                      tataCAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAg                                          <  1:3555192/100‑1 (MQ=255)
                                                      tataCAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAg                                          <  2:2585819/100‑1 (MQ=255)
                                                      tataCAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAg                                          <  2:2859528/100‑1 (MQ=255)
                                                       ataCAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGt                                         <  1:425190/100‑1 (MQ=255)
                                                        taCAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGtt                                        >  1:1140550/1‑100 (MQ=255)
                                                        taCAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGtt                                        >  2:3311869/1‑100 (MQ=255)
                                                          cAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTaa                                      >  1:919627/1‑100 (MQ=255)
                                                          cAGGGGTTAAATAATCCTCAAGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTaa                                      <  1:1114592/100‑1 (MQ=255)
                                                           aGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTaaa                                     <  1:1497086/100‑1 (MQ=255)
                                                            ggggTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAAt                                    >  1:1841449/1‑100 (MQ=255)
                                                            ggggTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAAt                                    >  1:3407358/1‑100 (MQ=255)
                                                             gggTTAAATAATCCTCAGGAGTTAAATAATCCGCAGAATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATg                                   <  2:32353/100‑1 (MQ=255)
                                                             gcgTTAACTCATCCTCCGGAGTTAAATAAATCGCAGCAATTACATGATTCTCAGGAGTTAAACAACTCGCAGGAATTAAATAGTCCACAGGAGCTAAATg                                   <  2:3391840/98‑1 (MQ=255)
                                                              ggTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGa                                  >  2:1806596/1‑100 (MQ=255)
                                                              ggTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGa                                  <  2:1755821/100‑1 (MQ=255)
                                                              ggTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGa                                  <  1:225111/100‑1 (MQ=255)
                                                               gTTAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTGAATGATTCCCAGAAGTTAAATAACTCAAAGGAATTTAATAGCCCACAGGACTTAAATGaa                                 >  2:3072662/1‑99 (MQ=25)
                                                                ttAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGAtg                                >  1:4080891/1‑99 (MQ=255)
                                                                 tAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATcc                               <  2:1376361/100‑1 (MQ=255)
                                                                 tAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATcc                               <  2:870731/100‑1 (MQ=255)
                                                                 tAAATAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATcc                               <  1:1942386/100‑1 (MQ=255)
                                                                  aaaTAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCg                              <  1:1659145/100‑1 (MQ=255)
                                                                  aaaTAATCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCg                              >  2:679971/1‑100 (MQ=255)
                                                                   aataatCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGc                             >  2:906733/1‑100 (MQ=255)
                                                                      aatCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAgg                          <  1:318297/100‑1 (MQ=255)
                                                                       atCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGa                         <  1:1661337/100‑1 (MQ=255)
                                                                        tcctccGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAg                        <  2:3171241/100‑1 (MQ=255)
                                                                        tCCTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAg                        <  1:1497909/100‑1 (MQ=255)
                                                                          cTCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGtt                      >  2:483145/1‑100 (MQ=255)
                                                                           tCAGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTa                     <  2:2823938/100‑1 (MQ=255)
                                                                             aGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTaaa                   >  1:1676251/1‑100 (MQ=255)
                                                                             aGGAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTaaa                   <  2:3162514/100‑1 (MQ=255)
                                                                              ggAGTTAAATAATCCGCAGCATTTAAATGATTNTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTaaan                  <  2:3013976/100‑2 (MQ=255)
                                                                              ggAGTTAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAaat                  <  1:3664407/100‑1 (MQ=255)
                                                                                   tAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCGACAGGAGTTATATGATCCGCAGGAGTTAAATAATtc             >  2:3786117/1‑100 (MQ=255)
                                                                                   tAAATAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATtc             >  1:47216/1‑100 (MQ=255)
                                                                                    aaaTAATCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATtct            >  2:1533301/1‑100 (MQ=255)
                                                                                      ataatCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCa          <  2:910224/100‑1 (MQ=255)
                                                                                       taatCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAg         <  2:1558266/100‑1 (MQ=255)
                                                                                       taatCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAg         <  2:1808894/100‑1 (MQ=255)
                                                                                        aatCCGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAgg        <  1:3229007/100‑1 (MQ=255)
                                                                                            cGCAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAGGACtt    >  2:1172638/1‑100 (MQ=255)
                                                                                             gcagcaTTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAGGACTTa   >  2:536405/1‑100 (MQ=255)
                                                                                             gcagcaTTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAGGACTTa   <  2:3527695/100‑1 (MQ=255)
                                                                                              cAGCATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAGGACTTaa  >  2:3217202/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TCTACAGAAATTAAATAATCTACAGAAGTTAAATAATCTACTGAAGTTAAATAATATACAGGGGTTAAATAATCCTCAGGAGTTAAATAATCCGCAGAATTTAAATGATTCTCAGGAGTTAAATAACTCGCAGGAATTAAATAGTCCACAGGAGTTAAATGATCCGCAGGAGTTAAATAATTCTCAGGACTTAA  >  NC_000913/4476018‑4476211

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