Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,993,671 C→A T840K (ACG→AAG)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,993,6710CA99.2% 414.0 / NA 118T840K (ACG→AAG) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (69/47);  minor base G (0/1);  total (70/48)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

TACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGG  >  NC_000913/3993534‑3993803
                                                                                                                                         |                                                                                                                                    
tACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAagc                                                                                                                                    <  1:794120/140‑1 (MQ=255)
tACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAagc                                                                                                                                    >  2:410891/1‑140 (MQ=255)
 aCCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAagcc                                                                                                                                   <  1:2651756/140‑1 (MQ=255)
     aCGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTa                                                                                                                               >  1:1772404/1‑140 (MQ=255)
       ggCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAagcc                                                                                                                                   <  1:524709/134‑1 (MQ=255)
       ggCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTAtt                                                                                                                             >  1:410879/1‑140 (MQ=255)
        gCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTa                                                                                                                            <  2:2647736/140‑1 (MQ=255)
         cTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTAc                                                                                                                           >  2:2639570/1‑140 (MQ=255)
             aGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCa                                                                                                                       >  2:1232743/1‑140 (MQ=255)
                 tcatcaACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAatat                                                                                                                   <  1:410891/140‑1 (MQ=255)
                  catcaACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAagc                                                                                                                                    >  2:788904/1‑122 (MQ=255)
                  catcaACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATAtt                                                                                                                  >  2:1905814/1‑140 (MQ=255)
                   atcaACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATAttt                                                                                                                 <  1:2639570/140‑1 (MQ=255)
                         ttCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtg                                                                                                           >  2:121421/1‑140 (MQ=255)
                          tCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGtga                                                                                                          >  1:1675702/1‑140 (MQ=255)
                            aaCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTCCAGCAATATTTTTCGtgatg                                                                                                        >  2:1230773/1‑140 (MQ=255)
                            aaCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGtgatg                                                                                                        >  1:2119185/1‑140 (MQ=255)
                             aCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGtgatga                                                                                                       <  1:459878/140‑1 (MQ=255)
                                  ccGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTg                                                                                                  >  1:1720707/1‑140 (MQ=255)
                                  ccGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTg                                                                                                  >  2:1323546/1‑140 (MQ=255)
                                   cGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGc                                                                                                 >  1:1550231/1‑140 (MQ=255)
                                     cAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCg                                                                                               <  2:2444474/140‑1 (MQ=255)
                                     cAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCg                                                                                               <  2:1586577/140‑1 (MQ=255)
                                     cAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCg                                                                                               <  2:1369310/140‑1 (MQ=255)
                                      aGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAACCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCgg                                                                                              >  1:1002950/1‑140 (MQ=255)
                                        ttCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGaa                                                                                            >  2:1858976/1‑140 (MQ=255)
                                           tATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGc                                                                                         <  1:1193593/140‑1 (MQ=255)
                                           tATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGc                                                                                         <  1:1389272/140‑1 (MQ=255)
                                                  ttGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACGTGTCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCtt                                                                                  >  1:2821807/1‑140 (MQ=255)
                                                  ttGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAg                                                                                          <  2:2040102/132‑1 (MQ=255)
                                                  ttGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAg                                                                                          >  2:1844682/1‑132 (MQ=255)
                                                  ttGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAg                                                                                          >  1:2040102/1‑132 (MQ=255)
                                                  ttGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAg                                                                                          <  1:1844682/132‑1 (MQ=255)
                                                  ttGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCtt                                                                                  >  2:1224462/1‑140 (MQ=255)
                                                  ttGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCtt                                                                                  >  2:1457381/1‑140 (MQ=255)
                                                   tgtgAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTa                                                                                 >  1:1485720/1‑140 (MQ=255)
                                                    gtgAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTAt                                                                                >  2:622903/1‑140 (MQ=255)
                                                     tgAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATc                                                                               >  1:2393370/1‑140 (MQ=255)
                                                     tgAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATc                                                                               >  2:2275286/1‑140 (MQ=255)
                                                       aaGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCg                                                                             <  1:1230773/140‑1 (MQ=255)
                                                       aaGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCg                                                                             <  2:410879/140‑1 (MQ=255)
                                                            tGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATg                                                                        >  2:1509976/1‑140 (MQ=255)
                                                             gATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGc                                                                       >  1:1645164/1‑140 (MQ=255)
                                                             gATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGc                                                                       >  1:2364536/1‑140 (MQ=255)
                                                               tGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACATGCCGGAAAGCGAGGCTTATCCGGCATGCaa                                                                     >  1:1421038/1‑140 (MQ=255)
                                                                   cGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCtt                                                                 <  2:445324/140‑1 (MQ=255)
                                                                    gTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTa                                                                >  2:1173644/1‑140 (MQ=255)
                                                                     tGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAg                                                               <  2:1720707/140‑1 (MQ=255)
                                                                       aaCAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCg                                                             >  2:2843583/1‑140 (MQ=255)
                                                                       aaCAGGCGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCg                                                             >  2:2361856/1‑140 (MQ=255)
                                                                        aCAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCgg                                                            <  2:2814050/140‑1 (MQ=255)
                                                                        aCAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCgg                                                            <  1:622903/140‑1 (MQ=255)
                                                                           ggTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGaaa                                                         <  1:2729694/140‑1 (MQ=255)
                                                                            gTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGaa                                                          <  1:659556/138‑1 (MQ=255)
                                                                            gTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGaa                                                          >  2:659556/1‑138 (MQ=255)
                                                                            gTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAac                                                        >  2:1091393/1‑140 (MQ=255)
                                                                            gTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAac                                                        >  1:1645330/1‑140 (MQ=255)
                                                                            gTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAac                                                        >  1:470009/1‑140 (MQ=255)
                                                                              gATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAg                                                               >  2:793350/1‑131 (MQ=255)
                                                                              gATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAg                                                               <  1:793350/131‑1 (MQ=255)
                                                                              gATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAcctg                                                      <  2:1333181/140‑1 (MQ=255)
                                                                              gATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAACCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactg                                                      <  2:1002950/140‑1 (MQ=255)
                                                                                  ccgttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgactg                                                  <  1:1905814/140‑1 (MQ=255)
                                                                                    gttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAg                                                               >  2:925517/1‑125 (MQ=255)
                                                                                    gttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAg                                                               <  1:925517/125‑1 (MQ=255)
                                                                                     ttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGttt                                               <  2:1645164/140‑1 (MQ=255)
                                                                                     ttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGttt                                               >  1:41989/1‑140 (MQ=255)
                                                                                      tccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTc                                              >  2:16419/1‑140 (MQ=255)
                                                                                      tccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTc                                              >  2:2318179/1‑140 (MQ=255)
                                                                                      tccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTc                                              >  1:2146533/1‑140 (MQ=255)
                                                                                       ccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCa                                             <  1:1317680/140‑1 (MQ=255)
                                                                                          cacaAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAccc                                          >  1:185882/1‑140 (MQ=255)
                                                                                            caAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGc                                        <  1:824683/140‑1 (MQ=255)
                                                                                                 tCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgct                                   <  1:852703/140‑1 (MQ=255)
                                                                                                  cTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgctg                                  <  2:2504659/140‑1 (MQ=255)
                                                                                                     tCGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGt                               <  1:874953/140‑1 (MQ=255)
                                                                                                      cGGTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTc                              <  1:2843583/140‑1 (MQ=255)
                                                                                                       ggTAACATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTcg                             <  1:2318179/140‑1 (MQ=255)
                                                                                                            cATGCCGCCTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCgccgcc                        <  2:432627/140‑1 (MQ=255)
                                                                                                               gccgccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGt                     >  1:2145567/1‑140 (MQ=255)
                                                                                                               gccgccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGt                     >  2:1540671/1‑140 (MQ=255)
                                                                                                               gccgccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGt                     >  1:1724243/1‑140 (MQ=255)
                                                                                                               gccgccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGt                     >  2:2035037/1‑140 (MQ=255)
                                                                                                                 cgccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTc                   <  2:2746733/140‑1 (MQ=255)
                                                                                                                   ccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGTGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCa                 >  1:966123/1‑140 (MQ=255)
                                                                                                                   ccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCa                 >  2:2291970/1‑140 (MQ=255)
                                                                                                                   ccTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCa                 >  1:1122679/1‑140 (MQ=255)
                                                                                                                    cTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAg                <  2:470009/140‑1 (MQ=255)
                                                                                                                    cTGCCAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAg                <  1:1759282/140‑1 (MQ=255)
                                                                                                                      gCCAATCAGTATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCa              <  2:547644/140‑1 (MQ=255)
                                                                                                                      gCCAATCAGGATCACGATAGGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCa              <  2:406746/140‑1 (MQ=255)
                                                                                                                       ccAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTcgcc                           >  2:2301275/1‑126 (MQ=255)
                                                                                                                       ccAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTcgcc                           <  1:2301275/126‑1 (MQ=255)
                                                                                                                       ccAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCaa             >  1:2069547/1‑140 (MQ=255)
                                                                                                                       ccAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCaa             <  1:1870844/140‑1 (MQ=255)
                                                                                                                       ccAATCAGGATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCaa             <  1:1457381/140‑1 (MQ=255)
                                                                                                                               gATCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGGCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGa     <  2:326090/140‑1 (MQ=255)
                                                                                                                                aTCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGaa    >  1:972927/1‑140 (MQ=255)
                                                                                                                                 tCACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAg      <  2:2115669/137‑1 (MQ=255)
                                                                                                                                  cACGATAAGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAgg  >  1:299138/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                    
TACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGG  >  NC_000913/3993534‑3993803

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