Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,289,064 A→G E182G (GAG→GGG)  narP → nitrate/nitrite response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,289,0640AG88.3% 702.4 / 76.5 274E182G (GAG→GGG) narPnitrate/nitrite response regulator
Reads supporting (aligned to +/- strand):  ref base A (3/29);  new base G (110/132);  total (113/161)
Fisher's exact test for biased strand distribution p-value = 4.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 276 total aligned reads displayed.

CAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCG  >  NC_000913/2288972‑2289156
                                                                                            |                                                                                            
cAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGa                                                                                              <  1:1781158/93‑1 (MQ=255)
  gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGc                                                                                            >  2:5592254/1‑93 (MQ=255)
     tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGa                                                                                         >  1:1562709/1‑93 (MQ=255)
     tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGa                                                                                         >  1:7076416/1‑93 (MQ=255)
     tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGa                                                                                         >  1:4069323/1‑93 (MQ=255)
     tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGa                                                                                         <  2:3150138/93‑1 (MQ=255)
        tGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAg                                                                                      <  2:141465/93‑1 (MQ=255)
        tGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAg                                                                                      <  2:711264/93‑1 (MQ=255)
         gACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGt                                                                                     <  2:6865886/93‑1 (MQ=255)
            ggAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTaaa                                                                                  <  2:807636/93‑1 (MQ=255)
             gAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTcaaa                                                                                 <  1:6975642/93‑1 (MQ=255)
             gAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTaaaa                                                                                 >  2:6606034/1‑93 (MQ=255)
             gAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTaaaa                                                                                 <  2:1258817/93‑1 (MQ=255)
              aGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAg                                                                                <  2:492499/93‑1 (MQ=255)
               gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGt                                                                               <  1:2052922/93‑1 (MQ=255)
                  cgAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACa                                                                            <  2:3978953/93‑1 (MQ=255)
                  cgAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACa                                                                            <  2:2791085/93‑1 (MQ=255)
                  cgAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACa                                                                            <  2:6448276/93‑1 (MQ=255)
                   gAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACat                                                                           <  1:327506/93‑1 (MQ=255)
                   gAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACat                                                                           <  1:196128/93‑1 (MQ=255)
                   gAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACat                                                                           <  1:5584699/93‑1 (MQ=255)
                    aGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACata                                                                          <  2:5424984/93‑1 (MQ=255)
                     gCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACatat                                                                         <  1:7247404/93‑1 (MQ=255)
                      cTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATAtt                                                                        <  1:1162004/93‑1 (MQ=255)
                       tGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTc                                                                       >  1:4647488/1‑93 (MQ=255)
                       tGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTc                                                                       >  1:934444/1‑93 (MQ=255)
                       tGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTc                                                                       >  2:2336192/1‑93 (MQ=255)
                       nngATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTc                                                                       <  1:87107/91‑1 (MQ=255)
                           tGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCaa                                                                   <  2:2026788/93‑1 (MQ=255)
                           tGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCaa                                                                   >  2:6228798/1‑93 (MQ=255)
                             ttCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATc                                                                 <  1:2959421/93‑1 (MQ=255)
                             ttCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATc                                                                 >  2:4911148/1‑93 (MQ=255)
                              tCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATct                                                                <  2:6925949/93‑1 (MQ=255)
                              tCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATct                                                                >  1:5917721/1‑93 (MQ=255)
                               cTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATctg                                                               <  1:6458177/93‑1 (MQ=255)
                                tGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATctgc                                                              <  2:3773117/93‑1 (MQ=255)
                                tGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATctgc                                                              <  2:4578554/93‑1 (MQ=255)
                                 gCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATctgct                                                             <  1:3440249/93‑1 (MQ=255)
                                 gCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATctgct                                                             >  1:4384095/1‑93 (MQ=255)
                                  cACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATctgctg                                                            <  1:3985756/93‑1 (MQ=255)
                                      aGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTa                                                        >  2:4756151/1‑93 (MQ=255)
                                       gCTGGCACAGGGGCTGTCAAATAAACAGATCGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                                       >  1:364911/1‑93 (MQ=255)
                                        cTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaaa                                                      <  2:3388504/93‑1 (MQ=255)
                                         tGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAAc                                                     <  1:5723008/93‑1 (MQ=255)
                                         tGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAAc                                                     <  1:3743320/93‑1 (MQ=255)
                                          ggCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACt                                                    <  2:2091691/93‑1 (MQ=255)
                                             acaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                                 <  2:4057622/93‑1 (MQ=255)
                                             acaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                                 <  2:2915561/93‑1 (MQ=255)
                                              caGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAAt                                                >  1:4496890/1‑93 (MQ=255)
                                               aGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATg                                               >  2:4291453/1‑93 (MQ=255)
                                                ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGt                                              <  1:4838513/93‑1 (MQ=255)
                                                 gggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTc                                             <  1:679704/93‑1 (MQ=255)
                                                  ggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTcc                                            <  2:7106830/93‑1 (MQ=255)
                                                    cTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGc                                          <  2:4261962/93‑1 (MQ=255)
                                                     tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCt                                         <  1:4543585/93‑1 (MQ=255)
                                                      gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAAt                                                                  >  2:2700306/1‑67 (MQ=255)
                                                      gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTc                                        >  1:996225/1‑93 (MQ=255)
                                                      gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTc                                        >  1:7148461/1‑93 (MQ=255)
                                                      gTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTc                                        >  1:2245154/1‑93 (MQ=255)
                                                       tCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCa                                       >  2:5727375/1‑93 (MQ=255)
                                                        cAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAc                                      <  2:5982586/93‑1 (MQ=255)
                                                         aaaTAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcg                                     >  1:6348743/1‑93 (MQ=255)
                                                          aaTAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTc                                             <  1:2822285/84‑1 (MQ=255)
                                                          aaTAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCAcgc                                    <  1:7180064/93‑1 (MQ=255)
                                                            tAAACAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGt                                  <  1:6750202/93‑1 (MQ=255)
                                                              aaCAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTgg                                >  2:6381868/1‑93 (MQ=255)
                                                              aaCAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTgg                                <  2:898904/93‑1 (MQ=255)
                                                               aCAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                               >  1:5639275/1‑93 (MQ=255)
                                                                cAGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcg                              >  2:5268353/1‑93 (MQ=255)
                                                                 aGATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcgg                             >  2:5545007/1‑93 (MQ=255)
                                                                  gATTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcggc                            <  1:5368236/93‑1 (MQ=255)
                                                                   aTTGCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcc                           >  2:2492426/1‑93 (MQ=255)
                                                                      gCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacc                        >  1:5811737/1‑93 (MQ=255)
                                                                      gCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacc                        >  1:249038/1‑93 (MQ=255)
                                                                      gCCTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacc                        >  1:2648638/1‑93 (MQ=255)
                                                                       ccTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacca                       <  1:4905664/93‑1 (MQ=255)
                                                                        cTCGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCAt                      >  1:807385/1‑93 (MQ=255)
                                                                          cGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTc                    <  2:7348845/93‑1 (MQ=255)
                                                                          cGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTc                    <  2:2476489/93‑1 (MQ=255)
                                                                          cGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTc                    <  1:5438108/93‑1 (MQ=255)
                                                                          cGGTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTc                    <  2:5886633/93‑1 (MQ=255)
                                                                           ggTGTTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCt                   <  2:5098857/93‑1 (MQ=255)
                                                                            gtgtTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTcnn                  >  1:135544/1‑91 (MQ=255)
                                                                              gtTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGtt                <  2:5035420/93‑1 (MQ=255)
                                                                              gtTGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGtt                >  1:2310735/1‑93 (MQ=255)
                                                                                tGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc              <  2:3504570/93‑1 (MQ=255)
                                                                                tGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc              <  2:1591708/93‑1 (MQ=255)
                                                                                tGAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTcc              >  1:1654736/1‑93 (MQ=255)
                                                                                 gAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCt             >  1:5882937/1‑93 (MQ=255)
                                                                                 gAATATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCt             >  2:1827277/1‑93 (MQ=255)
                                                                                  aaTATTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTg            >  2:3251415/1‑93 (MQ=255)
                                                                                     atTTCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaa         >  1:513911/1‑93 (MQ=255)
                                                                                     atTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaa         <  1:1099627/93‑1 (MQ=255)
                                                                                      tttCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaac        >  1:6686071/1‑93 (MQ=255)
                                                                                      tttCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaac        <  2:2173728/93‑1 (MQ=255)
                                                                                      tttCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaac        <  1:7357680/93‑1 (MQ=255)
                                                                                       ttCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca       >  2:3569677/1‑93 (MQ=255)
                                                                                       ttCCGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGcaaca       <  2:4945522/93‑1 (MQ=255)
                                                                                         ccGGGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGAGTGGCGGCCACCATTCTGTTCCTGCaacaac     >  2:5791406/1‑93 (MQ=255)
                                                                                           gAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAAcgc   <  2:2538047/93‑1 (MQ=255)
                                                                                            |                                                                                            
CAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCG  >  NC_000913/2288972‑2289156

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