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

GCGTTCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAGGTGGTGCAGGCTTCTGGCGGCGG  >  NC_000913/4430974‑4431449
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
ggggtCACTATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGTTTATCCTTGGCCCGTCCCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACTCGCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTggg                                                                                                                                                                                                                                   <  2:35624/247‑1 (MQ=255)
   ttCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGCTTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATATACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGCCACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGtg                                                                                                                                                                                                                                >  1:80232/1‑250 (MQ=255)
    tCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGCTACCCCGCTGTGGTTTATCGCGTTGGGGTTg                                                                                                                                                                                                                               >  1:22843/1‑251 (MQ=255)
     cACGATGATTACCACCGTTTCCGCGATTCTGGTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTg                                                                                                                                                                                                                               <  2:215907/250‑1 (MQ=255)
      aCGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCgttggggttggg                                                                                                                                                                                                                             >  1:189238/1‑251 (MQ=255)
          tgatTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctgg                                                                                                                                                                                                                         >  1:224627/1‑251 (MQ=255)
           gatTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctgg                                                                                                                                                                                                                         <  1:76303/250‑1 (MQ=255)
               ccAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGt                                                                                                                                                                                                                     <  2:188162/249‑1 (MQ=255)
               aCAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGtt                                                                                                                                                                                                                    >  1:197105/1‑251 (MQ=255)
                cAACCGTTTCCGCGATTCTGTCTATGTTCGTCTGGACGATTATCCTTTGCTCGTCCCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGttt                                                                                                                                                                                                                   <  1:385669/251‑1 (MQ=255)
                  aCCGTTTCCGCGCCTCTGTTTCTGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGCGAAGTCGATCTACAAGCTGCCGCTCGGCAAGCCGGTGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACCGTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGGTTATCGCGTTGGGGCTGGACGGGCTGTTTAt                                                                                                                                                                                                                 >  2:470970/1‑251 (MQ=255)
                            cgATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCGGGCTGTTTATTGGTaagaag                                                                                                                                                                                                       <  2:296211/251‑1 (MQ=255)
                              aTTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCg                                                                                                                                                                                                     <  1:100466/251‑1 (MQ=255)
                                   gTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTTCCTTGTGTATCGCAAAGAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTTCTGGGTATGTATGGCGTGCTTTGTGTGCGTGGTCGTGTTGCTGACACTTGAAGATGACACTCGCCAGGCGCTGCTGGTTGCCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGgct                                                                                                                                                                                                 <  1:328828/250‑1 (MQ=255)
                                      tATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGgctgct                                                                                                                                                                                              <  2:413965/250‑1 (MQ=255)
                                      tATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGATCGGCAAGCTCATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGAGGgctgctt                                                                                                                                                                                             >  2:398766/1‑250 (MQ=255)
                                         gTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAAc                                                                                                                                                                                          >  1:66780/1‑251 (MQ=255)
                                            cGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTg                                                                                                                                                                                        <  1:233883/250‑1 (MQ=255)
                                            cGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGAGCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTg                                                                                                                                                                                        <  2:239550/250‑1 (MQ=255)
                                                 ggACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGCAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGCCACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGCACTGCGGaa                                                                                                                                                                                   >  2:86911/1‑250 (MQ=255)
                                                  gACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGCGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGCTGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGAGGGCTGCTGAACTGCGGAAAt                                                                                                                                                                                 >  1:414013/1‑251 (MQ=255)
                                                     gATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAAc                                                                                                                                                                              >  1:295344/1‑251 (MQ=255)
                                                      aTTATCCTTTGCGCGTACCTTGTGTATCGCAAACAGCGTCCTCAGCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACGGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTGTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAAcc                                                                                                                                                                             <  1:171858/251‑1 (MQ=255)
                                                               tGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCAGCTACATGAGAAGTCGATCTACAAGATGCCGTTCGGCAAGCTGAGGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTTGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATc                                                                                                                                                                     <  2:486213/250‑1 (MQ=255)
                                                                gCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCAt                                                                                                                                                                   <  2:224627/251‑1 (MQ=255)
                                                                gCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCaa                                                                                                                                                                   <  1:236325/251‑2 (MQ=255)
                                                                   cGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCt                                                                                                                                                                <  1:470392/251‑1 (MQ=255)
                                                                        cTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGAATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACAATCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTGTTGGTAAGAAGCGGGCTGCTGCACTGCGGGAATAACCGCATTAt                                                                                                                                                                      >  2:486671/1‑240 (MQ=255)
                                                                        cTTGTGTATAGCAACCAGCGTCCTCATCTACATGCGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTCTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTAt                                                                                                                                                                      <  1:486671/240‑1 (MQ=255)
                                                                             gtATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCa                                                                                                                                                      <  1:51012/251‑1 (MQ=255)
                                                                                 cGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGc                                                                                                                                                  >  1:245745/1‑251 (MQ=255)
                                                                                   cAAACAGCGTCCGCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCa                                                                                                                                                <  2:197105/251‑1 (MQ=255)
                                                                                    aaaCAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCAt                                                                                                                                               >  1:436383/1‑251 (MQ=255)
                                                                                      aCAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCACGCTGATGTGCTGGGAATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAATCGGTCGGCTGAACTGCGTATATAACCGCGTTATCATGCTGGATAGCGCAATGCCaacc                                                                                                                                             >  2:110639/1‑248 (MQ=255)
                                                                                      aCAGAGTCCTCATCTCCATGGGACGTCGATCTTCAAGATGCAGCTCGTCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACCGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGGTTAGCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCCTTATCATGCTGTAGGGCGCAATGCCATcc                                                                                                                                             <  2:76058/251‑1 (MQ=255)
                                                                                                cATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTa                                                                                                                                     <  1:497780/249‑1 (MQ=255)
                                                                                                     aCATTAGAAGTCGATCTTCAAGTTGCCGCTCGTCAAGCTGATGTGTTAGGTATGTATGGCGTACTTTGTGTGCGTGGTCGTGTTGCTGACACTTGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCACGTTGGGGTTGGGCTGGCTTTTTATTGGTAAGAGGCGGGCTGCTTAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGAtcac                                                                                                                               <  1:375939/250‑1 (MQ=255)
                                                                                                     aCATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACGCTGGAAGATGACACTCGCCAGGCGCTGCTGGTGACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCGGCTGTTAGAtc                                                                                                                                 <  2:28999/248‑1 (MQ=255)
                                                                                                       aTGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCCTCCAGCTTTTAGAtcactc                                                                                                                             >  1:379429/1‑250 (MQ=255)
                                                                                                          agaAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCTTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCAccc                                                                                                                         <  1:302488/251‑1 (MQ=255)
                                                                                                           gaAGTCGATCTACAAGATGGCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCg                                                                                                                        <  2:309973/251‑1 (MQ=255)
                                                                                                              gTCGATCTACAAGATGCCGCTCGGCAAGCTGGTGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctggct                                                                                                                                                                                                                       >  2:466617/1‑153 (MQ=255)
                                                                                                              gTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGTTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGggctggct                                                                                                                                                                                                                       <  1:466617/153‑1 (MQ=255)
                                                                                                                     tACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACAATCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGcgcgcg                                                                                                              >  2:419019/1‑251 (MQ=255)
                                                                                                                      aCAAGATGCCGCGCGGCAAGCGTATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTTTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAAGGCCATCCAGCTTTTAGATCACTCCCCCGCCAGCGCGCGt                                                                                                             <  1:233777/251‑1 (MQ=255)
                                                                                                                        aaGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCCGCTTTTAGATCCCTCACCCGCCAGCGCGCGTgg                                                                                                           >  2:65097/1‑251 (MQ=255)
                                                                                                                              ccGCTCGGCAGGCGGATGTGCTGGGTATGTAGGGTGTTCTTTGTGTTCGTTGGCGGTGGGCTGACCCGGGAAGATGCCCCGCGCCAGGCGCTGCTGGTTCCCCCGCGGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGAGGGCGCAATGCCATCCAGCTTTTAGATCCCTCCCCCGCCAGCGCGCGTGGGAACAg                                                                                                     <  2:414013/251‑1 (MQ=255)
                                                                                                                                   cGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGGCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACa                                                                                                 <  1:368310/250‑1 (MQ=255)
                                                                                                                                            gatggggtGGGTATGTATGGCTTTCTTTGTGTTGGTGGTCGTGTTGCTGACATTGGAAGATGCCACTCGCCAGGCGCTGCTGTTTACCCCGCTGTGGGTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGTTTTTAGATCACACCCCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCa                                                                                       <  2:83078/244‑1 (MQ=255)
                                                                                                                                               gtgCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGGACAGTACATTGTTCTCCAGAc                                                                                    >  2:501096/1‑251 (MQ=255)
                                                                                                                                                        atgtatgGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGCACAGTACATTGTTTTCCAGACTGATGTGAt                                                                           >  2:437686/1‑251 (MQ=255)
                                                                                                                                                          gtatgGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGGTGACACTCGCCAGGCGCTGCTGGTTACCCCGCTTTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGGACTGCGGAAATAACCGCATTATCATGCTGGCTGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCACCAGCGCGCGTGGGCACAGTACATGGTTTTCCAGCCTgctttgatac                                                                         >  2:270799/1‑242 (MQ=255)
                                                                                                                                                                  tttttgttgtTCGTGGTCGTGTGTCTGACACTGGAAGATGACACTCTCCGGGCGGTGTTGGTTACCCCGCTGTGGTTTAGCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAGGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCGGGTTGGCGCAATGCCATCCAGTTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGt                                                                  <  2:45250/243‑1 (MQ=255)
                                                                                                                                                                   tCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGATTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtc                                                                 >  1:348500/1‑250 (MQ=255)
                                                                                                                                                                     tttgtgtgtgggggcgttttGCTGACAGTGGAAGATGACACTCGCCGGGCGCGGCTGTTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGACCTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcat                                                               <  2:436383/232‑1 (MQ=255)
                                                                                                                                                                     tttGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc                                                              >  1:217679/1‑251 (MQ=255)
                                                                                                                                                                     tttGTGTTCGTGGTCGTGGGGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGGGGTTTATCGCGTTGGGTTTGGGCGGGCTGTTTATGGGTAAGAAGGGGGCTGCTGAACTGCGGAAATAACCGCATTTTCATGTTGGATGGCGCAATGCCATCCGGCTTTGAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcat                                                               <  1:251306/250‑1 (MQ=255)
                                                                                                                                                                          tgtcgtgggggtgtTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGGTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCatcaatc                                                          <  2:189238/240‑1 (MQ=255)
                                                                                                                                                                                     tgtTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATg                                              >  1:372905/1‑251 (MQ=255)
                                                                                                                                                                                         gCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGTTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCGGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGt                                          <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
GCGTTCACGATGATTACAACCGTTTCCGCGATTCTGTTTATGTTCGTCTGGACGATTATCCTTTGCTCGTACCTTGTGTATCGCAAACAGCGTCCTCATCTACATGAGAAGTCGATCTACAAGATGCCGCTCGGCAAGCTGATGTGCTGGGTATGTATGGCGTTCTTTGTGTTCGTGGTCGTGTTGCTGACACTGGAAGATGACACTCGCCAGGCGCTGCTGGTTACCCCGCTGTGGTTTATCGCGCTGGGGTTGGGCTGGCTGTTTATTGGTAAGAAGCGGGCTGCTGAACTGCGGAAATAACCGCATTATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATTAATGCCGTTATACATCGCTTTCCAGGTGGTGCAGGCTTCTGGCGGCGG  >  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: