Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 621,281 A→C I303M (ATT→ATG)  fepD ← ferrienterobactin ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,2810AC94.6% 465.4 / 10.1 149I303M (ATT→ATG) fepDferrienterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base C (65/76);  total (69/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 151 total aligned reads displayed.

AATCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACG  >  NC_000913/621183‑621379
                                                                                                  |                                                                                                  
aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACaa                                                                                                   >  2:2065120/1‑100 (MQ=255)
 ntCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTCACCACAGAAACGCGCAGTCCGCCGGGCACCAt                                                                                                  >  1:3246282/2‑100 (MQ=255)
    aTGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCAc                                                                                               >  2:852065/1‑100 (MQ=255)
    aTGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCAc                                                                                               >  1:1100589/1‑100 (MQ=255)
     ttcaccacCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACTACAGAAACGCGCAGTTCGCCGGGCACCAGCAcc                                                                                              >  2:182003/3‑100 (MQ=255)
      gCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCGccc                                                                                             <  1:1654929/100‑1 (MQ=255)
      gCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCAccc                                                                                             >  2:2879665/1‑100 (MQ=255)
        accaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCAcccgc                                                                                           >  1:230610/1‑100 (MQ=255)
        accaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTACACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCAcccgc                                                                                           <  1:1857943/100‑1 (MQ=255)
         ccaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCAcccgcc                                                                                          <  1:3455353/100‑1 (MQ=255)
          caccTCGCGTATTACGTCGCACGAGGAAGATCAGCACCGGGGCGCCAATAAACGCACTGACCACAGAAACGCGCAGTACGCCGGGCACCATCAcctgttt                                                                                         <  2:3598410/100‑6 (MQ=255)
           accTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCCCTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCAcccgcccg                                                                                        <  2:513207/100‑1 (MQ=255)
             cTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCgat                                                                                      >  2:2742392/1‑100 (MQ=255)
             cTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCgat                                                                                      >  2:2270527/1‑100 (MQ=255)
                  gTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGatat                                                                                 >  1:955180/1‑100 (MQ=255)
                  gTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGatat                                                                                 <  2:1281650/100‑1 (MQ=255)
                    tttACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCg                                                                               <  1:1076135/100‑1 (MQ=255)
                      tACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGc                                                                             >  2:1547818/1‑100 (MQ=255)
                      tACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGc                                                                             >  1:641889/1‑100 (MQ=255)
                        cGTCGCACGAGGGAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAATTCGCCGGGCACCATCACCCGCCCGATGATATCGGCaa                                                                           >  2:3850456/1‑100 (MQ=255)
                        cGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCaa                                                                           <  1:602782/100‑1 (MQ=255)
                        cGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCaa                                                                           >  1:2715488/1‑100 (MQ=255)
                        cGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCaa                                                                           <  1:2612294/100‑1 (MQ=255)
                          tCGCACGAGGAAGATCTGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCAGGCACCATCTCCCGCCCGATGATATCGGCAAAc                                                                         <  2:3346026/100‑1 (MQ=255)
                           cGCACGCGGAAGATCAGCACCCGTGCACCAATAAACGCACTGACCACAGAAACGGGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAACca                                                                        >  2:3810589/1‑100 (MQ=255)
                            gCACGAGGAAGATCAGCACAGGTGCACAAATAAACGCCCTGACCACAGAAACCCACAGGTAGCGTGTCACAATCCCCGGCGACATGCTATCGACAATcag                                                                       >  2:3305880/1‑100 (MQ=255)
                             cACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAAcagc                                                                      >  2:2524258/1‑100 (MQ=255)
                               cGAGGAAGATCAGCACCGGTGCGCCAATAAACGCACTGACCACCGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAAcagcag                                                                    <  2:1482354/100‑1 (MQ=255)
                                 aGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAAcagcagca                                                                  >  2:2278848/1‑100 (MQ=255)
                                  ggAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAAcagcagcag                                                                 <  1:2923493/100‑1 (MQ=255)
                                        tCAGCACCGGTGCACCAATAAAGGCACTGACCACAGAAACGCGCAGGGCGCCGGGCACCAGCACCCGCCCGATGATATCGGCAAACAGCAgcagggcagg                                                           <  2:785864/100‑1 (MQ=255)
                                        tCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAgcagggcagg                                                           >  1:3216921/1‑100 (MQ=255)
                                         cAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGcagggcaggg                                                          <  1:2563536/100‑1 (MQ=255)
                                          aGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAgggg                                                         >  1:825535/1‑100 (MQ=255)
                                           gCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGt                                                        >  1:1118994/1‑100 (MQ=255)
                                           gCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGt                                                        <  1:497295/100‑1 (MQ=255)
                                           gCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGAAATCGGCAAACAGCAGCAGGGCAGGGGt                                                        <  1:3444076/100‑1 (MQ=255)
                                           gCACCGGTGCACCAAAAAACGCACTGAGCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGt                                                        <  2:3373134/100‑1 (MQ=255)
                                             aCCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATCTCGGCAAACAGCGGCAGGGCGGGGGGAg                                                      >  2:749470/1‑100 (MQ=255)
                                              ccGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGc                                                     <  1:1637489/100‑1 (MQ=255)
                                               cgttccACCAATAAACTCACCGACCACATAAACTCGCAAGTCGCCAGGCACCATGACTCGCCCGATGAGAGCGCCAAACAGCAGCAGGGCAGGGGTAGCa                                                    >  2:3993867/6‑100 (MQ=255)
                                               cGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCa                                                    <  2:2166820/100‑1 (MQ=255)
                                                ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCaa                                                   >  2:3556786/1‑100 (MQ=255)
                                                 gggCACCAATAAACGGACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAg                                                  <  2:2574435/98‑1 (MQ=255)
                                                    cACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAg                                               <  2:5372/100‑1 (MQ=255)
                                                    cACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAg                                               >  2:424886/1‑100 (MQ=255)
                                                    cACCAATAAACGCACTGACCACAGAAACGCGCAGTTCACCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAg                                               >  2:3912964/1‑100 (MQ=255)
                                                     aCCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGc                                              <  1:1377559/100‑1 (MQ=255)
                                                     aCCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGc                                              >  2:3701949/1‑100 (MQ=255)
                                                     aCCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGc                                              >  1:2255530/1‑100 (MQ=255)
                                                      ccAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCg                                             >  1:3149288/1‑100 (MQ=255)
                                                      ccAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCg                                             <  1:2924047/100‑1 (MQ=255)
                                                      ccAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCg                                             <  1:953008/100‑1 (MQ=255)
                                                       cAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGt                                            <  1:1084135/100‑1 (MQ=255)
                                                          tagaCGCACAGACCACAGAAACGCGCAGTGCGCCGGGCACCATCACCCGGCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGAc                                         <  2:787267/97‑1 (MQ=255)
                                                          tAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGAc                                         <  2:2940226/100‑1 (MQ=255)
                                                          tAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGAc                                         <  1:3520602/100‑1 (MQ=255)
                                                           aaaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACg                                        >  1:813819/1‑100 (MQ=255)
                                                            caCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACgg                                       <  2:750544/99‑1 (MQ=255)
                                                             aCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACggg                                      <  2:608483/100‑1 (MQ=255)
                                                              cGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGc                                     <  1:3441416/100‑1 (MQ=255)
                                                               gCACTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCa                                    >  2:3014688/1‑100 (MQ=255)
                                                                 aCTGACCACAGAAACGCGCATTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGc                                  >  1:3183339/1‑100 (MQ=255)
                                                                 aCTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGc                                  <  1:1216395/100‑1 (MQ=255)
                                                                  cTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCg                                 >  1:2528557/1‑100 (MQ=255)
                                                                  cTGACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCg                                 <  2:1954167/100‑1 (MQ=255)
                                                                    gACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGAc                               >  1:1551410/1‑100 (MQ=255)
                                                                    gACCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGAc                               >  1:837281/1‑100 (MQ=255)
                                                                     aCCACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGAcc                              >  1:2610247/1‑100 (MQ=255)
                                                                      ccACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCa                             >  2:393833/1‑100 (MQ=255)
                                                                      ccACAGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCa                             >  1:2083212/1‑100 (MQ=255)
                                                                       cacaGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAg                            >  1:1300918/1‑100 (MQ=255)
                                                                       cacaGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAg                            >  1:917102/1‑100 (MQ=255)
                                                                       cacaGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAg                            >  2:3859946/1‑100 (MQ=255)
                                                                        acaGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGc                           <  1:3405053/100‑1 (MQ=255)
                                                                        acaGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGc                           >  2:665089/1‑100 (MQ=255)
                                                                          aGAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCga                         >  1:3649902/1‑100 (MQ=255)
                                                                           gAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCgat                        >  2:664058/1‑100 (MQ=255)
                                                                           gAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCgat                        >  2:1006922/1‑100 (MQ=255)
                                                                           gAAACGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCgat                        >  1:2284637/1‑100 (MQ=255)
                                                                              aCGCGCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCgatgat                     <  1:3611280/100‑1 (MQ=255)
                                                                               cgcgCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATc                    >  1:769752/1‑100 (MQ=255)
                                                                               cgcgCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATc                    >  1:3321437/1‑100 (MQ=255)
                                                                               cgcgCAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAACAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATc                    <  1:2336891/100‑1 (MQ=255)
                                                                                 cgcAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCgg                  <  1:3901565/100‑1 (MQ=255)
                                                                                   cAGTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGca                <  1:4081881/100‑1 (MQ=255)
                                                                                     gTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGcacc              <  1:1743387/100‑1 (MQ=255)
                                                                                     gTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGcacc              >  1:77972/1‑100 (MQ=255)
                                                                                     gTTCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGcacc              >  1:1996296/1‑100 (MQ=255)
                                                                                       tCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGcaccca            <  1:2691342/100‑1 (MQ=255)
                                                                                       tCGCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGcaccca            <  2:61467/100‑1 (MQ=255)
                                                                                         gCCGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGcacccacc          <  1:399684/100‑1 (MQ=255)
                                                                                          ccGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACccacca         <  1:433019/100‑1 (MQ=255)
                                                                                           cGGGCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAg        >  1:1043961/1‑100 (MQ=255)
                                                                                            gggCACCCTCGCCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGc       <  1:1128938/100‑1 (MQ=255)
                                                                                            gggCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGc       <  2:1495631/100‑1 (MQ=255)
                                                                                             ggCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGcc      >  2:327569/1‑100 (MQ=255)
                                                                                             ggCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGcc      <  1:3855330/100‑1 (MQ=255)
                                                                                              gCACCATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCa     <  2:1993191/100‑1 (MQ=255)
                                                                                                 caATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAAcg  <  2:1354447/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
AATCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACG  >  NC_000913/621183‑621379

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