Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,430,584 C→T P241S (CCA→TCA)  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,430,5840CT92.4% 176.8 / 3.6 66P241S (CCA→TCA) cycAD‑serine/alanine/glycine/:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (32/29);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

GAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGT  >  NC_000913/4430338‑4430830
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gagaTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCaga                                                                                                                                                                                                                                                    <  1:282623/251‑1 (MQ=255)
   aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGAGTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGATGGGAGGTCCAAAAGCTGCGGAAACCAAAGATTCaggga                                                                                                                                                                                                                                                  >  2:319636/1‑247 (MQ=255)
      gAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGGGAAATc                                                                                                                                                                                                                                              >  2:192912/1‑251 (MQ=255)
      gAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATc                                                                                                                                                                                                                                              <  1:319636/251‑1 (MQ=255)
             ggTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGCGAAATCACTGcc                                                                                                                                                                                                                                        >  1:661700/1‑250 (MQ=255)
             ggTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTg                                                                                                                                                                                                                                          <  1:149290/248‑1 (MQ=255)
                   cGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCAcgcgcg                                                                                                                                                                                                                                 >  1:505861/1‑251 (MQ=255)
                       gatCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGTTTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAATTCACTGCCACGCGCGAtt                                                                                                                                                                                                                              >  2:238626/1‑250 (MQ=255)
                             aaTCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACCACAGCTGCGGGAACCACAAATTCAAAGAAATCACTGCCACGCGCGATTAACTcac                                                                                                                                                                                                                       >  2:650042/1‑249 (MQ=255)
                              aTCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGAGATTAACTCCAt                                                                                                                                                                                                                      >  1:592197/1‑251 (MQ=255)
                                cgtcgtCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCAtt                                                                                                                                                                                                                     <  1:501133/250‑1 (MQ=255)
                                      cATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATcc                                                                                                                                                                                                              >  1:120951/1‑251 (MQ=255)
                                       aTTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCg                                                                                                                                                                                                             >  1:464379/1‑251 (MQ=255)
                                              cGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCAtt                                                                                                                                                                                                      <  1:177720/251‑1 (MQ=255)
                                               gCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCAttt                                                                                                                                                                                                     >  1:483490/1‑250 (MQ=255)
                                               gCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTa                                                                                                                                                                                                     <  1:82706/251‑1 (MQ=255)
                                                 tGATTGTCGTCGGCCTGGTCATGGTGGCGATGCAATTTCAGTCACCGATTGGAGTGGCAGAGTCATTCGCGCATTTGTGGAATGCCGGCGGCAGGTTCCCGCAAGGTTTAAGTGGCTTCTTTGCCGGATTCAGGATAGAGGTTTTCTCTTTCGTGAGTATAGAACTGGTAGGTACAACAGCTGCGGCAACCAAAGGTTCAGAGAAATCACTACCACGAGCGATTAACTCCATTCCGATACGTCACATTATg                                                                                                                                                                                                   >  2:238391/1‑251 (MQ=255)
                                                  gATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTAt                                                                                                                                                                                                    <  1:566636/249‑1 (MQ=255)
                                                        cgtcgGCCTTGTCATGGTTGCGATGCACTTTCAGGCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTTGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACg                                                                                                                                                                                            <  1:218758/251‑1 (MQ=255)
                                                        cgtcgGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCGGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACg                                                                                                                                                                                            >  1:101489/1‑251 (MQ=255)
                                                        cgtcgACCTGGTCATGGTGGCGATGCACTTTCAGTCACCGCCTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGGTTAAGTGGATTCTTTGCCGGATTCCAGATAGCTGTTTTCGCTTTCGTGGGGGTTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCAATTCCGCTCCGTATCATTATGTTCtcg                                                                                                                                                                                             >  2:247675/1‑248 (MQ=255)
                                                             gCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTc                                                                                                                                                                                       <  1:343229/251‑1 (MQ=255)
                                                                 ggTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGCTTGAGCTGGTAGGTAAAACAGCTGCGGAAACAAAAGATTCAGAGACATCACTGCCACGCGCGCTTGACTCCATTCCGATCCCTATCATTATGTTCTACGTCTTCacgc                                                                                                                                                                                   >  2:579691/1‑251 (MQ=255)
                                                                 ggTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTcgcg                                                                                                                                                                                    <  1:505793/250‑1 (MQ=255)
                                                                             gATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGGTTAAGTGGCTTCTTTGCCGGATTCCATATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGGGGGAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATAATTATGTTCTACGTCTTCGCGCTGATTGTGAttt                                                                                                                                                                       >  1:3909/1‑250 (MQ=255)
                                                                                     ttCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCg                                                                                                                                                                 <  1:442135/249‑1 (MQ=255)
                                                                                       cAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCTTTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGAc                                                                                                                                                             <  1:171894/251‑1 (MQ=255)
                                                                                       cAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGAc                                                                                                                                                             <  1:247174/251‑1 (MQ=255)
                                                                                              cGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTgg                                                                                                                                                      >  1:8017/1‑251 (MQ=255)
                                                                                               gACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGATTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGAAGCCGTggg                                                                                                                                                     >  2:318170/1‑250 (MQ=255)
                                                                                                    gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGGTGATTGTGATTATGTCCGTGACGCCGTGGAGGTCg                                                                                                                                                >  2:187953/1‑251 (MQ=255)
                                                                                                      gtgGAAGCGTCATGCGCGCATTTGCGGAATGACGGCGGCTGGTTCCCGACAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGGTTTCGCTTTCATGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAACGATTCAGAGAAATCCCTGCCACGCGCGATTAACTCCGTTCCGATCCGTATCATTCTGTTCTACGTCTTCGCGCTGATTGGGTTTATGTCCGTGACGCCGTGGAGTTCGGt                                                                                                                                              >  2:233628/1‑251 (MQ=255)
                                                                                                        ggAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAg                                                                                                                                            <  1:192912/251‑1 (MQ=255)
                                                                                                          aaGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGCTCCGTATCATTATGTTTTACGTCTTCGCGCTGATTGTGTTTATGTCCGTGACGCCGTGGGGTTCGGTAGTc                                                                                                                                          >  2:260934/1‑251 (MQ=255)
                                                                                                           aGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTcc                                                                                                                                         <  1:195733/251‑1 (MQ=255)
                                                                                                                      cgcATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGc                                                                                                                              <  1:290892/251‑1 (MQ=255)
                                                                                                                        cATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCGGATAGCGGTTTTAGCTTTCGTGGGGATTGATATGGTAGGTACAACAGCTGCGGAAACCAAAGCTTCAGAGAAATCACTGCCACGCGCGATTAACGCCATTCCGCTCCGTATCATTATGTTCTACGTCTCCGAGCTGCTTGTGATTATGACCGTGACGCCGTGGAGTTCGGAAGGCCCGGTGAAAAGccc                                                                                                                            >  2:428344/1‑251 (MQ=255)
                                                                                                                         aTTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGCTTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCt                                                                                                                                                                                         >  1:661960/1‑189 (MQ=255)
                                                                                                                          tttGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGt                                                                                                                          >  1:387125/1‑251 (MQ=255)
                                                                                                                               ggAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATGCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCt                                                                                                                                                                                  <  1:632045/190‑1 (MQ=255)
                                                                                                                               ggAATGACGGAGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGACGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTACGATCCGTATCATTATGTTCTACGTCTTCGCGCt                                                                                                                                                                                  >  2:632045/1‑190 (MQ=255)
                                                                                                                                          ggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTttg                                                                                                                       >  1:388547/1‑238 (MQ=255)
                                                                                                                                          ggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTttg                                                                                                                       <  2:388547/238‑1 (MQ=255)
                                                                                                                                             tGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCg                                                                                                                           <  2:162802/231‑1 (MQ=255)
                                                                                                                                             tGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCg                                                                                                                           >  1:162802/1‑231 (MQ=255)
                                                                                                                                                     gAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGGTTTCGCTTTCGTGGGGAGTGAGCTGGTAGGTACAACAGCTGCTGAACCCAAAGATTCAGAGAAATCACTGCCACGCTCGATTAACTCCATGCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACTCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGc                                                                                               <  2:230125/251‑1 (MQ=255)
                                                                                                                                                        aGGTTTAAGTGGCTTCTTTTCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCTTTGATTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGctgc                                                                                            <  1:387347/251‑1 (MQ=255)
                                                                                                                                                        aGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGctgc                                                                                            <  1:56834/251‑1 (MQ=255)
                                                                                                                                                         gttttAAGAGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGTACTGTTCGTGTTGGTAGGGctgcc                                                                                           >  2:191370/3‑251 (MQ=255)
                                                                                                                                                                 ggtcttcgtTGCCGGATTCCAGATAGCGGTTTTCGCTGTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACTCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCg                                                                                   <  1:247675/248‑1 (MQ=255)
                                                                                                                                                                  ggCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGc                                                                                  >  2:204303/1‑251 (MQ=255)
                                                                                                                                                                  ggCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGc                                                                                  >  1:304562/1‑251 (MQ=255)
                                                                                                                                                                      tcttTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGc                                                                              >  1:145343/1‑251 (MQ=255)
                                                                                                                                                                               gattccctagaGCGGTTTTCGCTTTCGTGGGGAGTGAGCGTGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTTCCGCAAGCGTGATCAAc                                                                     <  2:44920/241‑1 (MQ=255)
                                                                                                                                                                                            ggTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTgc                                                        >  1:385735/1‑250 (MQ=255)
                                                                                                                                                                                                 tCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACctct                                                   >  2:141987/1‑251 (MQ=255)
                                                                                                                                                                                                       gcGTGGGGATTGGGCTGGTCGGTACAACAGCCGCGGAAACCAAAGACTCAGAGAAATCACGGACACGCGCGATTACCTCCACGCCGAGCCGTATCATTATGTTCTACGTCCTCGCGCTGGTTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGCCCCGGTGAAAAGCCCGGTTGTTGAACTGTTCGGGGTGGTAGGGCTGCCGGCTGCCGCAAGCGTTATCAACTTTGTGGTGCTGATCTCTgcggc                                              <  2:385735/249‑1 (MQ=255)
                                                                                                                                                                                                                  gAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGAc                                                                                                                                                             <  1:386160/128‑1 (MQ=255)
                                                                                                                                                                                                                  gAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGAGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGAc                                                                                                                                                             >  2:386160/1‑128 (MQ=255)
                                                                                                                                                                                                                        gtaggtaCAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTTAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAgcgg                            <  2:8017/251‑1 (MQ=255)
                                                                                                                                                                                                                          aggtaCAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAgcggcg                          <  1:580694/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGc            >  1:305701/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGc            >  1:141856/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ccAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATg      <  1:204303/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ccAAAGATTAAGAGAAATCCCTGCCTCGCGCGAGTAACTCCATTCCGATCCGTATCATTATTTTCTACGTCTTCGCGCTGATTGGGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTATTCTCTACCAGCCGTATg      <  2:304562/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGCCCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGt  >  1:448218/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGT  >  NC_000913/4430338‑4430830

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