Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,407 C→T S357S (AGC→AGT)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,4070CT97.1% 209.1 / NA 71S357S (AGC→AGT) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (37/31);  minor base G (1/1);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATC  >  NC_000913/1169‑1653
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
tCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATTCAGCATCAGTTTCTGc                                                                                                                                                                                                                                                   >  1:422321/1‑251 (MQ=255)
   cTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGGGTGTTTGCAGGGATGTCACGGGCCCGTATTTCCGGGGGGGTGATTTAGAGATAATCTTCCGAATACAGTATCAGTTTCTGCggt                                                                                                                                                                                                                                                >  2:268555/1‑249 (MQ=255)
      gCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGGACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCACGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTTAAGCGATGTCACGCGCCCGTATTTCAGTGGTGCTGATTACACAATCATCTTCCTAAACCAGTATCCGTTTCTGCTTTccg                                                                                                                                                                                                                                             >  1:117153/1‑250 (MQ=255)
             tAAAAATACCGGAAATCCTCAAGCTCCAGGTACGATCATTGGCGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGAGCGCGTGTTTGCAGCGAGGTCACGAGCCCGTATTTCCGTGGAGCTGATTACGCCATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCg                                                                                                                                                                                                                                      >  1:523969/1‑251 (MQ=255)
              aaaaaTACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCg                                                                                                                                                                                                                                      >  1:543366/1‑250 (MQ=255)
                 aaTACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCCtg                                                                                                                                                                                                                                  >  1:241550/1‑251 (MQ=255)
                             ccTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGAGCGCGTGCTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAAACCTCTTCCGAATACAGTATCAGTTTCTACGGTTCGCAAAGCGACTGTGTGCgggccg                                                                                                                                                                                                                       >  2:501096/1‑245 (MQ=255)
                                       cAGGTACGCTCATGGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTGCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATg                                                                                                                                                                                                            <  2:284096/251‑1 (MQ=255)
                                        aGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc                                                                                                                                                                                                           >  1:321352/1‑251 (MQ=255)
                                                     ggtggctgcctgtattaaGACGAACTACCGGGCAAGGGCATTTCCAAGCTGAATAATATGCAAATGTTCAGCGTGTCCTGTCCGGGGATGAAAGGGATGGTTGGCTGGGCGGCGCGCTGGGTTACATCAAGGTCACGCTCCCGTATTACCGTGGTGCTGATTACGCAATCATCTTCCGACTACAGGATCAGTTTCTGCTTTCCGCAAAGCGACTGGCTGCGAGCTGAACGGGCAATGCAGCAagtggtcta                                                                                                                                                                                              <  2:250280/235‑8 (MQ=255)
                                                             ccgggatgaatgcgcatgacctgTCAAGGGCTTTGCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCTTGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCTCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGGGCGAGCTGAACGGGCAAGGCAGGAAGAGTTCTTCCTGGAAc                                                                                                                                                                                      <  2:551203/229‑1 (MQ=255)
                                                                tgatgaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGAGTGTTTTCAGCGGTGTCCCGCGCCCGGATTTCCGTGGGGCTGGTTACGCAATCATCTTCCGGATACAGGATATGTTTTTGCGTTCCGCAAAGCGACTGTGTGGGCGCTAGACGGGGAATGCCGGAagggatctcaccgggactg                                                                                                                                                                                    >  2:13668/1‑233 (MQ=255)
                                                                            aTTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCTGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTAc                                                                                                                                                                       <  1:245272/251‑1 (MQ=255)
                                                                               acctgTCAAGGGCATTTCCAAGCAGAATAATATGGCAATGTGCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCGTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTg                                                                                                                                                                     <  2:239904/246‑1 (MQ=255)
                                                                               aCCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGACAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGGGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGACGAGTTCTACCTGGAACTGCaagacagcttcatgg                                                                                                                                                                    >  2:35666/1‑240 (MQ=255)
                                                                               aCCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCCATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTgg                                                                                                                                                                    >  1:324984/1‑251 (MQ=255)
                                                                                      aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTGCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCt                                                                                                                                                             <  1:368464/251‑1 (MQ=255)
                                                                                        gggCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACggg                                                                                                                                                                                                                 >  1:644821/1‑197 (MQ=255)
                                                                                        gggCAGTTCCAAGCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACggg                                                                                                                                                                                                                 <  2:644821/197‑1 (MQ=255)
                                                                                              ttCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTgg                                                                                                                                                     >  1:552754/1‑250 (MQ=255)
                                                                                               tCCAATCTGAATAATATGGCAAGGTTCATCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTTCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGa                                                                                                                                                     <  1:394236/250‑1 (MQ=255)
                                                                                                ccAATCTTAATAATATGTCAATGTGCAGCGTTTCCGGCCCGGGGATTAAAGGGATGGTTTGCATGGGGGCGCGCGTGATTTCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCACTCATCATCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGAAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTTGCGGTGAc                                                                                                                                                    <  2:594500/250‑1 (MQ=255)
                                                                                                  aaTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAC‑ggg                                                                                                                                                >  1:340465/1‑250 (MQ=255)
                                                                                                  aaTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGGTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCa                                                                                                                                                                                                               >  1:246934/1‑189 (MQ=255)
                                                                                                   aTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGACACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCCATGCAGGAAGAGGTTTAACTGGAACTGAAAGAAGGCTTACAGGGGCCGCTGGcgggggggg                                                                                                                                                  >  2:87618/1‑241 (MQ=255)
                                                                                                   aTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTgg                                                                                                                                                                                         <  1:643848/210‑1 (MQ=255)
                                                                                                   aTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAACGACTGTGTGCGAGCTGAACGGGCACTGCAGGAAAAGTTCTACCTgg                                                                                                                                                                                         >  2:643848/1‑210 (MQ=255)
                                                                                                         ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTg                                                                                                                                         <  1:35666/251‑1 (MQ=255)
                                                                                                         ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTg                                                                                                                                         <  2:266950/251‑1 (MQ=255)
                                                                                                            ataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGcc                                                                                                                                      >  1:251315/3‑251 (MQ=255)
                                                                                                             taTGGCAATGTTCAGCGTTTCCGGCCCGGGGCTGCACGGGCTGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCCCGCGCCCGTATTTCCGTGGTGCTGGTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGGACTTAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACTGGAACgcttgccc                                                                                                                                      >  1:254445/2‑244 (MQ=255)
                                                                                                                ggCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCAtt                                                                                                                                   <  1:604702/250‑1 (MQ=255)
                                                                                                                  cAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGC‑GGAAAGGcggggcactttt                                                                                                                                >  2:237356/1‑240 (MQ=255)
                                                                                                                     tGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGTTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGAAGGAAGAGTTCTACCAGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGGC‑GGGACGGCTgggcacttatt                                                                                                                               >  2:570280/1‑240 (MQ=255)
                                                                                                                         cAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAGATGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGGACGGCAGGCCATTATCTCgaggg                                                                                                                         >  1:24906/1‑247 (MQ=255)
                                                                                                                          agcgtgtccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCggtggt                                                                                                                        <  2:302104/239‑1 (MQ=255)
                                                                                                                                 ccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGAt                                                                                                               <  1:123992/246‑1 (MQ=255)
                                                                                                                                   ggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGGCATTATTTCGGTGGG‑AGG‑GGATg‑g                                                                                                            >  1:374563/4‑251 (MQ=255)
                                                                                                                                     cccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGAGAC‑GGAACGGCTGGCCATTCACTCGGTGGT‑AGG‑GGATg‑gg                                                                                                           >  1:198460/2‑249 (MQ=255)
                                                                                                                                      ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTAt                                                                                                         <  1:3999/251‑1 (MQ=255)
                                                                                                                                      ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTAt                                                                                                         <  1:323048/251‑1 (MQ=255)
                                                                                                                                      ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAC‑TGGACGGCTGGCCCTTATCTCGGTGGT‑AGG‑GGATG‑GTAt                                                                                                         >  1:319608/1‑251 (MQ=255)
                                                                                                                                              gAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTTTCTCGGGGGT‑AGG‑TGATG‑GTATGCGCACCt                                                                                                 >  1:606164/1‑251 (MQ=255)
                                                                                                                                                       gggtggCAGGGCGGCGCGCGTGGTTGCTGAGATATGACGCGTCCGTATTGCCGTGGTGCTGATTGCGCAATCATCTTCCGAATACAGTATCAGGTTCCGCGTTCCGCAAAGCGACTGGGTGCTAGCTGAACGGGCCATGCAGGAGGAGTTCTACTGGGAACTGAAAGAAGGCTTACTGGAGCCGCGGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGa                                                                                        <  2:162740/247‑1 (MQ=255)
                                                                                                                                                           ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGAtct                                                                                     <  1:570280/250‑1 (MQ=255)
                                                                                                                                                           ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCACGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGAtctc                                                                                    <  1:343346/251‑1 (MQ=255)
                                                                                                                                                                  cggcgCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑aaa                                                                            <  1:196652/251‑1 (MQ=255)
                                                                                                                                                                        gcgTGTTTGCAGCGATGTCACGCGCTCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGTATCAGTTTCTGCGTTCCGCAAAGCGACCGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCttt                                                                      <  1:643932/251‑1 (MQ=255)
                                                                                                                                                                            gtttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGACTACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGc                                                                                                                                                          >  2:232861/2‑168 (MQ=255)
                                                                                                                                                                            gtttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGc                                                                                                                                                          <  1:232861/167‑1 (MQ=255)
                                                                                                                                                                            gtttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGGACGGATGGCCATTAACACGG‑GGTGAGGGTG‑TT‑GTATGCGCACCATGGGTGGGATTTCGGCGCAAATT‑TTTGccc                                                                  >  2:268561/2‑250 (MQ=255)
                                                                                                                                                                             tttGCAGCGATGTCACGCGCCCGTATGTCCTTGGTGGTGATTACGCAAGCGTCTACCGAATGCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCGGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCg                                                                  <  2:164807/250‑1 (MQ=255)
                                                                                                                                                                                  aGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCt                                                                                                                                                                                                                        >  1:257070/1‑100 (MQ=255)
                                                                                                                                                                                  aGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTGTTCGAGCt                                                                                                                                                                                                                        <  2:257070/100‑1 (MQ=255)
                                                                                                                                                                                           aCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATG‑GTATGCGGAACTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGccc                                                   >  1:10537/1‑250 (MQ=255)
                                                                                                                                                                                                   gTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAc                                           <  1:150850/251‑1 (MQ=255)
                                                                                                                                                                                                       ttCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GGATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTg                                       >  1:93488/1‑251 (MQ=255)
                                                                                                                                                                                                         ccGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATGTGT‑TGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCACCATTGTc                                     >  2:110947/1‑251 (MQ=255)
                                                                                                                                                                                                           gtggtgCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAaaga                                                                                                                                                                               <  1:461598/116‑1 (MQ=255)
                                                                                                                                                                                                           gtggtgCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAaaga                                                                                                                                                                               >  2:461598/1‑116 (MQ=255)
                                                                                                                                                                                                           gtggtgCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCAACTTGCGTGGGGTCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTTTC‑Gc                                  >  1:655808/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAgg‑g                     <  1:8296/251‑1 (MQ=255)
                                                                                                                                                                                                                        cGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAgg‑g                     >  1:288402/1‑250 (MQ=255)
                                                                                                                                                                                                                        cGCAATCATCTGCCGAATACGGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAGG‑Ga                    <  2:340465/251‑1 (MQ=255)
                                                                                                                                                                                                                          cAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGG‑AGG‑GGATG‑GTCAGGGCACCCTTCGGGGGATCTCGGCG‑AAAATCTTTGCCGCCCTGGCCCCCCCCAAGAACCACATTGTCCGCC‑TTGCTCAgg‑ggt                   >  2:290569/1‑248 (MQ=255)
                                                                                                                                                                                                                               atcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATG‑GGATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTC‑GCCCTTGCTCAGG‑GATCTTCTg             >  1:400836/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ctCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCCACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAAc          >  1:193722/2‑251 (MQ=255)
                                                                                                                                                                                                                                   tCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGa            <  1:352537/248‑1 (MQ=255)
                                                                                                                                                                                                                                    ccGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACg         <  1:274611/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aaTACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGATCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGG‑AGG‑GGAGG‑GGAAGAGAACCCTGAGTGGGCTATCGGCG‑AAAATCTTTGCCGCGCTGGCCCGCGACAAAATAAACAATGTC‑GCACTTGCTCAGGAGA‑CTTCTGAAcccacc     >  2:439439/1‑246 (MQ=255)
                                                                                                                                                                                                                                          aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAAtc  <  1:237356/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
TCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATC  >  NC_000913/1169‑1653

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