Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,496,885 G→A intergenic (+396/‑1387) leuX → / → insC6 tRNA‑Leu/IS2 insertion element repressor InsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,496,8850GA97.6% 125.2 / NA 42intergenic (+396/‑1387)leuX/insC6tRNA‑Leu/IS2 insertion element repressor InsA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (18/23);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAACTCGAGATTTATTAGCACCTATTAAGGCCGT  >  NC_000913/4496638‑4497129
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ggTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCAAACGGATAAGCAATATTAGCTGACTGATGGTGGCTGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGTGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGacgc                                                                                                                                                                                                                                                   <  1:371423/251‑1 (MQ=255)
      aTAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTAAGTCGTGACGAGACGCGTAAgc                                                                                                                                                                                                                                             >  2:423078/1‑251 (MQ=255)
               gTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCg                                                                                                                                                                                                                                     <  1:249466/250‑1 (MQ=255)
               gTACCAACAGGGAGGGAAAACGATGGAATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCAGCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCg                                                                                                                                                                                                                                     <  2:62769/250‑1 (MQ=255)
                 aCCAACAGGGAGGGAATACGATGGCATTAACAGATATCAACGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCgcat                                                                                                                                                                                                                                  >  2:287684/1‑248 (MQ=255)
                   cAACAGGGAGGGAATACGATTGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATgg                                                                                                                                                                                                                                <  1:321602/251‑1 (MQ=255)
                       agggagggAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTc                                                                                                                                                                                                                            <  1:468234/251‑1 (MQ=255)
                               aaTACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCggg                                                                                                                                                                                                                    <  2:195406/251‑1 (MQ=255)
                                  aCGATGGCATTAACAGATAGCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGgag                                                                                                                                                                                                                  <  2:413761/250‑1 (MQ=255)
                                     aTGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACaa                                                                                                                                                                                                              <  1:122335/251‑1 (MQ=255)
                                      tGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGAATGCATATGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGAGTCGATCCGGGAGACaa                                                                                                                                                                                                              >  1:25425/1‑250 (MQ=255)
                                       ggCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACaaaa                                                                                                                                                                                                            <  2:153749/251‑1 (MQ=255)
                                               cAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCGGACTGATGGTGGCGGGATGCATCTGCTTGTCCATCAAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAAt                                                                                                                                                                                                    <  2:376752/251‑1 (MQ=255)
                                                               agcagcCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACaga                                                                                                                                                                                    <  2:82003/251‑1 (MQ=255)
                                                                 cagcCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGt                                                                                                                                                                                  >  1:90133/1‑251 (MQ=255)
                                                                  agccaagccaACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGAAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGCCGAGACGCGGAAGCTGCTTGCGAATGGCGTCGATCCGGGGGACAAAAAGAACAATGATAAGGGTGAACagaga                                                                                                                                                                                  >  1:445347/1‑249 (MQ=255)
                                                                                 tAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCttt                                                                                                                                                                   >  1:263150/1‑250 (MQ=255)
                                                                                      aaTATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGCGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGGGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAaaga                                                                                                                                                              >  1:447387/1‑250 (MQ=255)
                                                                                       atatTCGCTGACTGATGGTGGCGGCATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCCAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGCGGCGAGTAAACTGCTTGCGAATGGCGTCGATACGGGGGACAAAAAGAAAAATGATAAGGTTGAACAGTGTAAAGCACGAACCTTTAaaga                                                                                                                                                              >  2:469579/1‑249 (MQ=255)
                                                                                           tAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCAACCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGTTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACAAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTcg                                                                                                                                                         <  1:236897/250‑1 (MQ=255)
                                                                                                 actgaTGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGCCAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTg                                                                                                                                                   >  1:226547/1‑250 (MQ=255)
                                                                                                 actgaTGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGa                                                                                                                                                  <  1:423078/251‑1 (MQ=255)
                                                                                                   tgaTGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAgg                                                                                                                                                >  1:258254/1‑250 (MQ=255)
                                                                                                      tggtggCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGGACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGTGACAAAAAGAAAAATGATAAGGTTGAACAGAGAAAAGCCCGAACCTTTAAAGAAGTAGAGATTGAGtggc                                                                                                                                             >  2:465351/1‑251 (MQ=255)
                                                                                                             ggTATGCATCTGCTGTTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACTCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGATTGGCATGGCAc                                                                                                                                      <  2:438190/251‑1 (MQ=255)
                                                                                                                                              aaGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAAGGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCAt                                                                                                       <  2:290238/249‑1 (MQ=255)
                                                                                                                                                   cTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGc                                                                                                >  1:415034/1‑251 (MQ=255)
                                                                                                                                                        gTTTGCAGTACCGTTATGAGGGAAAGCAACAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTaaa                                                                                            >  1:353066/1‑250 (MQ=255)
                                                                                                                                                           tGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAg                                                                                         <  2:472536/250‑1 (MQ=255)
                                                                                                                                                                     ttATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGAt                                                                               >  1:351428/1‑250 (MQ=255)
                                                                                                                                                                           gggAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGaaa                                                                                                                                                                                                       <  1:434563/124‑1 (MQ=255)
                                                                                                                                                                           gggAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGaaa                                                                                                                                                                                                       >  2:434563/1‑124 (MQ=255)
                                                                                                                                                                                          gCTGGCACTTGGGGTTTATCCTGAAATCACACGAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATTGCGTCGATCCGTGAGACAAAAAGAAAAATTATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCTATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTTCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGt                                                          <  1:463601/250‑1 (MQ=255)
                                                                                                                                                                                                  ttGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGACGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTATTGAAGATAATCTTTTTGCAGCGCTAGGTGAACGTAat                                                 >  1:237891/1‑251 (MQ=255)
                                                                                                                                                                                                   tGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTaa                                                  <  2:263150/249‑1 (MQ=255)
                                                                                                                                                                                                                   aTCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAAc                                >  1:496715/1‑251 (MQ=255)
                                                                                                                                                                                                                      acacTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAGTGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAACTCg                             <  1:287684/251‑1 (MQ=255)
                                                                                                                                                                                                                      acacTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAAATCg                             >  1:27578/1‑251 (MQ=255)
                                                                                                                                                                                                                       cacTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGTTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAACTCga                            <  1:357004/251‑1 (MQ=255)
                                                                                                                                                                                                                        actagcgtagacacatggaCGTCGTGACGAGACGGGTAAGCTTCTTGCGAATGGAGTCGATCCGGGAGACAAAAATAAAAATGATAAGGATGAACAGACGAAAGCACGAACCTTTAAAGAAGTCTCGATTGAGTGGCATGGCACCAATAAAAAGTGTTCTTAAGATCACGCTCACCGTGGTCTAAAAATCCTTGAAGATAATCTTCTTGCAGCGCTAGGTGAACGTAATAGCCCGTATTTAAAAACTCgag                           <  2:455799/233‑1 (MQ=255)
                                                                                                                                                                                                                           aGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAACTCGAGAt                         <  1:282413/250‑1 (MQ=255)
                                                                                                                                                                                                                                     ggagTACGTCGTGACGAGACTAGTACGCTGCTTGCGAATGGCGTCGACCCGGGAGACAAAAAGAAAAATGATAAGGTTTATCAGAGTAAAGCACGAACCTGTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAATATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAACCCGAGATTTATTAGCAcc              <  2:351852/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   cGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGACTTGAAGCTAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAACTCGAGATTTATTAGCACATATTAAGGCCga  >  2:323951/1‑248 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTAAAAACTCGAGATTTATTAGCACCTATTAAGGCCGT  >  NC_000913/4496638‑4497129

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