Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,618,489 C→T T87T (ACC→ACT deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,618,4890CT98.1% 162.7 / NA 52T87T (ACC→ACTdeoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/24);  minor base G (0/1);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCG  >  NC_000913/4618244‑4618735
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
gAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACtggtgg                                                                                                                                                                                                                                                        <  1:311109/251‑1 (MQ=255)
        tCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCg                                                                                                                                                                                                                                                >  1:409343/1‑251 (MQ=255)
            aaaaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGACTCAGGAACCGTTCTCGCCTGGAAAAGCCTGCATCTGAATGGCCCGCTTGTTGATAAACATTCCACTGGTGGCGTGGGCGAtgg                                                                                                                                                                                                                                            >  1:379735/1‑250 (MQ=255)
              aaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGAtgtg                                                                                                                                                                                                                                           >  1:290734/1‑250 (MQ=255)
                           cATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGtt                                                                                                                                                                                                                             >  1:193234/1‑251 (MQ=255)
                              gctcggAGCGATGAAGAAATTCGTTTCTTTATCAACGGCACTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGGTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTgg                                                                                                                                                                                                                           <  2:346753/245‑1 (MQ=255)
                                  tGAGCGCTGAAGAAATTAGTTACTTTATCTACGGCATTAGCGACAACACTATCTCAGAAGGGCAGATTGCCGCCCTCGCGATGCCCATTTTCTTACACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCTATGCGAGCTTCAGGAACCGTTCTCGACTGGAAAAGCCTGAATCTGAATGGCCCGATTGTTGATAAACATTACCCTGGAGGCGTGGGCGATATGACTTCGCAGTTGTTggagc                                                                                                                                                                                                                        >  2:284241/1‑246 (MQ=255)
                                            agaaATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAg                                                                                                                                                                                                            <  1:416897/251‑1 (MQ=255)
                                              aaaTTCGTTTCTTTAACAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGcc                                                                                                                                                                                                          <  1:433874/251‑1 (MQ=255)
                                                      ttCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGCGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGcggcgg                                                                                                                                                                                                  >  2:424468/1‑251 (MQ=255)
                                                       tCTTTATCAACGGCATTCGCGACAACACTATCTCCGACGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGCTAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCggcggc                                                                                                                                                                                                 >  2:66913/1‑251 (MQ=255)
                                                         tttATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGcggcgggc                                                                                                                                                                                                >  2:365809/1‑248 (MQ=255)
                                                               aaCGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGCCTTCGATGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTcc                                                                                                                                                                                         >  1:468739/1‑251 (MQ=255)
                                                               aaCGGCATTCGCGAAAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTcc                                                                                                                                                                                         <  1:246726/251‑1 (MQ=255)
                                                                           gACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGTTCTCTgg                                                                                                                                                                             >  1:389216/1‑251 (MQ=255)
                                                                            acaacaCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGGTGGTCGCAGCCTGCGGCGGCCATATTCCgcttgtttctggg                                                                                                                                                                            >  2:51560/1‑239 (MQ=255)
                                                                                cacTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTcgcg                                                                                                                                                                        >  1:390536/1‑251 (MQ=255)
                                                                                        ccGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTc                                                                                                                                                                   <  1:418127/248‑1 (MQ=255)
                                                                                         cGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGACCCGTTCTCGCCTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGACGCAGGCTGCGGCGGCTATATTCCGATGGTCACTGGTCGCGGCCCCGGTc                                                                                                                                                               >  2:279764/1‑251 (MQ=255)
                                                                                         cGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCAGATGATCGCTGGGCGCGGCCTCGGt                                                                                                                                                                >  2:406816/1‑250 (MQ=255)
                                                                                               gCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTg                                                                                                                                                         <  2:193234/251‑1 (MQ=255)
                                                                                                                cgATGACCATTTTCTTCCCCGATATGACAATGCCTGAGCGGGTCTCGCTGACCAGGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACaaa                                                                                                                                       <  2:390536/251‑1 (MQ=255)
                                                                                                                 gATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGCTATCTGGACGCGACCTCCGTCATACTGGCGGTACGCTCGACGAAc                                                                                                                                      >  1:59595/1‑251 (MQ=255)
                                                                                                                 gATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCAACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGACGCCTGCGGAGGCCATATTCCGATGAGCTTTGGTCGCGGCCTCAGTTCTTCTGGCGGCTCGCTCGAccca                                                                                                                                       >  2:331936/1‑247 (MQ=255)
                                                                                                                    gACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGCTGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGTCCTCGGCCATCCAGGCGGTACGccggacacaccgg                                                                                                                                   >  2:368634/1‑239 (MQ=255)
                                                                                                                         ttttCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAatc                                                                                                                               <  1:292293/250‑1 (MQ=255)
                                                                                                                               tcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAatccatcc                                                                                                                         <  1:147167/249‑1 (MQ=255)
                                                                                                                               tcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAatccatc                                                                                                                          <  1:376472/248‑1 (MQ=255)
                                                                                                                               tcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATccc                                                                                                                        <  1:66913/250‑1 (MQ=255)
                                                                                                                                ccACGATATGGCAATGCCTGAGCGTGTCTCGCTGACCATGTCGAGGCGTGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAAGGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCt                                                                                                                       <  2:381645/250‑1 (MQ=255)
                                                                                                                                 cACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAACGGAATTCATCCCTgg                                                                                                                     >  1:207518/1‑251 (MQ=255)
                                                                                                                                     agatgccAATGCCTGGGCGTGTCTCGCAGACCATGGCGATGCGTGATGCAGGAACCGTTCTCGACCGGAAAAGCCTGCATCGTAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCTTCGGCGATGTGGCTTCGCTGATGTTGGGTCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGGCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTc                                                                                                                 <  2:190817/245‑1 (MQ=255)
                                                                                                                                     aaatGACAATGCCGGAGCGTGGCTCGCTTACCATGGCGATGCGAGAGTCAGGAACCGTTCTCGACGTGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCTCTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGAGGATCTCTGGTCGCTGCCGCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTc                                                                                                                 <  2:480424/249‑1 (MQ=255)
                                                                                                                                       atGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGCCTGGACAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGTTGATCTCTTGGCGGCGCCTCGGTCCTACTGGCGGATCCGTCGCCCAACTGGAAtccctccccgggtttac                                                                                                               >  2:103966/1‑242 (MQ=255)
                                                                                                                                           cAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCtata                                                                                                                                                                                             >  2:75763/1‑171 (MQ=255)
                                                                                                                                           cAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCtata                                                                                                                                                                                             <  1:75763/171‑1 (MQ=255)
                                                                                                                                                  tGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGCAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGGTGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGGACGCTCGACaaaa                                                                                                                                      >  1:410522/1‑217 (MQ=255)
                                                                                                                                                  tGAGCGTGTCTCGCTGACCATGGCGATGCGAGATCCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAc                                                                                                                                      <  2:410522/218‑1 (MQ=255)
                                                                                                                                                            ccGCTGACCATGGCGATGTGAGAGTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCg                                                                                        <  1:260854/250‑1 (MQ=255)
                                                                                                                                                             cGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAAcc                                                                                         <  2:379233/249‑1 (MQ=255)
                                                                                                                                                             cGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGGCATACTGGCGGTACGCTCGGCAAACTGGGATCCATCCCTGGCTTCGGCAATTTCCCGGATGACaaacgg                                                                                       >  1:200148/1‑247 (MQ=255)
                                                                                                                                                                    cATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTg                                                                                                                                                                              <  2:22990/161‑1 (MQ=255)
                                                                                                                                                                    cATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTg                                                                                                                                                                              >  1:22990/1‑161 (MQ=255)
                                                                                                                                                                    cATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGAc                                                                                                              <  2:134497/223‑1 (MQ=255)
                                                                                                                                                                                               tctcGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGGTCTCTGGTCGCGGCCTCGGTCATACTGGCGGGACGCTCGGCAAACTGGAATCCATCCCCGGCTTCGCCCTTTTTCCCGGGGGCAAACGTTTCCCCGAAAATTTTAAAAACATCGGCggggggg                                                     >  2:454073/1‑245 (MQ=255)
                                                                                                                                                                                                  cGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGGCGCGGCCTCGGTCATACTGGCGGTAATCCCGACAAAATGGAATCCCTTCCCGGGCTCGACATTTTCCCCGGTGACAACCGTTTCCCCGGAATTTTTAAAGACGTCGGCGTgggggttt                                                  >  2:169502/1‑245 (MQ=255)
                                                                                                                                                                                                       gaaaaatcttgcttctaaaTGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTAACTGCTATAATTTTGGGGCCTATGGTCGCAGCCTGCGGCGGCTATATTTCGATGATCTCTGGTCGCGGCCTCGGTAATCCTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCAACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCgg                                             <  1:3408/233‑1 (MQ=255)
                                                                                                                                                                                                                 gCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAgc                                   <  1:243161/250‑2 (MQ=255)
                                                                                                                                                                                                                     cTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGt                                                                <  1:209821/218‑1 (MQ=255)
                                                                                                                                                                                                                             cccGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCgg                      <  2:344859/251‑1 (MQ=255)
                                                                                                                                                                                                                                gATTGTTGATAAACATTCCACGGGGGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGAGGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGcgg                   <  2:409343/251‑3 (MQ=255)
                                                                                                                                                                                                                                     ttgATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGGCATTTTCCCGGATGCCAACCGTTTCCGCGAAATTATTAAAAACGTCGGCGGGGCGGATTTCGGTCAGACCaccccaccggctcccgcgggtaa               >  1:274578/1‑228 (MQ=255)
                                                                                                                                                                                                                                                 ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGGTCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGGAATCCATCCCTCGCCTCGCACTTTTTTCGGGTGGCAACCCTTTTCACGAAATTTTTTAAGACGTCGGCGGGGGGGTTTTCGGTTAGACCAACTCACCGGGTCCGGCCGGGAAACGTTTTTAcgcc  >  2:162307/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
GAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCG  >  NC_000913/4618244‑4618735

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