Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,431,220 C→T L453L (CTG→TTG)  cycA → D‑serine/alanine/glycine/:H(+)symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,431,2200CT95.3% 184.5 / ‑0.4 65L453L (CTG→TTG) cycAD‑serine/alanine/glycine/:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (27/34);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

GCGTTCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAGGTGGTGCAGGCTTCTGGCGGCGG  >  NC_000913/4430974‑4431449
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
ggggtCACTATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGTTTATCCTTGGCCCGTCCCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACTCGCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTggg                                                                                                                                                                                                                                   <  2:35624/247‑1 (MQ=255)
   ttCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGCTTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATATACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGCCACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGtg                                                                                                                                                                                                                                >  1:80232/1‑250 (MQ=255)
    tCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGCTACCCCGCTGTGGTTTATCGCGTTGGGGTTg                                                                                                                                                                                                                               >  1:22843/1‑251 (MQ=255)
     cACGATGATTACCACCGTTTCCGCGATTCTGGTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTg                                                                                                                                                                                                                               <  2:215907/250‑1 (MQ=255)
      aCGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCgttggggttggg                                                                                                                                                                                                                             >  1:189238/1‑251 (MQ=255)
          tgatTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctgg                                                                                                                                                                                                                         >  1:224627/1‑251 (MQ=255)
           gatTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctgg                                                                                                                                                                                                                         <  1:76303/250‑1 (MQ=255)
               ccAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGt                                                                                                                                                                                                                     <  2:188162/249‑1 (MQ=255)
               aCAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGtt                                                                                                                                                                                                                    >  1:197105/1‑251 (MQ=255)
                cAACCGTTTCCGCGATTCTGTCTATGTTCGTCTGGACGATTATCCTTTGCTCGTCCCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGttt                                                                                                                                                                                                                   <  1:385669/251‑1 (MQ=255)
                  aCCGTTTCCGCGCCTCTGTTTCTGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGCGAAGTCGATCTACAAGCTGCCGCTCGGCAAGCCGGTGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACCGTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGGTTATCGCGTTGGGGCTGGACGGGCTGTTTAt                                                                                                                                                                                                                 >  2:470970/1‑251 (MQ=255)
                            cgATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCGGGCTGTTTATTGGTaagaag                                                                                                                                                                                                       <  2:296211/251‑1 (MQ=255)
                              aTTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCg                                                                                                                                                                                                     <  1:100466/251‑1 (MQ=255)
                                   gTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTTCCTTGTGTATCGCAAAGAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTTCTGGGTATGTATGGCGTGCTTTGTGTGCGTGGTCGTGTTG‑CTGACACTTGAAGATGACACTCGCCAGGCGCTGCTGGTTGCCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGgct                                                                                                                                                                                                 <  1:328828/250‑1 (MQ=255)
                                      tATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGgctgct                                                                                                                                                                                              <  2:413965/250‑1 (MQ=255)
                                      tATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGATCGGCAAGCTCATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGAGGgctgctt                                                                                                                                                                                             >  2:398766/1‑250 (MQ=255)
                                         gTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAAc                                                                                                                                                                                          >  1:66780/1‑251 (MQ=255)
                                            cGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTg                                                                                                                                                                                        <  1:233883/250‑1 (MQ=255)
                                            cGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGAGCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTg                                                                                                                                                                                        <  2:239550/250‑1 (MQ=255)
                                                 ggACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGCAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGCCACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGCACTGCGGaa                                                                                                                                                                                   >  2:86911/1‑250 (MQ=255)
                                                  gACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGCGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGCTGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGAGGGCTGCTGAACTGCGGAAAt                                                                                                                                                                                 >  1:414013/1‑251 (MQ=255)
                                                     gATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAAc                                                                                                                                                                              >  1:295344/1‑251 (MQ=255)
                                                      aTTATCCTTTGCGCGTACCTTGTGTATCGCAAACAGCGTCCTCAGCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACGGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTGTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAAcc                                                                                                                                                                             <  1:171858/251‑1 (MQ=255)
                                                               tGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCAGCTACATGAGAAGTCGATCTACAAGATGCCGTTCGGCAAGCTGAGGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTTGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATc                                                                                                                                                                     <  2:486213/250‑1 (MQ=255)
                                                                gCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCAt                                                                                                                                                                   <  2:224627/251‑1 (MQ=255)
                                                                gCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCaa                                                                                                                                                                   <  1:236325/251‑2 (MQ=255)
                                                                   cGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCt                                                                                                                                                                <  1:470392/251‑1 (MQ=255)
                                                                        cTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGAATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACAATCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTGTTGGTAAGAAGCGGGCTGCTGCACTGCGGGAATAACCGCATTAt                                                                                                                                                                      >  2:486671/1‑240 (MQ=255)
                                                                        cTTGTGTATAGCAACCAGCGTCCTCATCTACATGCGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTCTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTAt                                                                                                                                                                      <  1:486671/240‑1 (MQ=255)
                                                                             gtATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCa                                                                                                                                                      <  1:51012/251‑1 (MQ=255)
                                                                                 cGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGc                                                                                                                                                  >  1:245745/1‑251 (MQ=255)
                                                                                   cAAACAGCGTCCGCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCa                                                                                                                                                <  2:197105/251‑1 (MQ=255)
                                                                                    aaaCAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCAt                                                                                                                                               >  1:436383/1‑251 (MQ=255)
                                                                                      aCAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCACGCTGATGTGCTGGGAATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAATCGGTCGGCTGAACTGCGTATATAACCGCGTTATCATGCTGGATAGCGCAATGCCaacc                                                                                                                                             >  2:110639/1‑248 (MQ=255)
                                                                                      aCAGAGTCCTCATCTCCATGGGACGTCGATCTTCAAGATGCAGCTCGTCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACCGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGGTTAGCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCCTTATCATGCTGTAGGGCGCAATGCCATcc                                                                                                                                             <  2:76058/251‑1 (MQ=255)
                                                                                                cATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTa                                                                                                                                     <  1:497780/249‑1 (MQ=255)
                                                                                                     aCATTAGAAGTCGATCTTCAAGTTGCCGCTCGTCAAGCTGATGTGTTAGGTATGTATGGCGTACTTTGTGTGCGTGGTCGTGTTG‑CTGACACTTGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCACGTTGGGGTTGGGCTGGCTTTTTATTGGTAAGAGGCGGGCTGCTTAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGAtcac                                                                                                                               <  1:375939/250‑1 (MQ=255)
                                                                                                     aCATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACGCTGGAAGATGACACTCGCCAGGCGCTGCTGGTGACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCGGCTGTTAGAtc                                                                                                                                 <  2:28999/248‑1 (MQ=255)
                                                                                                       aTGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCCTCCAGCTTTTAGAtcactc                                                                                                                             >  1:379429/1‑250 (MQ=255)
                                                                                                          agaAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCTTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCAccc                                                                                                                         <  1:302488/251‑1 (MQ=255)
                                                                                                           gaAGTCGATCTACAAGATGGCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCg                                                                                                                        <  2:309973/251‑1 (MQ=255)
                                                                                                              gTCGATCTACAAGATGCCGCTCGGCAAGCTGGTGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctggct                                                                                                                                                                                                                       >  2:466617/1‑153 (MQ=255)
                                                                                                              gTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGTTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctggct                                                                                                                                                                                                                       <  1:466617/153‑1 (MQ=255)
                                                                                                                     tACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACAATCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGcgcgcg                                                                                                              >  2:419019/1‑251 (MQ=255)
                                                                                                                      aCAAGATGCCGCGCGGCAAGCGTATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTTTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAAGGCCATCCAGCTTTTAGATCACTCCCCCGCCAGCGCGCGt                                                                                                             <  1:233777/251‑1 (MQ=255)
                                                                                                                        aaGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCCGCTTTTAGATCCCTCACCCGCCAGCGCGCGTgg                                                                                                           >  2:65097/1‑251 (MQ=255)
                                                                                                                              ccGCTCGGCAGGCGGATGTGCTGGGTATGTAGGGTGTTCTTTGTGTTCGTTGGCG‑GTGGGCTGACCCGGGAAGATGCCCCGCGCCAGGCGCTGCTGGTTCCCCCGCGGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGAGGGCGCAATGCCATCCAGCTTTTAGATCCCTCCCCCGCCAGCGCGCGTGGGAACAg                                                                                                     <  2:414013/251‑1 (MQ=255)
                                                                                                                                   cGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGGCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACa                                                                                                 <  1:368310/250‑1 (MQ=255)
                                                                                                                                            gatggggtGGGTATGTATGGCTTTCTTTGTGTTGGTGGTCGTGTTG‑CTGACATTGGAAGATGCCACTCGCCAGGCGCTGCTGTTTACCCCGCTGTGGGTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGTTTTTAGATCACACCCCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCa                                                                                       <  2:83078/244‑1 (MQ=255)
                                                                                                                                               gtgCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGGACAGTACATTGTTCTCCAGAc                                                                                    >  2:501096/1‑251 (MQ=255)
                                                                                                                                                        atgtatgGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGCACAGTACATTGTTTTCCAGACTGATGTGAt                                                                           >  2:437686/1‑251 (MQ=255)
                                                                                                                                                          gtatgGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGGTGACACTCGCCAGGCGCTGCTGGTTACCCCGCTTTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGGACTGCGGAAATAACCGCATTATCATGCTGGCTGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCACCAGCGCGCGTGGGCACAGTACATGGTTTTCCAGCCTgctttgatac                                                                         >  2:270799/1‑242 (MQ=255)
                                                                                                                                                                  tttttgttgtTCGTGGTCGTGTGT‑CTGACACTGGAAGATGACACTCTCCGGGCGGTGTTGGTTACCCCGCTGTGGTTTAGCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAGGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCGGGTTGGCGCAATGCCATCCAGTTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGt                                                                  <  2:45250/243‑1 (MQ=255)
                                                                                                                                                                   tCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGATTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtc                                                                 >  1:348500/1‑250 (MQ=255)
                                                                                                                                                                     tttgtgtgtgggggcgttttG‑CTGACAGTGGAAGATGACACTCGCCGGGCGCGGCTGTTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGACCTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcat                                                               <  2:436383/232‑1 (MQ=255)
                                                                                                                                                                     tttGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc                                                              >  1:217679/1‑251 (MQ=255)
                                                                                                                                                                     tttGTGTTCGTGGTCGTGGGG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGGGGTTTATCGCGTTGGGTTTGGGCGGGCTGTTTATGGGTAAGAAGGGGGCTGCTGAACTGCGGAAATAACCGCATTTTCATGTTGGATGGCGCAATGCCATCCGGCTTTGAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcat                                                               <  1:251306/250‑1 (MQ=255)
                                                                                                                                                                          tgtcgtgggggtgtTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGGTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCatcaatc                                                          <  2:189238/240‑1 (MQ=255)
                                                                                                                                                                                     tgtTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATg                                              >  1:372905/1‑251 (MQ=255)
                                                                                                                                                                                         g‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCGGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGt                                          <  1:462512/251‑1 (MQ=255)
                                                                                                                                                                                                        aTGACACTCGCCAGGCGCTGCTGGTTCCCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTc                            <  1:86911/251‑1 (MQ=255)
                                                                                                                                                                                                              cTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAggtg                      >  1:277337/1‑251 (MQ=255)
                                                                                                                                                                                                               ccGCCAGGCCCGGCTGGTTACCCCGCTGTGGTTTAGCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCGGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGTTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAggt                       <  2:134822/248‑1 (MQ=255)
                                                                                                                                                                                                                 gCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGGACAGTACATTGTTTTCCAGCCTGATGTGATCCATCAGGTAATCAATAAGTTCATTAATGCCGTTATACATCGCTTTCCAGggggtg                   >  2:440494/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGGCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAGGTGGTGc                  >  2:391145/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cccccGCTGTGGTTTAGCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAGGTGGTGCAGGCTTCTGGcggcgg  <  1:255845/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
GCGTTCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTG‑CTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAGGTGGTGCAGGCTTCTGGCGGCGG  >  NC_000913/4430974‑4431449

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