Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,875,213 T→C intergenic (‑120/+120) b3698 ← / ← gyrB hypothetical protein/DNA gyrase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,875,2130TC98.9% 596.8 / ‑3.1 188intergenic (‑120/+120)b3698/gyrBhypothetical protein/DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (81/105);  total (82/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 200 total aligned reads displayed.

ACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGC  >  NC_000913/3875123‑3875302
                                                                                          |                                                                                           
aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGc                                                                                           <  1:3121424/93‑1 (MQ=255)
 cccAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGGTGATCTCCGCCCTTTTCTTATAGCTTCTCGcc                                                                                          >  1:533544/1‑93 (MQ=255)
   cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCgg                                                                                        >  2:953311/1‑93 (MQ=255)
   cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCgg                                                                                        >  2:1488618/1‑93 (MQ=255)
    aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGa                                                                                       >  2:515251/1‑93 (MQ=255)
    aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGa                                                                                       <  2:3295051/93‑1 (MQ=255)
     aaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGAt                                                                                      >  1:3309792/1‑93 (MQ=255)
       aaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgc                                                                                    <  1:3741996/93‑1 (MQ=255)
       aaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgc                                                                                    <  2:4170965/93‑1 (MQ=255)
       aaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgc                                                                                    <  1:1974615/93‑1 (MQ=255)
       aaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgc                                                                                    >  1:1209030/1‑93 (MQ=255)
         gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATCAAAAGGGCGGAGATGCACTCCGCCCTTTTCTTATAGCTTCTCGCCCGGTgcgg                                                                                  >  2:2523731/1‑93 (MQ=255)
         gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgcgg                                                                                  <  2:1661801/93‑1 (MQ=255)
         gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgcgg                                                                                  >  1:2815358/1‑93 (MQ=255)
          tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGGCGGCTgcggg                                                                                 >  2:1895515/1‑92 (MQ=39)
           ataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgcggcg                                                                                <  2:3842006/93‑1 (MQ=255)
            taGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGGTgcgggcg                                                                               >  1:984218/1‑90 (MQ=255)
            taGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGt                                                                               >  2:799892/1‑93 (MQ=255)
            taGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGt                                                                               <  2:4379420/93‑1 (MQ=255)
              gATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGa                                                                             >  1:4094100/1‑93 (MQ=255)
               aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGaa                                                                            <  2:115996/93‑1 (MQ=255)
               aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGaa                                                                            <  1:739174/93‑1 (MQ=255)
                tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAAc                                                                           <  1:2469912/93‑1 (MQ=255)
                tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAAc                                                                           >  2:1505818/1‑93 (MQ=255)
                   gCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGcc                                                                       >  1:1182001/1‑93 (MQ=255)
                     gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCtt                                                                    <  2:2205770/93‑1 (MQ=255)
                     gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCtt                                                                    <  1:2127402/93‑1 (MQ=255)
                     gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCtt                                                                    <  2:2347847/93‑1 (MQ=255)
                     gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCtt                                                                    >  1:286736/1‑93 (MQ=255)
                      ggATGATAATTGCGGATTGCATAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCCGATGCGGCGTGCACGCCttt                                                                   >  1:2616606/1‑92 (MQ=38)
                      ggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTa                                                                   <  2:1428441/93‑1 (MQ=255)
                       gatgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTAt                                                                  >  1:2046708/1‑93 (MQ=255)
                         tgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATcc                                                                >  1:363434/1‑93 (MQ=255)
                          gatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCg                                                               <  1:3802645/93‑1 (MQ=255)
                            tAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                             >  1:2444164/1‑93 (MQ=255)
                            tAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                             >  2:1525886/1‑93 (MQ=255)
                             aaTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                            <  1:1019627/93‑1 (MQ=255)
                               ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGAGGCGGCGTGAAACCCTTATCCCGCCtc                                                          >  1:1801611/1‑92 (MQ=38)
                                tGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                         <  1:2011776/93‑1 (MQ=255)
                                tGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCCGCCTaa                                                         >  1:611952/1‑92 (MQ=255)
                                  cGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                                                       <  2:2283571/93‑1 (MQ=255)
                                       tGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                  <  2:2072850/93‑1 (MQ=255)
                                         cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATg                                                <  2:1118208/93‑1 (MQ=255)
                                          ataaaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATTCGGGCCACCAAAACATGa                                               >  2:2088368/3‑93 (MQ=37)
                                          aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCCGATGCGGCGTGGACTCCCTTTTCCGGCCACAAAAACATGa                                               >  2:1953097/1‑93 (MQ=255)
                                           gAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGaa                                              >  1:834208/1‑93 (MQ=255)
                                           gAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACAAAATTATGaa                                              >  1:587635/1‑93 (MQ=255)
                                            aaaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCCGCGTGAACGCCTTTTCCCCGCTACAAAAACAttaca                                             >  2:3552259/1‑89 (MQ=37)
                                               tAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAAAGCCCTATTCGGGCTAAAAAAACATGaaaaat                                          >  2:3703721/1‑91 (MQ=25)
                                                aaaaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACAAAATCATGAAAAATc                                         >  1:4373329/1‑93 (MQ=255)
                                                   aaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAAt                                      >  1:3161711/1‑93 (MQ=255)
                                                    aGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTTTTCCGGCCTCAAAAAAAATAAAAAttaaca                                      >  2:3533623/1‑89 (MQ=25)
                                                    aGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATa                                     <  1:2721342/93‑1 (MQ=255)
                                                     ggtgggAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATAc                                    <  2:1460956/89‑1 (MQ=255)
                                                     gggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGGACCCCCTATTCCGGCTACCAAAACATTAAAAATCaaaaa                                    >  2:2636679/1‑90 (MQ=255)
                                                     gggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATAc                                    >  1:1275853/1‑93 (MQ=255)
                                                       gCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAt                                  <  1:420137/93‑1 (MQ=255)
                                                        cGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                 <  2:2691182/93‑1 (MQ=255)
                                                        cGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                 <  2:3999254/93‑1 (MQ=255)
                                                        cGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                 <  1:1957388/93‑1 (MQ=255)
                                                         ggAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                <  1:3409981/93‑1 (MQ=255)
                                                            gatgatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                             >  1:413950/1‑93 (MQ=255)
                                                            gatgatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                             <  2:2472910/93‑1 (MQ=255)
                                                               gatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAt                          <  2:2815230/93‑1 (MQ=255)
                                                                atCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtt                         >  2:734512/1‑93 (MQ=255)
                                                                 tctcCGCTCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAttt                        <  1:455691/93‑1 (MQ=255)
                                                                 tctcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAttt                        <  1:2241782/93‑1 (MQ=255)
                                                                  ctcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                       >  1:719920/1‑93 (MQ=255)
                                                                  ctcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                       >  1:4661363/1‑93 (MQ=255)
                                                                   tcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                      <  1:267342/93‑1 (MQ=255)
                                                                    ccGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCg                     >  1:4001214/1‑93 (MQ=255)
                                                                     cGCCCTTTTCTTATAGCTTCTCGCCGGCTGCGGCGTGAACCCCTTATCCCGCCTACAAAATCAGGAAAATTCCATACATTGCCAGATTTTCGt                    >  2:4628233/1‑93 (MQ=37)
                                                                     cGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGt                    <  2:2389271/93‑1 (MQ=255)
                                                                      gaccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:3423641/91‑1 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGGTGCGGCGTGGAAGCCCTATTCCGCCCACAAAAACATGAAAAATTCATAAAATTCAAGAATTTcttt                   >  2:89175/1‑90 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   >  1:2958987/1‑93 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  2:2314877/93‑1 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:3251060/93‑1 (MQ=255)
                                                                       cccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAg                  <  1:1593412/93‑1 (MQ=255)
                                                                        ccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAgg                 >  1:4034654/1‑93 (MQ=255)
                                                                         cTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTTGtttgc                >  2:3779664/1‑89 (MQ=21)
                                                                         cTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                <  2:2749878/93‑1 (MQ=255)
                                                                         cTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTACCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                <  1:3225357/93‑1 (MQ=255)
                                                                          ttttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGcc               <  1:1333641/93‑1 (MQ=255)
                                                                          ttttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGcc               <  2:4681196/93‑1 (MQ=255)
                                                                          ttttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGcc               <  1:778505/93‑1 (MQ=255)
                                                                          ttttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGcc               <  2:4423725/93‑1 (MQ=255)
                                                                           tttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCt              <  1:4066982/93‑1 (MQ=255)
                                                                           tttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCt              <  1:4362818/93‑1 (MQ=255)
                                                                           tttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCt              <  2:4248598/93‑1 (MQ=255)
                                                                           tttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGcca              <  2:2851423/93‑2 (MQ=255)
                                                                              cTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGAt           <  1:1131087/93‑1 (MQ=255)
                                                                               ttATAGCTTCTCGCCGGATGCGGCGTGGACGCCTTATCCGGGCTACAAAATCATGAAAAATCAATACATTGCAAGATTTTCTTAGGCCTGATa          >  1:2997452/1‑93 (MQ=18)
                                                                               ttATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATa          <  2:4313681/93‑1 (MQ=255)
                                                                                tataGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATaa         >  2:3225896/1‑93 (MQ=255)
                                                                                tataGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATaa         <  1:3944202/93‑1 (MQ=255)
                                                                                  taGCTTCTCGCCGGATGCGGCGTGAACGCCTTATTCGGGCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTTGTAGGCCTGATAAGc       >  2:3024594/1‑93 (MQ=18)
                                                                                   aGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACAAAATCATGAAAACTCAATACATTGCAAGATTTTCGTAGGCATGAGCAGCg      >  2:1560369/1‑93 (MQ=18)
                                                                                   aGCTTCTCGCCGGGTGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCg      >  2:2101472/1‑93 (MQ=17)
                                                                                   aGCTTCTCGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCg      >  2:2556066/1‑93 (MQ=255)
                                                                                          |                                                                                           
ACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGC  >  NC_000913/3875123‑3875302

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