Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 169,953 G→A R59H (CGT→CAT)  fhuC → iron‑hydroxamate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_000091169,9530GA100.0% 500.7 / NA 127R59H (CGT→CAT) fhuCiron‑hydroxamate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (60/67);  total (60/67)

Only 100 of 128 total aligned reads displayed.

TGTCGTTAACCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCGTCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGG  >  AC_000091/169866‑170039
                                                                                       |                                                                                      
tGTCGTTAACCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCAtc                                                                                      <  2:796850/90‑1 (MQ=255)
 gTCGTTAACCTTTCCTGCCTGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCAtca                                                                                     <  1:2181961/90‑1 (MQ=255)
 gTCGTTAACCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCAtca                                                                                     >  1:784187/1‑90 (MQ=255)
   cGTTAACCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCatcatc                                                                                   <  1:152587/90‑1 (MQ=255)
    gTTAACCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCAtcatca                                                                                  <  1:86988/90‑1 (MQ=255)
    gTTAACCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCAtcatca                                                                                  <  1:577561/90‑1 (MQ=255)
        aCCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAgccg                                                                              >  1:1682396/1‑90 (MQ=255)
        aCCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAgccg                                                                              <  1:1927501/90‑1 (MQ=255)
        aCCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAgccg                                                                              <  2:132046/90‑1 (MQ=255)
          cTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAgccgcc                                                                            <  1:3329752/90‑1 (MQ=255)
          cTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAgccgcc                                                                            <  1:3414148/90‑1 (MQ=255)
          cTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAgccgcc                                                                            <  1:100191/90‑1 (MQ=255)
           tttCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGccgccg                                                                           <  1:2102558/90‑1 (MQ=255)
              ccTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCg                                                                        <  1:1036521/90‑1 (MQ=255)
              ccTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCg                                                                        <  2:2804530/90‑1 (MQ=255)
              ccTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGATCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCg                                                                        >  2:1802862/1‑90 (MQ=255)
               cTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCgg                                                                       >  2:3297877/1‑90 (MQ=255)
               cTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCgg                                                                       >  2:159920/1‑90 (MQ=255)
               cTGCCGGGAAAGGGACCGGTCTGATTGGGCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCgg                                                                       >  2:1960736/1‑90 (MQ=255)
                tGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGa                                                                      <  2:1580736/90‑1 (MQ=255)
                tGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGa                                                                      >  1:1095097/1‑90 (MQ=255)
                tGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGa                                                                      <  1:2978818/90‑1 (MQ=255)
                tGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGa                                                                      >  2:340025/1‑90 (MQ=255)
                 gCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGaa                                                                     <  2:1615092/90‑1 (MQ=255)
                 gCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGaa                                                                     <  2:2622916/90‑1 (MQ=255)
                    gggAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAggg                                                                  >  2:660007/1‑90 (MQ=255)
                    gggAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAggg                                                                  >  1:3175157/1‑90 (MQ=255)
                     ggAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAgggg                                                                 >  2:2200056/1‑90 (MQ=255)
                     ggAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAgggg                                                                 <  1:1797963/90‑1 (MQ=255)
                     ggAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAgggg                                                                 <  1:49480/90‑1 (MQ=255)
                      gAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGga                                                                >  2:2118451/1‑90 (MQ=255)
                      gAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGga                                                                >  2:1218149/1‑90 (MQ=255)
                          gTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGAtt                                                            <  1:988811/90‑1 (MQ=255)
                          gTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGAtt                                                            <  1:2904998/90‑1 (MQ=255)
                            gACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGAttct                                                          >  1:3293228/1‑90 (MQ=255)
                            gACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGAttct                                                          <  1:2871517/90‑1 (MQ=255)
                             aCCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGAttctt                                                         <  1:1938054/90‑1 (MQ=255)
                             aCCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGAttctt                                                         >  2:2813566/1‑90 (MQ=255)
                              ccGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGAttcttc                                                        <  2:2818535/90‑1 (MQ=255)
                               cGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATtcttct                                                       <  1:2520095/90‑1 (MQ=255)
                                 gTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCTCCGTCGGAAGGGGAGATTCTTCTTg                                                     <  1:233505/90‑1 (MQ=255)
                                      aTTGGTCACAACGGTTCTGGTAAATCCACTCCGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGcc                                                <  1:1927354/90‑1 (MQ=255)
                                       ttGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGccc                                               <  2:466808/90‑1 (MQ=255)
                                         ggTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCaa                                             >  1:1695286/1‑90 (MQ=255)
                                          gTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAAc                                            >  1:2199496/1‑90 (MQ=255)
                                          gTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAAc                                            >  1:261622/1‑90 (MQ=255)
                                           tCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAAcc                                           <  1:588868/90‑1 (MQ=255)
                                           tCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAAcc                                           <  2:1951081/90‑1 (MQ=255)
                                            cacaACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCg                                          >  2:1521955/1‑90 (MQ=255)
                                            cacaACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCg                                          >  2:2306585/1‑90 (MQ=255)
                                               aaCGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTg                                       >  2:1810081/1‑90 (MQ=255)
                                                 cGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGa                                     >  2:2221128/1‑90 (MQ=255)
                                                   gTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGaaa                                   >  2:1215802/1‑90 (MQ=255)
                                                     tCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGc                                 >  1:1450863/1‑90 (MQ=255)
                                                      cTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCt                                <  2:3064515/90‑1 (MQ=255)
                                                      cTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCt                                >  2:1296782/1‑90 (MQ=255)
                                                      cTGGTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCt                                >  1:3083784/1‑90 (MQ=255)
                                                         gTAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGa                             <  1:2200056/90‑1 (MQ=255)
                                                          tAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGag                            <  1:2611290/90‑1 (MQ=255)
                                                          tAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGag                            <  2:1497068/90‑1 (MQ=255)
                                                          tAAATCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGag                            <  2:412605/90‑1 (MQ=255)
                                                            aaTCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGggca                          <  2:107624/90‑1 (MQ=255)
                                                            aaTCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGagca                          >  2:2685956/1‑90 (MQ=255)
                                                             aTCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGagcag                         >  2:2020409/1‑90 (MQ=255)
                                                             aTCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGagcag                         <  1:3037983/90‑1 (MQ=255)
                                                             aTCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGagcag                         <  1:3421202/90‑1 (MQ=255)
                                                             aTCCACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGagcag                         <  2:1313564/90‑1 (MQ=255)
                                                               ccACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAgcagca                       <  1:3368860/90‑1 (MQ=255)
                                                                cACTCTGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGGAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCaa                      >  2:2641909/1‑90 (MQ=255)
                                                                   tctGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGc                   <  1:412244/90‑1 (MQ=255)
                                                                   tctGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGc                   <  1:711645/90‑1 (MQ=255)
                                                                   tctGCTCAAAATCCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGc                   >  2:2523311/1‑90 (MQ=255)
                                                                    ctGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCg                  >  1:2019469/1‑90 (MQ=255)
                                                                    ctGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCg                  >  1:3254985/1‑90 (MQ=255)
                                                                    ctGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCg                  <  1:972235/90‑1 (MQ=255)
                                                                     tGCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGt                 >  2:3275484/1‑90 (MQ=255)
                                                                      gCTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGtt                <  1:216034/90‑1 (MQ=255)
                                                                       gtCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGttt               >  2:1810785/2‑90 (MQ=255)
                                                                       cTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGttt               <  1:566839/90‑1 (MQ=255)
                                                                       cTCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGttt               <  1:354481/90‑1 (MQ=255)
                                                                        tCAAAATGCTTGGCCATCATCAGGCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTg              >  2:2554917/1‑90 (MQ=255)
                                                                        tCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTg              <  2:680627/90‑1 (MQ=255)
                                                                        tCAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTg              <  1:1711354/90‑1 (MQ=255)
                                                                         cAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGc             <  2:1419920/90‑1 (MQ=255)
                                                                         cAAAATGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGc             <  2:694095/90‑1 (MQ=255)
                                                                            aaTGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCg          <  2:2347486/90‑1 (MQ=255)
                                                                             aTGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGc         >  1:926610/1‑90 (MQ=255)
                                                                             aTGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGc         <  1:592162/90‑1 (MQ=255)
                                                                              tGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGGCCGCa        >  1:1497669/1‑90 (MQ=255)
                                                                              tGCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCa        >  1:770227/1‑90 (MQ=255)
                                                                               gCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCaa       >  2:1568267/1‑90 (MQ=255)
                                                                               gCTTGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCaa       <  1:2038063/90‑1 (MQ=255)
                                                                                 ttGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAg     >  2:610039/1‑90 (MQ=255)
                                                                                  tGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGt    <  1:3232755/90‑1 (MQ=255)
                                                                                  tGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGt    >  1:1329606/1‑90 (MQ=255)
                                                                                  tGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGt    <  2:1544701/90‑1 (MQ=255)
                                                                                  tGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGt    >  1:2494588/1‑90 (MQ=255)
                                                                                  tGGCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGt    <  2:222458/90‑1 (MQ=255)
                                                                                   ggCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTg   >  1:1299841/1‑90 (MQ=255)
                                                                                    gCCATCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTgg  >  2:1058119/1‑90 (MQ=255)
                                                                                       |                                                                                      
TGTCGTTAACCTTTCCTGCCGGGAAAGTGACCGGTCTGATTGGTCACAACGGTTCTGGTAAATCCACTCTGCTCAAAATGCTTGGCCGTCATCAGCCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGG  >  AC_000091/169866‑170039

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