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,1210AC99.4% 621.3 / ‑4.3 172S458R (AGT→AGGfepAferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (80/91);  total (81/91)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 172 total aligned reads displayed.

CGCTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGA  >  NC_000913/611025‑611216
                                                                                                |                                                                                               
cgcTTTATAAGCACGGGGGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCtg                                                                                              <  2:4658478/100‑1 (MQ=255)
 gcTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCtgt                                                                                             >  2:5321030/1‑100 (MQ=255)
  cTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACGATCCtgtg                                                                                            <  2:5272898/100‑1 (MQ=255)
  cTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCtgtg                                                                                            >  1:94348/1‑100 (MQ=255)
   tttATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtga                                                                                           <  1:4284613/100‑1 (MQ=255)
   tttATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCGTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtga                                                                                           <  2:1134446/100‑1 (MQ=255)
    ttATAAGCACGGGCGATGCCCATTTTCCGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtgat                                                                                          <  1:5455634/100‑1 (MQ=255)
     tataAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGGTATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtgatg                                                                                         >  1:501819/1‑100 (MQ=255)
     tataAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtgatg                                                                                         <  1:3462678/100‑1 (MQ=255)
     tataAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAAACTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGtgatg                                                                                         >  1:5433932/1‑100 (MQ=255)
        aaGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACGATCCTGTGATGATc                                                                                      <  2:2672914/100‑1 (MQ=255)
        aaGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATc                                                                                      >  1:2492587/1‑100 (MQ=255)
         aGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCg                                                                                     >  1:2711236/1‑100 (MQ=255)
          gCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCgt                                                                                    <  2:2896962/100‑2 (MQ=255)
          gCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGa                                                                                    <  1:1320840/100‑1 (MQ=255)
               agCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGc                                                                               >  1:4240682/2‑100 (MQ=255)
                gCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCTATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCa                                                                              <  1:2146664/100‑1 (MQ=255)
                  gATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGc                                                                            >  1:26915/1‑100 (MQ=255)
                     gCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGcccc                                                                         >  1:3417980/1‑100 (MQ=255)
                      cccATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCg                                                                        >  2:5381439/1‑100 (MQ=255)
                       ccATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCgg                                                                       >  2:3703728/1‑100 (MQ=255)
                        cATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGc                                                                      >  2:3548532/1‑100 (MQ=255)
                        cATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGc                                                                      >  2:178844/1‑100 (MQ=255)
                         aTTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCg                                                                     >  1:2857985/1‑100 (MQ=255)
                         aTTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCg                                                                     >  1:1956541/1‑100 (MQ=255)
                             tCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCATCCCCGGCGTTAc                                                                 <  1:4785253/100‑1 (MQ=255)
                             tCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTAc                                                                 <  1:3570650/100‑1 (MQ=255)
                             tCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGGTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTAc                                                                 <  1:4598176/100‑1 (MQ=255)
                                gCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGAt                                                              >  2:3261990/1‑100 (MQ=255)
                                gCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGAt                                                              <  2:2251597/100‑1 (MQ=255)
                                  gTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATgg                                                            <  1:2375149/100‑1 (MQ=255)
                                    gAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTg                                                          >  1:699833/1‑100 (MQ=255)
                                     aaGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAAGCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGc                                                         <  2:689718/100‑1 (MQ=255)
                                      aGTCATCGCCTACACCGTGCGATAGTTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATTATCTAAACGCAGCCCCGGCGTTACGATGTTGcc                                                        >  2:3326311/1‑99 (MQ=255)
                                       gTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTg                                                       <  1:3407857/100‑1 (MQ=255)
                                       gTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTg                                                       <  1:3966223/100‑1 (MQ=255)
                                         cATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtc                                                     <  2:517580/100‑1 (MQ=255)
                                         cATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtc                                                     >  1:4729642/1‑100 (MQ=255)
                                          aTCGCCTAAACCTTGCGACATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtca                                                    >  2:4970571/1‑100 (MQ=255)
                                           tCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcag                                                   >  1:1686347/1‑100 (MQ=255)
                                            cGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagt                                                  <  2:2165517/100‑1 (MQ=255)
                                            cGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagt                                                  >  2:2392510/1‑100 (MQ=255)
                                            cGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagt                                                  <  1:1829608/100‑1 (MQ=255)
                                             gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                 <  1:5081174/100‑1 (MQ=255)
                                             gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                 <  2:3046297/100‑1 (MQ=255)
                                             gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                 >  2:1283186/1‑100 (MQ=255)
                                             gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                 >  1:2175589/1‑100 (MQ=255)
                                             gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                 >  1:2423717/1‑100 (MQ=255)
                                             gCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtc                                                 >  2:3817031/1‑100 (MQ=255)
                                              ccTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTgtcagtca                                                <  1:1171217/100‑1 (MQ=255)
                                                 aaaCCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCt                                             <  1:1424318/100‑1 (MQ=255)
                                                    ccTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCa                                          >  1:2089915/1‑100 (MQ=255)
                                                     cTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAt                                         >  2:4374083/1‑100 (MQ=255)
                                                     cTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAt                                         <  2:3181572/100‑1 (MQ=255)
                                                     cTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAt                                         <  1:3838961/100‑1 (MQ=255)
                                                      ttGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtg                                        <  2:4535167/100‑1 (MQ=255)
                                                       tGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgt                                       <  2:2479810/100‑1 (MQ=255)
                                                         cGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttg                                     >  1:149691/1‑100 (MQ=255)
                                                         cGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttg                                     >  1:782814/1‑100 (MQ=255)
                                                          gATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGCTCGCAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttgt                                    <  2:1214305/100‑1 (MQ=255)
                                                          gATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttgt                                    >  2:1614642/1‑100 (MQ=255)
                                                          gATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttgt                                    <  1:4258509/100‑1 (MQ=255)
                                                          gATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCAtgttgt                                    >  2:2012145/1‑100 (MQ=255)
                                                               gTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCg                               >  2:1642432/1‑100 (MQ=255)
                                                               gTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCg                               >  1:800263/1‑100 (MQ=255)
                                                               gTTCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCg                               <  2:2118939/100‑1 (MQ=255)
                                                               gTTCAGCGCCGGACTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCg                               <  1:572691/100‑1 (MQ=255)
                                                                ttCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTACATGTTGTTTTCgg                              <  1:3581630/100‑1 (MQ=255)
                                                                 tCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGc                             <  1:1201741/100‑1 (MQ=255)
                                                                 tCAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGc                             <  1:1439831/100‑1 (MQ=255)
                                                                  cAGCGCCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCa                            >  2:4092410/1‑100 (MQ=255)
                                                                   aGCGCCGGGCTCCAGTTATTGCCGTCAATCCTGTGATGAGCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCaa                           >  2:5223185/1‑100 (MQ=255)
                                                                    gcgcCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCaaa                          >  2:1675428/1‑100 (MQ=255)
                                                                    gcgcCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCaaa                          >  2:4869105/1‑100 (MQ=255)
                                                                     cgcCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAAc                         <  1:3892479/100‑1 (MQ=255)
                                                                      gcCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACa                        >  1:2808262/1‑100 (MQ=255)
                                                                      gcCGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACa                        >  1:2808757/1‑100 (MQ=255)
                                                                       ccGGGCTCCAGTTATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAg                       >  1:2195129/1‑100 (MQ=255)
                                                                        cGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGc                      <  2:2378702/100‑1 (MQ=255)
                                                                        cGGGCTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGAAGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGc                      <  2:1998336/100‑1 (MQ=255)
                                                                            cTCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGATCCATGTTGTTTTCGGCAAAACGCgaga                  >  2:4667870/1‑100 (MQ=255)
                                                                             tCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGaa                 <  2:5393089/100‑1 (MQ=255)
                                                                             tCCAGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGaa                 <  2:345988/100‑1 (MQ=255)
                                                                                aGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAAt              >  2:3992173/1‑100 (MQ=255)
                                                                                aGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAAt              >  1:3293488/1‑100 (MQ=255)
                                                                                aGTTATTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAAt              <  1:2005666/100‑1 (MQ=255)
                                                                                 cgatctGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATCGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAAtt             <  1:4885249/95‑1 (MQ=255)
                                                                                    aTTGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCt          <  1:3738072/100‑1 (MQ=255)
                                                                                     ttGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTg         <  1:5393721/100‑1 (MQ=255)
                                                                                     ttGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTg         <  1:1892830/100‑1 (MQ=255)
                                                                                      tGCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGc        <  1:4069863/100‑1 (MQ=255)
                                                                                       gCCGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCt       >  1:5107806/1‑100 (MQ=255)
                                                                                        ccGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAACCAGCGAGAAAATTTCTGCtt      >  2:3909572/1‑100 (MQ=255)
                                                                                        ccGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCtt      >  1:2111299/1‑100 (MQ=255)
                                                                                         cGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCttt     <  2:4645960/100‑1 (MQ=255)
                                                                                         cGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCttt     >  2:359885/1‑100 (MQ=255)
                                                                                         cGACAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCttt     >  2:2831979/1‑100 (MQ=255)
                                                                                           aCAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTg   <  2:867815/100‑1 (MQ=255)
                                                                                           aCAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTg   <  2:1895191/100‑1 (MQ=255)
                                                                                           aCAATCCTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTg   <  2:2155882/100‑1 (MQ=255)
                                                                                                |                                                                                               
CGCTTTATAAGCACGGGCGATGCCCATTTTCAGCGTGAAGTCATCGCCTAAACCTTGCGATATGTTCAGCGCCGGGCTCCAGTTATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGA  >  NC_000913/611025‑611216

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: