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,2130TC99.5% 1233.8 / ‑4.4 395intergenic (‑120/+120)b3698/gyrBhypothetical protein/DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (189/203);  total (191/204)
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.93e-01

Only 100 of 429 total aligned reads displayed.

ACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGC  >  NC_000913/3875123‑3875302
                                                                                          |                                                                                         
   cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCgg                                                                                      >  1:685845/1‑93 (MQ=255)
    aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGa                                                                                     <  2:2807856/93‑1 (MQ=255)
     aaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGGGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGAt                                                                                    <  2:4637205/93‑1 (MQ=255)
     aaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCggtg                                                                                    >  2:4862157/1‑91 (MQ=255)
     aaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCCGAt                                                                                    >  2:2934663/1‑93 (MQ=255)
         gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgcgg                                                                                >  1:6075366/1‑93 (MQ=255)
         gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgcgg                                                                                <  2:2665623/93‑1 (MQ=255)
           ataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATgcggcg                                                                              >  1:7518068/1‑93 (MQ=255)
            taGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGGCGGGTgcggggg                                                                             >  2:511595/1‑90 (MQ=255)
               aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGaa                                                                          >  1:3311179/1‑93 (MQ=255)
               aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGaa                                                                          <  1:2989787/93‑1 (MQ=255)
               aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGaa                                                                          <  1:7773021/93‑1 (MQ=255)
               aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGaa                                                                          <  2:3129122/93‑1 (MQ=255)
                tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAAc                                                                         >  1:5199845/1‑93 (MQ=255)
                     gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGGACGccct                                                                    >  1:7241332/1‑91 (MQ=255)
                      ggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTa                                                                   >  1:1366547/1‑93 (MQ=255)
                       gatgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGGTGCGGCGTGGACGGCCTAt                                                                  >  2:4461085/1‑93 (MQ=255)
                       gatgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTAt                                                                  <  2:6080439/93‑1 (MQ=255)
                       gatgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTAt                                                                  >  2:3596401/1‑93 (MQ=255)
                       gatgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTAt                                                                  >  1:2612219/1‑93 (MQ=255)
                       gatgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAAAGcccttt                                                                  >  2:6032424/1‑89 (MQ=255)
                        atgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATc                                                                 <  1:7707238/93‑1 (MQ=255)
                        atgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATc                                                                 <  1:5758991/93‑1 (MQ=255)
                          gatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCg                                                               >  2:5307889/1‑93 (MQ=255)
                          gatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCg                                                               >  2:1188113/1‑93 (MQ=255)
                             aaTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                            <  2:376085/93‑1 (MQ=255)
                             aaTTGCGGATTGAAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                            <  2:1622692/93‑1 (MQ=255)
                              aTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                           >  1:4030350/1‑93 (MQ=255)
                                tGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                         >  1:6733216/1‑93 (MQ=255)
                                tGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                         <  1:646598/93‑1 (MQ=255)
                                 gCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                        <  1:2533316/93‑1 (MQ=255)
                                    gATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                     <  1:284217/93‑1 (MQ=255)
                                    gATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                     >  1:3860158/1‑93 (MQ=255)
                                    gATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGACCGATGCTGCGTGAACCACTTATCCCGCCTAccaga                                                     >  1:3521315/1‑89 (MQ=255)
                                     aTTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                    <  2:3487883/93‑1 (MQ=255)
                                     aTTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCCTATCCCGGCTACaaaaa                                                    >  2:522846/1‑92 (MQ=38)
                                      ttGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCCTATCCGGCCTACAAAATc                                                   >  1:5460566/1‑93 (MQ=255)
                                       tGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                  >  1:2748749/1‑93 (MQ=255)
                                       tGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                  <  1:5018332/93‑1 (MQ=255)
                                       tGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                  <  1:2192403/93‑1 (MQ=255)
                                         cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATg                                                >  2:2290720/1‑93 (MQ=255)
                                          aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGa                                               <  2:4725475/93‑1 (MQ=255)
                                          aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACAAAATCATGa                                               >  1:3804492/1‑93 (MQ=255)
                                           gAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGaa                                              >  2:1230141/1‑93 (MQ=255)
                                           gAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGaa                                              >  2:4200892/1‑93 (MQ=255)
                                             aaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGaaaa                                            >  1:1169842/1‑93 (MQ=255)
                                             aaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACAAAATAATGaaaa                                            >  1:4861086/1‑93 (MQ=255)
                                               tAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGACTTATCCGGCCTACAAAATCATGAAAAtt                                          >  2:6898088/1‑93 (MQ=255)
                                                aaaaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTc                                         <  2:3305743/93‑1 (MQ=255)
                                                 aaaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGGTGCGGCGTGGACCCCTTTTTCCGCCTTCAAAAACACGAAAATTCa                                        >  2:6182594/1‑93 (MQ=25)
                                                  aaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCaa                                       <  2:7566035/93‑1 (MQ=255)
                                                  aaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCaa                                       <  2:3538394/93‑1 (MQ=255)
                                                  aaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGAAAATTCaa                                       >  2:6290936/1‑93 (MQ=255)
                                                    aGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGGAAATTCAATa                                     >  2:6587699/1‑93 (MQ=255)
                                                      ggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAAACACCTATCCAGCCTCCAAAATCATGAAAATTCCATACa                                   >  2:853569/1‑93 (MQ=25)
                                                       gCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAt                                  <  1:5891657/93‑1 (MQ=255)
                                                        cGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                 <  2:3374658/93‑1 (MQ=255)
                                                         ggAGATGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAATGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                <  2:5508893/93‑1 (MQ=255)
                                                          gagaTGATCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGc                               <  2:6814994/93‑1 (MQ=255)
                                                            gatgatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                             <  2:1778002/93‑1 (MQ=255)
                                                             atgatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAg                            >  1:3454084/1‑93 (MQ=255)
                                                             atgatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAg                            <  1:5358703/93‑1 (MQ=255)
                                                               gatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACAAAATCATGAAAAATTAATACATTGCAAGAt                          >  1:4803348/1‑93 (MQ=25)
                                                               gatCTCCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAt                          >  1:1007350/1‑93 (MQ=255)
                                                                 tctcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAttt                        >  1:6197589/1‑93 (MQ=255)
                                                                  ctcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                       <  2:7086289/93‑1 (MQ=255)
                                                                   tcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                      <  2:6847464/93‑1 (MQ=255)
                                                                   tcCGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                      <  1:3738853/93‑1 (MQ=255)
                                                                     cGCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGt                    >  1:1287835/1‑93 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  2:2153306/93‑1 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:7171531/93‑1 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:4248051/93‑1 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:1706240/93‑1 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:2360776/93‑1 (MQ=255)
                                                                      gCCCTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:2210028/93‑1 (MQ=255)
                                                                      ccccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                   <  1:4375206/92‑1 (MQ=255)
                                                                       cccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAg                  <  2:7797240/93‑1 (MQ=255)
                                                                       cccTTTTCTTATAGCTTCTCGCCCGATGCGGCGTGAACCCCTTATTCCGGCTACAAAATCATTAAAATTCCATACAATGCAAGATTTTTGttg                  >  2:248596/1‑91 (MQ=21)
                                                                        ccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAgg                 >  2:1044402/1‑93 (MQ=255)
                                                                        ccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAgg                 <  2:7234169/93‑1 (MQ=255)
                                                                        ccTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAgg                 <  2:3822715/93‑1 (MQ=255)
                                                                         cTTTTCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                <  2:2538388/93‑1 (MQ=255)
                                                                           tttCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCt              >  1:2477756/1‑93 (MQ=255)
                                                                            ttCTTATAGCTTCTCGCCGGATGCGGCGTGAACACCTTATCCCGCCTACAAAATCATGAAAAATCAATACATTTCAAAATTTTCGTAGGcccg             >  2:5773912/1‑91 (MQ=21)
                                                                            ttCTTATACCTTCTCTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTg             <  2:5910447/93‑1 (MQ=255)
                                                                             tCTTATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGa            <  1:3627886/93‑1 (MQ=255)
                                                                               ttATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATa          >  1:2287586/1‑93 (MQ=255)
                                                                               ttATAGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATAAATTGCAAGATTTTCGTAGGCCTGATa          >  2:4648252/1‑93 (MQ=255)
                                                                                tataGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATaa         <  1:5554620/93‑1 (MQ=255)
                                                                                tataGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATaa         <  1:7896920/93‑1 (MQ=255)
                                                                                tataGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATaa         <  1:4512818/93‑1 (MQ=255)
                                                                                tataGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATaa         <  1:4724130/93‑1 (MQ=255)
                                                                                 ataGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAg        >  1:1989244/1‑93 (MQ=255)
                                                                                  taGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGc       <  1:3204496/93‑1 (MQ=255)
                                                                                   aGCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCg      <  2:7383335/93‑1 (MQ=255)
                                                                                    gCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt     >  2:3971159/1‑93 (MQ=255)
                                                                                    gCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt     >  1:4908660/1‑93 (MQ=255)
                                                                                    gCTTCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt     <  1:7464683/93‑1 (MQ=255)
                                                                                    gCTTCTCGCCGGATGCGGCGTGAACGCCCTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCgg     >  1:1881467/1‑92 (MQ=255)
                                                                                      ttctcGCCGGATGCGGCGTGAACCCCCTATCCCGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAg   >  2:1906245/1‑93 (MQ=17)
                                                                                          |                                                                                         
ACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGC  >  NC_000913/3875123‑3875302

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