Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,358,472 G→A G1250S (GGC→AGC)  gltB → glutamate synthase subunit GltB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,358,4720GA98.1% 182.0 / NA 53G1250S (GGC→AGC) gltBglutamate synthase subunit GltB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/26);  minor base C (0/1);  total (26/27)
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

CCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTA  >  NC_000913/3358225‑3358716
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
cctgtttaaaGAGCTGGCCGGTTTCTACGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGATACTGCCGAACCGCATCAAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCTATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCt                                                                                                                                                                                                                                                    <  1:260455/244‑1 (MQ=255)
  tgctgAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCtat                                                                                                                                                                                                                                                  >  1:248909/1‑251 (MQ=255)
      gAAAGAGCTGGACGGTTTCCCCGCCAACCAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATc                                                                                                                                                                                                                                               <  1:164276/250‑1 (MQ=255)
                 gcGGTTTCCCCGCCAAACAGAAGAAGCTGGCGCTGTCGAAGCTGCTGGAGGCTGCAGACCCGCATCCAGGTAAGGCACGCTACGGCCCCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCAc                                                                                                                                                                                                                                   <  2:206539/250‑1 (MQ=255)
                    gTTTCACCGCCAACAAGCAGAAGCTGGCGCTGTCGAAGCTGCTCGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTCCTGCACCGAAAACAAACCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGc                                                                                                                                                                                                                                >  2:110307/1‑251 (MQ=255)
                     tttCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCg                                                                                                                                                                                                                               >  2:151041/1‑251 (MQ=255)
                        cACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCCGTTGCTGCAACCGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCNATc                                                                                                                                                                                                                            >  2:255829/1‑251 (MQ=255)
                          ccGCCAAACAGCAGAAGCTGGCGCTTTCGAAGCTGCTGGGGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAg                                                                                                                                                                                                                          <  2:266498/251‑1 (MQ=255)
                              cAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGACTGCTGAAAGCACAGTTGCTGCAACAGGCGAAACCGTTTGGCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATACGCAACACCGCCCGTTCTGGCGGCGCGTCGCTTTCAAGCTATAACGCCCAGCCGCACGGCGATCAAGGGc                                                                                                                                                                                                                      >  1:162350/1‑251 (MQ=255)
                                      agAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGcc                                                                                                                                                                                                              <  1:57063/251‑1 (MQ=255)
                                      agAAGCTGGCGCTGTCGAAGCTGCGGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGcc                                                                                                                                                                                                              <  1:233277/251‑1 (MQ=255)
                                      aaaaGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGcc                                                                                                                                                                                                              <  1:110307/249‑1 (MQ=255)
                                            tGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCt                                                                                                                                                                                                        >  1:156560/1‑251 (MQ=255)
                                                  tGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCaa                                                                                                                                                                                                   >  1:251300/1‑250 (MQ=255)
                                                     cGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCg                                                                                                                                                                                               <  1:51725/251‑1 (MQ=255)
                                                      gAAGCTGCTGGAGACTGCCGACCCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGt                                                                                                                                                                                              <  2:248909/251‑1 (MQ=255)
                                                       aaGCTGCTGGAGCCTGCCGAACCGCATCCAGGTAAGTCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGTGCTGGCAGCCGTTCCTATCAAAGCGTa                                                                                                                                                                                             <  1:251051/251‑1 (MQ=255)
                                                              tgGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGCCGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGCACTTCAAc                                                                                                                                                                                      >  2:85468/1‑251 (MQ=255)
                                                                       ccGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCa                                                                                                                                                                             <  1:7793/251‑1 (MQ=255)
                                                                               cacccAGGAAAGGCACCCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCagag                                                                                                                                                                     <  2:114365/248‑1 (MQ=255)
                                                                                 tCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATAAGGGGCTGGCAGCCGATCCTATCAAAGAGTACTTCAACGGCACCGCAGGCCAGAGCt                                                                                                                                                                   >  2:177028/1‑251 (MQ=255)
                                                                                      gTAAGGCACTCTACTGCACCGAAAACAACTCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCTGACGCACGGCGATCAGGGGCTGGCAGACGATCCTATCAAAGCGTACTTCAACGGCACCCCAGGCCACAGCtgccgc                                                                                                                                                              >  1:188623/1‑246 (MQ=255)
                                                                                                               caacCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAAc                                                                                                                                     >  1:2088/1‑251 (MQ=255)
                                                                                                                  accgccgTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTCTCAAAGCGTACTTCAACGGCACAGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGCGGAACTGt                                                                                                                                  >  2:214763/2‑251 (MQ=255)
                                                                                                                      acgTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGATTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCCCAGGCCAGAGCTTCTGCGTGTGGTACGCGGGCGGCTGGGACCTGTACCt                                                                                                                              >  2:196395/2‑251 (MQ=255)
                                                                                                                           tGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCGTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGAcc                                                                                                                          <  1:100155/250‑1 (MQ=255)
                                                                                                                                aCGGCCTGGTGAACGTGCAGTGCCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGAt                                                                                                                    <  2:37385/251‑1 (MQ=255)
                                                                                                                                      tgctgAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAAc                                                                                                              >  2:224781/1‑251 (MQ=255)
                                                                                                                                                  gggggatgccatagccgtatccGTGTGTCGACGAGCGCGAGAGCACACGCTTCTGGTTCGATAACAGAAACACCGCCCGGTCTGTCGGAGCGTCGATTTCACGCAATATCGCCCAAACGAACGGCGCGAAGGGGCTGGCAGCAGAGCCTATCAAAGCGTACTTCAACGGCACCGCCGGCCAGAGCTGCGGCGTGTCGAACGCGGGGTGCGTTGAACTGTGCCTGACCGGTGATGCCAACGCCTATGTctgt                                                                                                  <  2:152338/231‑4 (MQ=255)
                                                                                                                                                      gctgcacaaggtgAAACCGTTTATCGATGAGACCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACCCCCACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGATGGCAGCCGCTCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAg                                                                                              <  1:140570/239‑1 (MQ=255)
                                                                                                                                                        tgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGc                                                                                            >  1:87145/1‑251 (MQ=255)
                                                                                                                                                          cAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCAt                                                                                          >  2:117192/1‑251 (MQ=255)
                                                                                                                                                          cAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGAGATGCCAACGACTATGTCGGTAAAGGCAt                                                                                          >  1:120816/1‑251 (MQ=255)
                                                                                                                                                              aGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATggcg                                                                                      <  2:156560/251‑1 (MQ=255)
                                                                                                                                                               ggCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGGCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATggcgg                                                                                     <  1:59002/251‑1 (MQ=255)
                                                                                                                                                                      ccGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGggcggcgt                                                                              >  1:59811/1‑249 (MQ=255)
                                                                                                                                                                       cGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGAGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGgcgggcgggtt                                                                              >  2:238705/1‑247 (MQ=255)
                                                                                                                                                                                            gCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGggcggccgtatcgaccatcatcctccgggt                                                        >  2:283233/1‑227 (MQ=255)
                                                                                                                                                                                                  ccTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTcc                                                  <  1:77412/251‑1 (MQ=255)
                                                                                                                                                                                                                atTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCa                                    <  1:60986/251‑1 (MQ=255)
                                                                                                                                                                                                                 ttCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCa                                    <  2:144314/250‑1 (MQ=255)
                                                                                                                                                                                                                    gcaccacCGACCGGTCTGTCGGCGCGTCACATTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGaa                                <  2:204779/247‑1 (MQ=255)
                                                                                                                                                                                                                    gCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCTTCCTCCCGGTGGGTCCGCCTTCCGCCGCCATGa                                 >  1:198316/1‑250 (MQ=255)
                                                                                                                                                                                                                     cAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGGTTAATAGCCATTCGTCCTCCGGTTTGGTCCGCCCTCCGCAGCCATGaa                                >  1:173847/1‑250 (MQ=255)
                                                                                                                                                                                                                        cacCGACCGGTCTGTCGGCTCGTCGCTTGCAAGCTATATCGCCCAGACGCACGGCCAGCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTTGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCGGCCGTCTGCAGCCATGatgcaa                           <  2:31923/251‑1 (MQ=255)
                                                                                                                                                                                                                                gTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGGTTAATCGCCATTCATCCTCCGGTTGGTTTCGCCTTCCGCAGCCCTGAAGAAAACATTATc                   >  1:149184/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTTCTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATTAAGCAAGCATTATc                   <  1:251704/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cTGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGGGGCTTAATCGCCCCTCCTCCTCCGGGTGGGTCCGCCCTCCTCAGCCATTAAGCCAAGCTTTAtgtg                 >  2:210263/1‑248 (MQ=255)
                                                                                                                                                                                                                                    tGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTaa                                                                                                >  2:128085/1‑171 (MQ=255)
                                                                                                                                                                                                                                    tGTCGGCGCGTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCAGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTaa                                                                                                <  1:128085/171‑1 (MQ=255)
                                                                                                                                                                                                                                            cgTCGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTCAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGaagcaa                           <  1:183544/231‑1 (MQ=255)
                                                                                                                                                                                                                                               cGCTTTCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGACAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGGTGGGTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACAcctccctg    >  1:132477/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ttCAAGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTa  <  1:151041/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCGCCCAGACGCACGGCGATCAGGGGCTGGCAGCCGATCCTATCAAAGCGTACTTCAACGGCACCGCAGGCCAGAGCTTCGGCGTGTGGAACGCGGGCGGCGTGGAACTGTACCTGACCGGTGATGCCAACGACTATGTCGGTAAAGGCATGGCGGGCGGCTTAATCGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTA  >  NC_000913/3358225‑3358716

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: