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% 701.2 / NA 208T29R (ACG→AGG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (89/119);  total (89/119)

Only 100 of 237 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803
                                                                                             |                                                                                      
tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                                         <  2:769867/93‑1 (MQ=255)
 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgc                                                                                        >  2:3719086/1‑92 (MQ=255)
  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                       <  1:2334993/93‑3 (MQ=255)
   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTACgctg                                                                                      >  1:5571483/1‑90 (MQ=255)
   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGCAGTACgctg                                                                                      >  1:2697908/1‑90 (MQ=255)
   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACgctg                                                                                      >  1:4502285/1‑90 (MQ=255)
    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTACgctgg                                                                                     >  1:2244818/1‑89 (MQ=255)
       caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcgc                                                                                  <  1:668225/93‑1 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                              >  2:1529801/1‑93 (MQ=255)
            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                             <  1:4074664/93‑1 (MQ=255)
             caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcc                                                                            >  1:2511515/1‑93 (MQ=255)
             caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAGTAAATCAACTTGAAGTACGCTGGCGCTATAcc                                                                            >  2:3280661/1‑93 (MQ=255)
              accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACGCTGGCGCTATAccc                                                                           >  1:5315777/1‑93 (MQ=255)
              accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATAAACTTGAAGTACGCTGGCGCTATAccc                                                                           >  2:1130891/1‑93 (MQ=255)
               ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcccc                                                                          <  2:4581840/93‑1 (MQ=255)
                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTACGCTTGCGCTATAAACCaa                                                                        >  1:3986185/1‑93 (MQ=255)
                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaa                                                                        <  2:1782113/93‑1 (MQ=255)
                   aaaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCATGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaat                                                                      <  2:4521770/93‑1 (MQ=255)
                   aaaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaat                                                                      >  2:4649552/1‑93 (MQ=255)
                    aaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATTAAATTTAAGGACGCTGGCGCTATACCCCcaata                                                                     >  1:1292443/1‑93 (MQ=255)
                    aaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaata                                                                     >  1:4227844/1‑93 (MQ=255)
                     aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataa                                                                    >  1:1852290/1‑93 (MQ=255)
                     aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataa                                                                    >  2:1238685/1‑93 (MQ=255)
                     aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataa                                                                    >  2:3637254/1‑93 (MQ=255)
                         aaaaGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaata                                                                >  1:4332770/1‑93 (MQ=255)
                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaata                                                                >  2:5313152/1‑93 (MQ=255)
                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCCaataaaataa                                                               >  1:1109869/1‑93 (MQ=255)
                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGc                                                             >  2:2760008/1‑93 (MQ=255)
                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTTGCGCTATACCCCAAATAAAATAAGCg                                                            >  2:235718/1‑93 (MQ=255)
                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                                            <  2:2821078/93‑1 (MQ=255)
                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                                            <  2:6202890/93‑1 (MQ=255)
                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGAAGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                                           >  1:2866687/1‑93 (MQ=255)
                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                                           >  1:2935924/1‑93 (MQ=255)
                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCggga                                                         >  1:2130998/1‑91 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        <  2:2493796/93‑1 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        <  1:344248/93‑1 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTAGGCTGGCGCTATACCCCCAATAAAATAAGAGACaa                                                        >  2:2895775/1‑93 (MQ=255)
                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                       >  2:2768634/1‑93 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTTGCGCTATACCCCAAATAAAATAAGCCGCAaga                                                      >  2:6423293/1‑93 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      >  2:5902706/1‑93 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  2:1318858/93‑1 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  1:3546099/93‑1 (MQ=255)
                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                     >  2:5647640/1‑93 (MQ=255)
                                     cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGGAagagg                                                    >  1:1664362/1‑92 (MQ=255)
                                     cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                    >  2:1085492/1‑93 (MQ=255)
                                     cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                    <  1:2801830/93‑1 (MQ=255)
                                       cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                  <  2:5695640/93‑1 (MQ=255)
                                       cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                  <  1:6086783/93‑1 (MQ=255)
                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                 <  1:791580/93‑1 (MQ=255)
                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                                >  2:1758716/1‑93 (MQ=255)
                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAgga                                               >  1:5621412/1‑91 (MQ=255)
                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             >  2:297208/1‑93 (MQ=255)
                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:6113743/93‑1 (MQ=255)
                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:889215/93‑1 (MQ=255)
                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATaag                                          >  1:3811734/1‑91 (MQ=255)
                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                          >  2:3273268/1‑93 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  2:2063268/93‑1 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  1:4687865/93‑1 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  2:1404864/93‑1 (MQ=255)
                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                      <  2:3391906/93‑1 (MQ=255)
                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                      <  2:3924470/93‑1 (MQ=255)
                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAATAATACGCCCt                                      >  1:573964/1‑93 (MQ=255)
                                                     aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGc                                    <  1:107507/93‑1 (MQ=255)
                                                         cATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgcg                                >  1:841393/1‑93 (MQ=255)
                                                          aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATCAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                               <  2:6305359/93‑1 (MQ=255)
                                                          aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                               <  1:5813957/93‑1 (MQ=255)
                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              <  1:1582714/93‑1 (MQ=255)
                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              >  2:4784692/1‑93 (MQ=255)
                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              <  1:5469392/93‑1 (MQ=255)
                                                             ggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                            <  1:5136956/93‑1 (MQ=255)
                                                               cggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                          <  1:706429/93‑1 (MQ=255)
                                                               cggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                          <  1:1006902/93‑1 (MQ=255)
                                                                gggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                         <  2:1294392/93‑1 (MQ=255)
                                                                 ggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATCCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                        <  2:1868041/93‑1 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     <  1:4826529/93‑1 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  1:2099285/1‑93 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  2:5126704/1‑93 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  1:4908919/1‑93 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  1:5259007/1‑93 (MQ=255)
                                                                     cGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                    >  1:3239725/1‑93 (MQ=255)
                                                                       tagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  <  1:241343/93‑1 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                <  1:1327457/93‑1 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                >  1:4565536/1‑93 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                <  2:4194247/93‑1 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                >  2:1967580/1‑93 (MQ=255)
                                                                          taAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc               >  1:1248442/1‑93 (MQ=255)
                                                                          taAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc               <  1:62217/93‑1 (MQ=255)
                                                                           aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg              <  1:2540025/93‑1 (MQ=255)
                                                                           aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg              >  1:1977997/1‑93 (MQ=255)
                                                                             aTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCggg            <  2:4981296/93‑1 (MQ=255)
                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         <  2:1435859/93‑1 (MQ=255)
                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         <  1:239913/93‑1 (MQ=255)
                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         <  2:6298092/93‑1 (MQ=255)
                                                                                  cTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       >  2:5396568/1‑93 (MQ=255)
                                                                                   ttGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACt      <  1:5303266/93‑1 (MQ=255)
                                                                                    tGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTc     >  1:4645502/1‑93 (MQ=255)
                                                                                     gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa    <  2:2583004/93‑1 (MQ=255)
                                                                                     gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa    <  2:3128251/93‑1 (MQ=255)
                                                                                     gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa    >  1:1909459/1‑93 (MQ=255)
                                                                                       aGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtt  >  1:6426501/1‑93 (MQ=255)
                                                                                             |                                                                                      
TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

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