Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,717 +C coding (86/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,7171.C100.0% 1217.4 / NA 387T29R (ACG→AGG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (188/199);  total (188/199)

Only 100 of 446 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAAT  >  NC_000913/4001625‑4001809
                                                                                               |                                                                                            
tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                                               <  1:1561059/93‑1 (MQ=255)
 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgc                                                                                              <  2:9002031/93‑2 (MQ=255)
 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgc                                                                                              >  1:5323135/1‑92 (MQ=255)
 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgc                                                                                              <  2:6551169/93‑2 (MQ=255)
 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAATACgc                                                                                              >  1:8309601/1‑92 (MQ=255)
  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                             <  1:9730714/93‑3 (MQ=255)
   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgctg                                                                                            <  1:9509232/93‑4 (MQ=255)
      tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                                         >  2:5888111/1‑93 (MQ=255)
       caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcgc                                                                                        <  2:5295275/93‑1 (MQ=255)
        aGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt                                                                                       >  2:8722881/1‑93 (MQ=255)
         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCta                                                                                      <  1:8068102/93‑1 (MQ=255)
          caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtat                                                                                     >  1:9219809/1‑93 (MQ=255)
          caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACGCTGGCGCtat                                                                                     >  2:1890105/1‑93 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                                    <  2:5675382/93‑1 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                                    >  2:7922249/1‑93 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                                    >  2:4181691/1‑93 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                                    <  1:1669352/93‑1 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                                    >  1:6098288/1‑93 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAATAAATCAAATTGAAAGACGATGGCGCtata                                                                                    >  1:10004813/1‑93 (MQ=255)
            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                                   <  1:6021078/93‑1 (MQ=255)
            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGGACGCTGGCGCTATAc                                                                                   >  2:8198527/1‑93 (MQ=255)
            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTCGAAGTACGCTGGCGCTATAc                                                                                   <  2:7017374/93‑1 (MQ=255)
            ccaccaGAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                                   <  1:1215288/93‑1 (MQ=255)
             caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtataac                                                                                  >  2:3974453/1‑91 (MQ=255)
               ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcccc                                                                                <  2:981177/93‑1 (MQ=255)
                 aTAAAGAAAAACGACCAGATCACCCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTACGCTGGCGCTATACCCCaa                                                                              <  2:7220375/93‑1 (MQ=255)
                   aaaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaaa                                                                            >  1:424022/1‑92 (MQ=255)
                       aaaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCacataaaa                                                                        >  2:4464733/1‑93 (MQ=255)
                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaata                                                                      >  1:8775750/1‑93 (MQ=255)
                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAaataaaataa                                                                     <  1:4359015/93‑1 (MQ=255)
                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGCGCACGTAGTAAATCAACTTGAAGTACGCGGGCGCTATACCCCAAATAAAATAAGc                                                                   <  1:5737955/93‑1 (MQ=255)
                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAGGTAGGAAATCAACTTGAATTACACTGGCTCTATACACCAaataaaataaaag                                                                  >  1:2009258/1‑90 (MQ=255)
                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                                                  >  1:290855/1‑93 (MQ=255)
                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                                                 >  1:5894275/1‑93 (MQ=255)
                               ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTTGCGCTATACCCCAAATAAAATAAGCGGc                                                                >  2:4593638/1‑93 (MQ=255)
                               ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGc                                                                <  2:8289713/93‑1 (MQ=255)
                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCa                                                               >  2:4885251/1‑93 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGTGCATGTAGTAAATCAACTTGAAGTACGCTGGCGCTATCCCCCAAATAAAATAAGCGGCaa                                                              <  1:5238681/93‑1 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGGGCTATACCCCAAATAAAATAAGCGGCaa                                                              >  1:9977484/1‑93 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                              <  1:3206510/93‑1 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                              >  1:2151765/1‑93 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                              >  2:1092281/1‑93 (MQ=255)
                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaaa                                                             >  1:145224/1‑92 (MQ=255)
                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                             >  2:2548234/1‑93 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                            <  2:9385827/93‑1 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                            <  2:243603/93‑1 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                            >  2:8066183/1‑93 (MQ=255)
                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCg                                                         <  1:1932531/93‑1 (MQ=255)
                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                       <  2:3271008/93‑1 (MQ=255)
                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                       >  1:890952/1‑93 (MQ=255)
                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATAACTTTAAGGACCGTTGGGCTATAACCCAAATAAAAAAAACCGCGAGAACGAgg                                                      >  2:6131923/1‑92 (MQ=255)
                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                                      <  1:723307/93‑1 (MQ=255)
                                         cgATGCGTCAGGATTTAATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAga                                                      >  2:2184366/1‑92 (MQ=255)
                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                                    >  2:22660/1‑93 (MQ=255)
                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                                    >  1:5053143/1‑93 (MQ=255)
                                           aTGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                                    <  1:556028/93‑1 (MQ=255)
                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                                   <  1:7702115/93‑1 (MQ=255)
                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                                  <  2:1213056/93‑1 (MQ=255)
                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                                  <  2:6257125/93‑1 (MQ=255)
                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                                 <  1:6421224/93‑1 (MQ=255)
                                                tCAGGATTTCATCGGCGGGCTCGCAGTAAATCAACTTGAAGTACGCCTGCGCTATACCCCTAAAAAAATAAGCTGCTAGAGCGAGGAATACGc                                               >  2:9264854/1‑93 (MQ=255)
                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                               <  2:15796/93‑1 (MQ=255)
                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                               >  2:7079535/1‑93 (MQ=255)
                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                               <  2:3737298/93‑1 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAAGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                             >  1:4007498/1‑93 (MQ=255)
                                                   ggATTTCATCGGCGGGCCCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGGAAGAGCGAGGAATAGGccgt                                            >  2:1777783/1‑91 (MQ=255)
                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                            >  2:4349375/1‑93 (MQ=255)
                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                            >  1:4669314/1‑93 (MQ=255)
                                                    gATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                                           >  2:884113/1‑93 (MQ=255)
                                                     aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGc                                          <  2:9089001/93‑1 (MQ=255)
                                                     aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGc                                          >  2:291922/1‑93 (MQ=255)
                                                     aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGc                                          >  2:1373953/1‑93 (MQ=255)
                                                       ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                        >  2:7936720/1‑93 (MQ=255)
                                                        tCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc                                       <  2:96548/93‑1 (MQ=255)
                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                                    <  1:6854728/93‑1 (MQ=255)
                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                                    >  1:2602917/1‑93 (MQ=255)
                                                            cggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                                   <  1:9655462/93‑1 (MQ=255)
                                                             ggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                                  <  2:7623181/93‑1 (MQ=255)
                                                                   cACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAAGCCCTGCCGCGTTTGTTTTGc                            <  1:1849370/93‑1 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                           <  1:9152872/93‑1 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                           <  1:7132533/93‑1 (MQ=255)
                                                                     cGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                          >  2:6342110/1‑93 (MQ=255)
                                                                     cGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                          <  1:2364653/93‑1 (MQ=255)
                                                                      gtagtaAATCAACTTGAAGTACGCTGGCGCTATCCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                         <  2:6289524/93‑1 (MQ=255)
                                                                      gtagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                         <  1:6278341/93‑1 (MQ=255)
                                                                      gtagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                         <  1:54347/93‑1 (MQ=255)
                                                                          taAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc                     >  1:7349227/1‑93 (MQ=255)
                                                                              tCAAGTTGAAGTACTCTGGCGCTATACCGCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa                 <  1:5724745/93‑1 (MQ=255)
                                                                              tCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa                 <  1:6671170/93‑1 (MQ=255)
                                                                               cAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGaa                <  1:3522289/93‑1 (MQ=255)
                                                                               cAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGaa                <  2:7403956/93‑1 (MQ=255)
                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt               >  1:8518309/1‑93 (MQ=255)
                                                                                 aCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa              >  2:8089816/1‑93 (MQ=255)
                                                                                 aCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa              >  2:8383167/1‑93 (MQ=255)
                                                                                 aCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa              <  2:3480055/93‑1 (MQ=255)
                                                                                 aCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAACAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa              <  2:2895390/93‑1 (MQ=255)
                                                                                  cTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc             <  2:4064470/93‑1 (MQ=255)
                                                                                  cTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc             >  1:9424774/1‑93 (MQ=255)
                                                                                      aaGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAt         >  1:3650492/1‑93 (MQ=255)
                                                                                       aGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtt        <  2:5491839/93‑1 (MQ=255)
                                                                                               |                                                                                            
TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAAT  >  NC_000913/4001625‑4001809

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