Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,488 G→A A471A (GCC→GCT yidE ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,4880GA100.0% 836.1 / NA 225A471A (GCC→GCTyidEputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (103/122);  total (103/122)

Only 100 of 228 total aligned reads displayed.

GCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGGGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGAA  >  NC_000913/3862399‑3862577
                                                                                         |                                                                                         
gCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCa                                                                                        <  2:2536368/93‑1 (MQ=255)
gCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCa                                                                                        <  2:5937318/93‑1 (MQ=255)
gCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCa                                                                                        <  2:4366392/93‑1 (MQ=255)
 ccAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCAc                                                                                       <  1:4170619/93‑1 (MQ=255)
  cAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCAcc                                                                                      >  2:3342739/1‑93 (MQ=255)
   aGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCa                                                                                     >  1:5259646/1‑93 (MQ=255)
   aGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCa                                                                                     >  2:5427522/1‑93 (MQ=255)
   aGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCa                                                                                     >  1:4834053/1‑93 (MQ=255)
   aGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGAACCa                                                                                     >  1:2884988/1‑93 (MQ=255)
    gCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCAt                                                                                    >  2:3926369/1‑93 (MQ=255)
     cATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATa                                                                                   >  2:4603592/1‑93 (MQ=255)
        cccGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAAcc                                                                                >  2:2262166/1‑93 (MQ=255)
          cGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCaa                                                                              <  1:37375/93‑1 (MQ=255)
          cGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCaa                                                                              >  1:1175097/1‑93 (MQ=255)
             acaTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATcc                                                                           >  2:3539173/1‑93 (MQ=255)
              caTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCa                                                                          <  2:4789713/93‑1 (MQ=255)
                tGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGc                                                                        >  1:1081150/1‑93 (MQ=255)
                tGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGc                                                                        >  2:2816924/1‑93 (MQ=255)
                 ggTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCt                                                                       <  1:4937964/93‑1 (MQ=255)
                  gTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCtt                                                                      >  1:877238/1‑93 (MQ=255)
                  gTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCtt                                                                      >  1:409750/1‑93 (MQ=255)
                   tCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTa                                                                     <  1:657802/93‑1 (MQ=255)
                   tCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCACCTTa                                                                     >  1:1863615/1‑93 (MQ=255)
                    cAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAg                                                                    <  1:5147165/93‑1 (MQ=255)
                     aGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGc                                                                   <  2:902743/93‑1 (MQ=255)
                      ggTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGcc                                                                  >  1:4356423/1‑93 (MQ=255)
                       gTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGccc                                                                 >  2:578705/1‑93 (MQ=255)
                        tAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCt                                                                >  1:2430854/1‑93 (MQ=255)
                        tAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCt                                                                <  1:1889951/93‑1 (MQ=255)
                          aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                              >  2:751561/1‑93 (MQ=255)
                          aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                              >  1:2180906/1‑93 (MQ=255)
                          aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGCAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                              >  2:2978573/1‑93 (MQ=255)
                             cATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGcc                                                           >  2:2663165/1‑93 (MQ=255)
                               ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                         <  1:2469738/93‑1 (MQ=255)
                               ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                         >  1:928889/1‑93 (MQ=255)
                                tttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAtt                                                        <  1:3720362/93‑1 (MQ=255)
                                 ttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTg                                                       <  2:1877962/93‑1 (MQ=255)
                                    gCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAcc                                                    <  2:3367925/93‑1 (MQ=255)
                                     cTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCa                                                   >  1:5756720/1‑93 (MQ=255)
                                     cTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCa                                                   <  1:1969551/93‑1 (MQ=255)
                                      tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                  <  1:6044966/93‑1 (MQ=255)
                                      tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                  <  2:4220332/93‑1 (MQ=255)
                                      tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                  <  2:1180601/93‑1 (MQ=255)
                                      tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                  <  2:2767510/93‑1 (MQ=255)
                                      tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                  <  1:2595028/93‑1 (MQ=255)
                                       aaCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAgg                                                 >  1:5498346/1‑93 (MQ=255)
                                       aaCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAgg                                                 <  1:3597909/93‑1 (MQ=255)
                                        aaaTCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAggg                                                >  2:4549949/3‑93 (MQ=255)
                                        aCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAggg                                                >  1:745375/1‑93 (MQ=255)
                                        aCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAggg                                                <  1:4561012/93‑1 (MQ=255)
                                           tCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTAt                                             >  2:179803/1‑93 (MQ=255)
                                             cgcgCCAGAATGCCAACAGTAATCAGCGGAACGGCGGCGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGCATTc                                           >  2:1704889/1‑93 (MQ=255)
                                             cgcgCCAGAATGCCAACAGTAATCAGCAGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTc                                           <  1:3991013/93‑1 (MQ=255)
                                              gcgcCAGAATGCCAACAGTAATCAGCGGATCGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTca                                          <  2:2138733/93‑1 (MQ=255)
                                                gcCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTcaca                                        >  2:3308902/1‑93 (MQ=255)
                                                 ccAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaa                                       <  1:1093570/93‑1 (MQ=255)
                                                 ccAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaa                                       >  1:4590982/1‑93 (MQ=255)
                                                  cAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaa                                      >  1:2148671/1‑93 (MQ=255)
                                                   aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                     >  1:1813718/1‑93 (MQ=255)
                                                   aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                     >  2:4181589/1‑93 (MQ=255)
                                                   aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                     >  1:5574995/1‑93 (MQ=255)
                                                    gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                    <  2:5931610/93‑1 (MQ=255)
                                                    gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                    >  2:4728704/1‑93 (MQ=255)
                                                    gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                    >  2:6066196/1‑93 (MQ=255)
                                                    gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                    >  2:2626729/1‑93 (MQ=255)
                                                     aaagCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATc                                   <  2:5518006/90‑1 (MQ=255)
                                                     aaTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATc                                   <  2:3599238/93‑1 (MQ=255)
                                                     aaTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATc                                   <  1:2296216/93‑1 (MQ=255)
                                                     aaTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATc                                   <  1:6278307/93‑1 (MQ=255)
                                                     aaTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATc                                   <  1:3777485/93‑1 (MQ=255)
                                                       tGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATccc                                 <  1:308670/93‑1 (MQ=255)
                                                        gCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATcccc                                <  2:1501129/93‑1 (MQ=255)
                                                        gCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATcccc                                >  1:2398161/1‑93 (MQ=255)
                                                         ccAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCcca                               <  1:2365111/93‑1 (MQ=255)
                                                         ccAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCcca                               <  1:1298403/93‑1 (MQ=255)
                                                          cAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccac                              <  2:1054990/93‑1 (MQ=255)
                                                          cAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccac                              <  1:5061668/93‑1 (MQ=255)
                                                          cAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccac                              <  1:5377543/93‑1 (MQ=255)
                                                            aCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacca                            >  2:3697279/1‑93 (MQ=255)
                                                            aCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacca                            >  2:5738177/1‑93 (MQ=255)
                                                            aCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacca                            <  1:4231012/93‑1 (MQ=255)
                                                             cAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAg                           >  1:3449311/1‑93 (MQ=255)
                                                             cAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAg                           >  2:3773566/1‑93 (MQ=255)
                                                              aGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGa                          >  1:3933961/1‑93 (MQ=255)
                                                               gTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAt                         >  1:5323221/1‑93 (MQ=255)
                                                                tAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAtt                        <  1:1648973/93‑1 (MQ=255)
                                                                  aTCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTc                      <  1:3029036/93‑1 (MQ=255)
                                                                   tCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCa                     <  2:1561813/93‑1 (MQ=255)
                                                                      gCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGAc                  >  1:4702359/1‑93 (MQ=255)
                                                                       cGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGAcc                 <  2:5283144/93‑1 (MQ=255)
                                                                        ggAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCa                <  1:1512898/93‑1 (MQ=255)
                                                                         gAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCaa               >  1:3373112/1‑93 (MQ=255)
                                                                            cggcggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGa            >  2:592316/1‑93 (MQ=255)
                                                                             ggcggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGAc           <  1:1388484/93‑1 (MQ=255)
                                                                              gcggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGAcc          >  2:2843323/1‑93 (MQ=255)
                                                                                ggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCga        >  1:5166237/1‑93 (MQ=255)
                                                                                ggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCga        <  2:4897663/93‑1 (MQ=255)
                                                                                ggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCga        >  2:3971399/1‑93 (MQ=255)
                                                                                ggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCga        <  1:3802984/93‑1 (MQ=255)
                                                                                   gATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGagag     <  1:4111916/93‑1 (MQ=255)
                                                                                         |                                                                                         
GCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGGGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGAA  >  NC_000913/3862399‑3862577

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