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

Only 100 of 252 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAT  >  NC_000913/4001625‑4001802
                                                                                             |                                                                                     
tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                                        >  1:2175927/1‑93 (MQ=255)
 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgc                                                                                       >  2:3326783/1‑92 (MQ=255)
  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                      >  2:2824202/1‑91 (MQ=255)
  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                      <  1:1341771/93‑3 (MQ=255)
  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACgct                                                                                      >  1:5465104/1‑91 (MQ=255)
   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgctg                                                                                     >  1:3945805/1‑90 (MQ=255)
   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGAACgctg                                                                                     >  1:4882795/1‑90 (MQ=255)
    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgctgg                                                                                    >  1:4076158/1‑89 (MQ=255)
        aGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt                                                                                <  1:6244454/93‑1 (MQ=255)
          caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtat                                                                              <  1:600101/93‑1 (MQ=255)
           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAATTTGAAGTACGCTGGCGCtata                                                                             <  2:4008949/93‑1 (MQ=255)
            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                            >  2:4559731/1‑93 (MQ=255)
            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                            <  2:1601262/93‑1 (MQ=255)
              accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAccc                                                                          <  2:1589871/93‑1 (MQ=255)
              accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAccc                                                                          >  1:5884042/1‑93 (MQ=255)
                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTCGTCAATCAACTTGAAGCAGGCTGGCGCTATACCCCa                                                                        >  1:5261649/1‑93 (MQ=255)
                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACGCTGGCGCTATACCCCa                                                                        >  2:2071528/1‑93 (MQ=255)
                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaa                                                                       >  2:2239886/1‑93 (MQ=255)
                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTACGCTGGCGCTATACCCCaaa                                                                      >  2:883195/1‑93 (MQ=255)
                    aaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaata                                                                    >  1:2930120/1‑93 (MQ=255)
                       aaaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaa                                                                 <  2:4349163/93‑1 (MQ=255)
                           aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                                             <  1:3747335/93‑1 (MQ=255)
                           aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                                             <  1:596674/93‑1 (MQ=255)
                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGc                                                            <  1:2623905/93‑1 (MQ=255)
                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGc                                                            <  1:3182348/93‑1 (MQ=255)
                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGc                                                            <  1:4175820/93‑1 (MQ=255)
                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                                           <  1:4769510/93‑1 (MQ=255)
                             caCCAGCTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                                           <  2:2875549/92‑1 (MQ=255)
                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGTAGTAAATCAACTTGAAGGACGCTTGCGCTATACCCCAAATAAAATAAgatg                                                          >  2:1143691/1‑90 (MQ=255)
                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAAGCTGGCGCTATACCCCAAATAAAAAAAGAgg                                                          >  1:880211/1‑93 (MQ=255)
                               ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGc                                                         <  2:1263059/93‑1 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                       <  1:4574604/93‑1 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                       >  1:6131571/1‑93 (MQ=255)
                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                       <  1:296887/93‑1 (MQ=255)
                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                      <  2:745848/93‑1 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                     <  1:4385216/93‑1 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                     >  1:1677336/1‑93 (MQ=255)
                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                     <  2:5156459/93‑1 (MQ=255)
                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                    <  2:6028138/93‑1 (MQ=255)
                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                    <  1:3852176/93‑1 (MQ=255)
                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                    <  2:1517323/93‑1 (MQ=255)
                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                    <  2:2321566/93‑1 (MQ=255)
                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                    <  1:5229805/93‑1 (MQ=255)
                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGAACGCTGGCGCTATACCCCAAATAAAATAAGCGGCCagag                                                    >  1:1111880/1‑93 (MQ=255)
                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCg                                                  >  2:4790543/1‑93 (MQ=255)
                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCg                                                  <  2:3405366/93‑1 (MQ=255)
                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCg                                                  >  2:2713790/1‑93 (MQ=255)
                                       cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                 <  1:1185854/93‑1 (MQ=255)
                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                <  2:4760080/93‑1 (MQ=255)
                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATTAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAACGGCAAGAGCGAGt                                               >  2:5367409/1‑93 (MQ=255)
                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                               <  2:5516121/93‑1 (MQ=255)
                                         cgATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAga                                               >  1:3316467/1‑92 (MQ=255)
                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                             <  2:5956425/93‑1 (MQ=255)
                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                             >  1:4990858/1‑93 (MQ=255)
                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAAGCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                           <  1:3985786/93‑1 (MQ=255)
                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                          <  2:5423424/93‑1 (MQ=255)
                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                         >  1:5622932/1‑93 (MQ=255)
                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                        <  2:4602204/93‑1 (MQ=255)
                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                        <  2:5894619/93‑1 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                      >  2:4456612/1‑93 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                      >  2:239511/1‑93 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                      <  1:5703704/93‑1 (MQ=255)
                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                      <  2:3112108/93‑1 (MQ=255)
                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                     <  2:3541685/93‑1 (MQ=255)
                                                    gATTTCATCGGCGGGGACGTAGTAAATCAACCTGGAGTACGCTGGGGGTATACCCCAAATAAAATAAACGGCAAGAGCGAGCAATACGCCCTg                                    >  2:3668330/1‑93 (MQ=255)
                                                     aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGc                                   <  1:509286/93‑1 (MQ=255)
                                                     aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGc                                   <  2:4458950/93‑1 (MQ=255)
                                                      tttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcc                                  >  2:3459155/1‑93 (MQ=255)
                                                          aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                              <  2:5471307/93‑1 (MQ=255)
                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                             <  2:4366726/93‑1 (MQ=255)
                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                             <  2:2550856/93‑1 (MQ=255)
                                                             ggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                           >  1:2742820/1‑93 (MQ=255)
                                                             ggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                           >  1:1000585/1‑93 (MQ=255)
                                                              gcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgt                          <  2:6221888/93‑1 (MQ=255)
                                                                gggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                        <  2:2424507/93‑1 (MQ=255)
                                                                gggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                        <  1:1566877/93‑1 (MQ=255)
                                                                 ggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                       <  2:137959/93‑1 (MQ=255)
                                                                 ggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                       >  2:3968062/1‑93 (MQ=255)
                                                                   cACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc                     <  2:951136/93‑1 (MQ=255)
                                                                   cACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc                     >  1:5910080/1‑93 (MQ=255)
                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                    <  1:3619060/93‑1 (MQ=255)
                                                                       tagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                 >  2:2955843/1‑93 (MQ=255)
                                                                        agtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                >  2:5704627/1‑93 (MQ=255)
                                                                        agtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                <  1:5495114/93‑1 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat               >  2:4175933/1‑93 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat               <  2:6255785/93‑1 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat               <  1:2314409/93‑1 (MQ=255)
                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat               >  1:4633673/1‑93 (MQ=255)
                                                                          taAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc              >  2:2473377/1‑93 (MQ=255)
                                                                           aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg             >  1:761367/1‑93 (MQ=255)
                                                                           aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg             >  1:3005048/1‑93 (MQ=255)
                                                                           aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg             >  2:2019096/1‑93 (MQ=255)
                                                                           aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg             >  1:1969386/1‑93 (MQ=255)
                                                                            aaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATACATCgg            <  1:1776864/93‑1 (MQ=255)
                                                                              tCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa          <  1:3423242/93‑1 (MQ=255)
                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt        <  2:5137965/93‑1 (MQ=255)
                                                                                 aCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa       <  2:1310191/93‑1 (MQ=255)
                                                                                    tGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTc    <  1:1859331/93‑1 (MQ=255)
                                                                                     gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa   <  1:1051849/93‑1 (MQ=255)
                                                                                      aaGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAt  <  1:4399208/93‑1 (MQ=255)
                                                                                             |                                                                                     
TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAT  >  NC_000913/4001625‑4001802

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