Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,291,084 A→C 27.4% K195N (AAA→AAC narP → response regulator in two‑component regulatory system with NarQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,291,0840AC27.4% 188.0 / 109.9 128K195N (AAA→AACnarPresponse regulator in two‑component regulatory system with NarQ
Reads supporting (aligned to +/- strand):  ref base A (51/42);  new base C (21/14);  total (72/56)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-02

Only 100 of 129 total aligned reads displayed.

CACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTT  >  NC_000913/2290986‑2291183
                                                                                                  |                                                                                                   
cACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAAc                                                                                                    <  1:2194330/100‑1 (MQ=255)
   gAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCa                                                                                                 <  2:1280249/100‑1 (MQ=255)
     gCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAAt                                                                                               >  2:1271546/1‑100 (MQ=255)
      cTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATg                                                                                              >  2:2077034/1‑100 (MQ=255)
       tGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGt                                                                                             >  2:6550/1‑100 (MQ=255)
         gCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTcc                                                                                           <  2:2622082/100‑1 (MQ=255)
          cacaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCg                                                                                          <  2:2407743/100‑1 (MQ=255)
            caGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCt                                                                                        <  1:2158746/100‑1 (MQ=255)
             aGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTc                                                                                       <  2:599614/100‑1 (MQ=255)
              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCa                                                                                      >  1:2150229/1‑100 (MQ=255)
              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCa                                                                                      >  2:3183654/1‑100 (MQ=255)
               gggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAc                                                                                     >  2:3261856/1‑100 (MQ=255)
                 gCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcgc                                                                                   >  1:1445231/1‑100 (MQ=255)
                  ttccgatctAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcgcg                                                                                  <  2:1427656/92‑1 (MQ=255)
                   tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGt                                                                                 >  2:1820160/1‑100 (MQ=255)
                   tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGt                                                                                 >  2:2206686/1‑100 (MQ=255)
                    gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCGCACGCGTg                                                                                >  2:2256249/1‑100 (MQ=255)
                     tCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTgg                                                                               >  2:1912571/1‑100 (MQ=255)
                     tCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTgg                                                                               >  1:175535/1‑100 (MQ=255)
                     tCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTgg                                                                               <  1:1209475/100‑1 (MQ=255)
                      cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGAAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                                                                              >  1:721615/1‑100 (MQ=255)
                       aaaTAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcg                                                                             >  1:3172497/1‑100 (MQ=255)
                        aaTAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcgg                                                                            <  2:1445237/100‑1 (MQ=255)
                          tAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGcc                                                                          >  2:716770/1‑100 (MQ=255)
                          tAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGcc                                                                          >  1:338618/1‑100 (MQ=255)
                          tAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGcc                                                                          <  2:2666182/100‑1 (MQ=255)
                           aaaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGcca                                                                         <  1:181230/100‑1 (MQ=255)
                           aaaCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcgcgagatcggaa                                                                         >  1:1427650/1‑91 (MQ=255)
                           aaaCAGATTGCCTCGGTGGTGAATATTTCCGAGAAAACAGTAAACGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcgc                                                                                   >  1:3195208/1‑90 (MQ=255)
                            agcAGAGTGCCTCGGTGTTGAATCTTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcc                                                                          <  1:1865182/96‑1 (MQ=255)
                            aaCGGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTGGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcca                                                                         >  1:1167915/1‑99 (MQ=255)
                             aCAGATTGCCTCGGTGTTGAATCTTTCCGAGCAGACAGTAAAAGTACCTATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacc                                                                       <  2:153640/100‑1 (MQ=255)
                             aCAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacc                                                                       >  2:3158178/1‑100 (MQ=255)
                              cAGATTGCCTCGGTGTTGAATATTTCCGCGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacca                                                                      >  2:398903/1‑100 (MQ=255)
                                gATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCAtt                                                                    >  2:2913667/1‑100 (MQ=255)
                                 aTTGCCTGGGTGTTGAATATGTCCGAGCAGACAGTAAAAGTACATATTCGGAATCGGCTGCGTAAACTCAATGTCCGCGCGCGCGGGGCGGACACCATTc                                                                   <  1:2358165/100‑1 (MQ=255)
                                  ttGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCt                                                                  <  2:1170111/100‑1 (MQ=255)
                                    gCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGt                                                                <  2:2573832/100‑1 (MQ=255)
                                       tCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc                                                             <  1:849777/100‑1 (MQ=255)
                                       tCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc                                                             <  2:2324104/100‑1 (MQ=255)
                                       tCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc                                                             <  2:123231/100‑1 (MQ=255)
                                         ggTGTTGATTATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGTGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTg                                                           <  1:518468/100‑1 (MQ=255)
                                         ggTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTg                                                           <  2:1077266/100‑1 (MQ=255)
                                           tgtTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGca                                                         <  1:2444620/100‑1 (MQ=255)
                                           tgtTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGca                                                         <  2:2273864/100‑1 (MQ=255)
                                             ttGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaac                                                       >  2:1491927/1‑100 (MQ=255)
                                             ttGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaac                                                       <  2:1009319/100‑1 (MQ=255)
                                             tggAATATTTCCGAGCAGACAGTAAATGTACATATTCGCAATCTGCTGCGTAAACTCCATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaac                                                       <  1:1860597/98‑1 (MQ=255)
                                              tGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca                                                      >  1:3222638/1‑100 (MQ=255)
                                              tGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca                                                      >  1:438493/1‑100 (MQ=255)
                                              tGAATATTTCCGAGCAAACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca                                                      >  2:2550867/1‑100 (MQ=255)
                                               gAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaacaa                                                     <  1:2565693/100‑1 (MQ=255)
                                                 atatTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAAcg                                                   >  2:3073896/1‑100 (MQ=255)
                                                  tatTTCCGAGCAGACAGTAAAAGTACATATTCACAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAAcgc                                                  >  2:102782/1‑100 (MQ=255)
                                                    tttCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCgg                                                >  2:1877650/1‑100 (MQ=255)
                                                    tttCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCgg                                                >  1:1147333/1‑100 (MQ=255)
                                                      tCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCgggg                                              <  2:473119/100‑1 (MQ=255)
                                                      tCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCgggg                                              <  2:2949922/100‑1 (MQ=255)
                                                       ccGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGc                                             <  2:1031509/100‑1 (MQ=255)
                                                         gAGCAGACAGTAAAAGTACATATTCGCAATCTGCTTCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGcac                                           <  2:690789/100‑1 (MQ=255)
                                                         gAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGcac                                           >  1:831804/1‑100 (MQ=255)
                                                             aGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATa                                       >  1:3328480/1‑100 (MQ=255)
                                                              gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaa                                      >  1:2420316/1‑100 (MQ=255)
                                                              gACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaa                                      >  1:2500819/1‑100 (MQ=255)
                                                               aCAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCGCAATaaa                                     <  2:1775059/100‑1 (MQ=255)
                                                               aCAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCCCTCACGCGTGGCGGCCACCATTCTGTCCCTGCAACAACGCGGGGCACAATaaa                                     <  2:2212407/100‑1 (MQ=255)
                                                                cAGTAAAAGTACATATTCGCCATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaaaa                                    >  1:2512648/1‑100 (MQ=255)
                                                                cAGTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaaaa                                    >  2:1404821/1‑100 (MQ=255)
                                                                 aGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATaaaaa                                   >  2:1134063/1‑100 (MQ=255)
                                                                  gTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                                  >  2:2109874/1‑100 (MQ=255)
                                                                  gTAAAAGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                                  >  1:303637/1‑100 (MQ=255)
                                                                  gTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                                  >  1:2967107/1‑100 (MQ=255)
                                                                  gTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                                  >  1:2814837/1‑100 (MQ=255)
                                                                  gTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAAt                                  >  2:1645588/1‑100 (MQ=255)
                                                                  gTAAAAGTACATATGCGGAATCTGGTGCGTAAACTCAATGACCGCTCAGACGTGGCGGCAACACTGCTGTTAATGCAACAACGCGGGGGAGAATaaaaa                                   >  1:1281807/1‑99 (MQ=255)
                                                                     aaaGTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGc                               <  1:336737/100‑1 (MQ=255)
                                                                      aaGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGcc                              <  1:2256233/100‑1 (MQ=255)
                                                                        gTACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTTTTCCTACAACAACGCGGGGCACAATAAAAATAGCCCg                            <  1:1813119/100‑1 (MQ=255)
                                                                        gTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCg                            <  2:2249476/100‑1 (MQ=255)
                                                                        gTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCg                            >  2:1761206/1‑100 (MQ=255)
                                                                         tACATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGa                           >  2:774696/1‑100 (MQ=255)
                                                                           cATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATg                         >  2:3039815/1‑100 (MQ=255)
                                                                           cATATTCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATg                         <  2:338619/100‑1 (MQ=255)
                                                                             tatTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGa                       <  2:120471/100‑1 (MQ=255)
                                                                               ttCGCAATCTGCTGCGTAACCTCAATGTCCGCTCACGCTTGGCGGCCACCATTCTGTTTCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAtt                     <  1:1846120/100‑1 (MQ=255)
                                                                               ttCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAtt                     <  2:3256673/100‑1 (MQ=255)
                                                                                tCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAttt                    >  1:696074/1‑100 (MQ=255)
                                                                                 cGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTa                   <  2:427088/100‑1 (MQ=255)
                                                                                   cAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTatc                 <  1:469000/100‑1 (MQ=255)
                                                                                       ctgctgCGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCg             >  1:100668/1‑100 (MQ=255)
                                                                                         gctgcGTAACCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCggg           >  1:1503986/1‑100 (MQ=255)
                                                                                         gctgcGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCggg           <  2:1605248/100‑1 (MQ=255)
                                                                                         gctgcGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCggg           >  1:3181070/1‑100 (MQ=255)
                                                                                          ctgcGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGc          <  2:457901/100‑1 (MQ=255)
                                                                                          ctgcGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGc          >  2:2601485/1‑100 (MQ=255)
                                                                                           tgcGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCt         >  2:547198/1‑100 (MQ=255)
                                                                                             cGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTga       <  2:2897305/100‑1 (MQ=255)
                                                                                              gTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTgag      >  2:177977/1‑100 (MQ=255)
                                                                                               tAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTgaga     >  1:2506656/1‑100 (MQ=255)
                                                                                                  aCTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGAttt  >  2:189240/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTT  >  NC_000913/2290986‑2291183

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