Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 611,121 A→C S458R (AGT→AGG fepA ← ferrienterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913611,1210AC94.4% 392.5 / 9.2 125S458R (AGT→AGGfepAferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base C (63/55);  total (65/60)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 125 total aligned reads displayed.

GCGCTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTG  >  NC_000913/611024‑611215
                                                                                                 |                                                                                              
gcgcTGTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTACACCTTGCTAGATGTTCAGCGCCGGGCTCCAGTTCTTGCCGACAATCCt                                                                                              >  2:1270392/1‑100 (MQ=255)
 cgcTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCtg                                                                                             <  1:1152137/100‑1 (MQ=255)
 cgcTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCtg                                                                                             <  2:3135161/100‑1 (MQ=255)
    tttATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTCGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAAACCTGtga                                                                                          >  2:3505844/1‑100 (MQ=255)
     ttATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtgat                                                                                         <  1:1161303/100‑1 (MQ=255)
     ttATAAGCACGGGCGATGCCCATTTTCAGCGTCAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtgat                                                                                         <  1:3734124/100‑1 (MQ=255)
     ttATAAGCACCGACGATGCCCATTTTCAGCGTGAATTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTCtgat                                                                                         <  2:3395916/100‑1 (MQ=255)
      tataAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtgatg                                                                                        >  2:123462/1‑100 (MQ=255)
      tataAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATACTGtgatg                                                                                        >  2:1692861/1‑100 (MQ=255)
        taAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTgatgat                                                                                      >  1:2247316/1‑100 (MQ=255)
        taAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGACTCCAGTTATTGCCGCCAATCCTGtgatgaa                                                                                      >  2:637106/1‑99 (MQ=255)
         aaGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATc                                                                                     >  1:1366353/1‑100 (MQ=255)
         aaGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATc                                                                                     <  1:1967180/100‑1 (MQ=255)
         aaGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATc                                                                                     <  2:281918/100‑1 (MQ=255)
          aGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCg                                                                                    >  1:2481878/1‑100 (MQ=255)
          aGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCg                                                                                    >  1:1883966/1‑100 (MQ=255)
           gCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGa                                                                                   >  1:3467238/1‑100 (MQ=255)
              cGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAAc                                                                                >  2:3691705/1‑100 (MQ=255)
              cGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAAc                                                                                <  2:3500366/100‑1 (MQ=255)
                   gATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGc                                                                           <  2:140268/100‑1 (MQ=255)
                     tGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGccc                                                                         >  1:3227248/1‑100 (MQ=255)
                      gCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGcccc                                                                        >  2:1744652/1‑100 (MQ=255)
                      gCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGcccc                                                                        <  1:1621482/100‑1 (MQ=255)
                       cccATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCg                                                                       >  1:3391328/1‑100 (MQ=255)
                       cccATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCg                                                                       >  2:1004665/1‑100 (MQ=255)
                        ccATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCgg                                                                      >  2:3701962/1‑100 (MQ=255)
                        ccATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCgg                                                                      >  1:4017472/1‑100 (MQ=255)
                         cATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGc                                                                     <  1:129603/100‑1 (MQ=255)
                          aTTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCg                                                                    <  1:1252335/100‑1 (MQ=255)
                          aTTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCg                                                                    <  1:1247980/100‑1 (MQ=255)
                              tCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTAc                                                                <  1:1024459/100‑1 (MQ=255)
                                aGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGa                                                              >  1:1049578/1‑100 (MQ=255)
                                 gCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGAt                                                             <  2:2937361/100‑1 (MQ=255)
                                    tGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGt                                                          >  2:2392063/1‑100 (MQ=255)
                                     gAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTg                                                         <  2:3236158/100‑1 (MQ=255)
                                     gAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTg                                                         <  1:2063890/100‑1 (MQ=255)
                                      aaGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGc                                                        <  2:106596/100‑1 (MQ=255)
                                          cATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtc                                                    <  1:760274/100‑1 (MQ=255)
                                          cATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtc                                                    <  1:3832841/100‑1 (MQ=255)
                                           aTCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtca                                                   >  2:2230308/1‑100 (MQ=255)
                                           aTCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtca                                                   >  1:549990/1‑100 (MQ=255)
                                             cgcgTACACCTTGCGATATGTTCAGCGCTGGGCTATAGCTTTTGCCGACAATCGTGTGATGATCGAAACCCAGCCACGGCGTTACGATGGTGCTgtcagt                                                 <  2:538017/100‑1 (MQ=255)
                                              gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGTTgtcagtc                                                <  2:3586859/100‑1 (MQ=255)
                                              gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                >  1:1315113/1‑100 (MQ=255)
                                              gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                >  1:3743250/1‑100 (MQ=255)
                                               ccTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtca                                               >  2:651129/1‑100 (MQ=255)
                                               ccTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtca                                               <  1:1039566/100‑1 (MQ=255)
                                               ccTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtca                                               <  1:2264593/100‑1 (MQ=255)
                                                 tAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGc                                             >  1:2526736/1‑100 (MQ=255)
                                                 tAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGc                                             <  1:543501/100‑1 (MQ=255)
                                                  aaaCCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCt                                            <  1:1888791/100‑1 (MQ=255)
                                                    aCCTTGCGATATGTTCAGCGCCGGGCTCCAGGTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTcc                                          <  2:2142328/100‑1 (MQ=255)
                                                     ccTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCa                                         <  1:3377610/100‑1 (MQ=255)
                                                     ccTTGCGAGATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCCCCa                                         >  2:1369121/1‑100 (MQ=255)
                                                      cTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAt                                        >  2:1258386/1‑100 (MQ=255)
                                                      cTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAt                                        >  2:1258723/1‑100 (MQ=255)
                                                         gCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgtt                                     >  2:138286/1‑100 (MQ=255)
                                                         gCGATATGTTCAGCGCCGCGCTCCAGTTAATGCCGACAATCCTGTGATGATCAAAACGCAGCCCCGGCGTTACGATGATGCTGTCAGTCAGCTCCAtgtt                                     <  2:2328333/100‑1 (MQ=255)
                                                          cGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttg                                    >  1:262611/1‑100 (MQ=255)
                                                          cGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttg                                    <  2:2487890/100‑1 (MQ=255)
                                                               tGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTc                               >  1:1928005/1‑100 (MQ=255)
                                                                gTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCg                              <  2:780265/100‑1 (MQ=255)
                                                                gTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCg                              >  1:3290178/1‑100 (MQ=255)
                                                                 ttCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCgg                             <  2:3558504/100‑1 (MQ=255)
                                                                 ttCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCgg                             <  1:511240/100‑1 (MQ=255)
                                                                 ttCAGCGCCGGGCTCCAGTCACTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCgg                             >  1:3483234/1‑100 (MQ=255)
                                                                   cAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCa                           <  1:1853692/100‑1 (MQ=255)
                                                                   cAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCa                           >  2:3594936/1‑100 (MQ=255)
                                                                   cAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAAAGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGTCa                           >  2:3263115/1‑100 (MQ=255)
                                                                       gcCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACa                       <  1:286028/100‑1 (MQ=255)
                                                                       gcCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACa                       >  1:862721/1‑100 (MQ=255)
                                                                        ccGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAg                      <  2:1972251/100‑1 (MQ=255)
                                                                        ccGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAg                      <  2:417660/100‑1 (MQ=255)
                                                                        ccGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAg                      >  2:1380565/1‑100 (MQ=255)
                                                                        ccGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAg                      >  1:467400/1‑100 (MQ=255)
                                                                          gggCTCCATTTATTGCCGACAAGCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCg                    <  2:813754/100‑1 (MQ=255)
                                                                          gggCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCg                    <  1:2936374/100‑1 (MQ=255)
                                                                           ggCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCga                   >  2:2751100/1‑100 (MQ=255)
                                                                             cTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCgaga                 >  2:1521584/1‑100 (MQ=255)
                                                                              tctaGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGaa                >  2:734488/4‑100 (MQ=255)
                                                                              tCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAGCAGGGAGaa                >  2:847015/1‑100 (MQ=255)
                                                                              tCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGaa                <  1:3185401/100‑1 (MQ=255)
                                                                               ccAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGaaa               >  1:2376694/1‑100 (MQ=255)
                                                                                cAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGaaaa              >  1:1942721/1‑100 (MQ=255)
                                                                                 aGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAAt             >  2:1975718/1‑100 (MQ=255)
                                                                                  nttATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAAtt            >  1:3408171/2‑100 (MQ=255)
                                                                                  gTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAAtt            >  1:2892751/1‑100 (MQ=255)
                                                                                    tATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTc          <  2:2050492/100‑1 (MQ=255)
                                                                                     nttGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCt         >  1:3513144/2‑100 (MQ=255)
                                                                                     aTTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAAGTTGTTTTCAGCAAACAGAGAGAAAATTTCt         >  1:1266716/1‑100 (MQ=255)
                                                                                      ttGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTg        <  2:3561884/100‑1 (MQ=255)
                                                                                       tGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGc       >  2:2744765/1‑100 (MQ=255)
                                                                                        gCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCt      >  1:1818519/1‑100 (MQ=255)
                                                                                        gCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCt      >  2:161335/1‑100 (MQ=255)
                                                                                        gCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGCTTTCGGCAAACAGCGAGAAAATTTCTGCt      >  2:3516126/1‑100 (MQ=255)
                                                                                          cGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCttt    >  2:3633974/1‑100 (MQ=255)
                                                                                          cGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCttt    <  2:1897309/100‑1 (MQ=255)
                                                                                          cGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCttt    <  1:1135643/100‑1 (MQ=255)
                                                                                            aCAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTTTGCTTTTg  <  2:3832640/100‑1 (MQ=255)
                                                                                            aCAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTg  <  1:2793374/100‑1 (MQ=255)
                                                                                                 |                                                                                              
GCGCTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTG  >  NC_000913/611024‑611215

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