Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,341 C→G V335V (GTC→GTG 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,3410CG98.2% 148.5 / ‑3.4 54V335V (GTC→GTGthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (28/25);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

AGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCG  >  NC_000913/1095‑1575
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
aGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGAATGGCGGCGCGCGTGTTTg                                                                                                                                                                                                                                        <  1:466144/251‑1 (MQ=255)
    aTGGATCTGTCCTACGTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTGTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGc                                                                                                                                                                                                                                    <  2:21671/251‑1 (MQ=255)
                 aCTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCAGCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATATTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGcc                                                                                                                                                                                                                       <  1:111699/251‑1 (MQ=255)
                      ggCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTAt                                                                                                                                                                                                                  >  1:432366/1‑251 (MQ=255)
                      ggCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGGGTTTGCAGCGCTGTCACGCGCCCGTAt                                                                                                                                                                                                                  >  2:189771/1‑251 (MQ=255)
                        cgcTACAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTTACGCGCCCGGAttt                                                                                                                                                                                                                >  2:458866/1‑251 (MQ=255)
                         gcTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCAGGCCCGGGGAGGAAAGGGATGGTTGGCATGGCGGCGCGCGTGGTTGGAGCGGTGTCACGCGCCCGGATTTc                                                                                                                                                                                                               >  1:86567/1‑251 (MQ=255)
                                ttcttcACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTAATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGGCACGCGCCCGTATTTCCgtggtg                                                                                                                                                                                                        >  1:49138/1‑251 (MQ=255)
                                    tcACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCt                                                                                                                                                                                                      <  1:63015/249‑1 (MQ=255)
                                      ccacccGCACCATTACCCCCATCGCCCAGGTCCAGATACCGTTCCTGATTAAAAATCCCGGAAATCCTCAAGCTCCAGATACGCGCATTGGTGCAAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATGTCCAATCGGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGTTGAtt                                                                                                                                                                                                  <  2:118772/248‑1 (MQ=255)
                                                 aTTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAatcatc                                                                                                                                                                                       <  1:414361/251‑1 (MQ=255)
                                                       cccATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTAGTTAAGCCATCATCTTccc                                                                                                                                                                                  >  1:78371/1‑249 (MQ=255)
                                                         cATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAAt                                                                                                                                                                               <  1:180192/251‑1 (MQ=255)
                                                               ccAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGAGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTACGCAATCATCTTCCGAATACAgt                                                                                                                                                                          >  1:124887/1‑249 (MQ=255)
                                                                           cccTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGTTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCg                                                                                                                                                             <  1:279987/251‑1 (MQ=255)
                                                                                       tAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGCTGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCCTAGTTCCGTGGTGCTGATTACGCAATCATCTTCCGACTACAGTATCCGTTTCTGCGTTc                                                                                                                                                          >  1:186718/1‑242 (MQ=255)
                                                                                       tAAAAATACCGGAAATCCACAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTc                                                                                                                                                          <  2:186718/242‑1 (MQ=255)
                                                                                        aaaaaTACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCg                                                                                                                                                 >  1:376125/1‑250 (MQ=255)
                                                                                          aaaTACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACt                                                                                                                                              <  1:458866/251‑1 (MQ=255)
                                                                                                                ccAGGTACGCTACTTGGTGCCAGCTGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGGTCAGCGTTTCCTGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCTCGCGCCCGTATTTCCGTGGTGCTGATTACGCAAGCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGATTGTGTGCGAGCTGAACGGGCa                                                                                                                          <  2:35869/249‑1 (MQ=255)
                                                                                                                  aGGTACGCTCCTTGGTGCCAGCCGTGATGATGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATGACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc                                                                                                                      <  2:432366/251‑1 (MQ=255)
                                                                                                                   ggTAAGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGCAAGGGTTGGTTGGCATGGCGGAGCGCGTGTTTTCAGCGGTGGCAGGCGCCCGTATTTCCGTGGAGCTGATTACGACATCACCTCCCGAATACAGATTCAGTTTCTGCCTTCCTCAAAGGGACTGGGTGCGCGCTGCACGGGCCATGcc                                                                                                                     >  2:351390/1‑250 (MQ=255)
                                                                                                                                       ccGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCGGACCCGGGGATGACAGGGATGGTTGGCATGGCGGCGCGCGTGGTTGCAGCGATGTCAGGCGCCCGGATTTCGGTGGTGGTGCATCCGCAATAATCTTCCGAATACAGTTTCAGGTTCTGCGGTCCGCAAAACGCCTGTGTGCGGGCTGTACGGGCAAAGAAGGCAGGGTTCCACCTGGaa                                                                                                  >  1:193211/1‑250 (MQ=255)
                                                                                                                                              gaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAaaga                                                                                          <  1:233208/251‑1 (MQ=255)
                                                                                                                                                    gAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCTGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGGTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCtt                                                                                    <  2:352514/251‑1 (MQ=255)
                                                                                                                                                    gAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCtt                                                                                    >  1:151720/1‑251 (MQ=255)
                                                                                                                                                            ggTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGGGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGGAGGCTTACTGGAGc                                                                            >  2:49764/1‑251 (MQ=255)
                                                                                                                                                                 aGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTg                                                                       <  1:189771/251‑1 (MQ=255)
                                                                                                                                                                     cATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGGGGTGCTGTTTACGCAATCATCTTCCGAATACAGTATCAGTTTCCGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGGGCCGCTGGcgg                                                                   >  2:405549/1‑249 (MQ=255)
                                                                                                                                                                      aTTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTAACGCGCCCGTATTTCCGTGGTGCTGTTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAGCTGCAAGACGGCTTACTGGAGCCGCTGGcggg                                                                  >  2:121379/1‑248 (MQ=255)
                                                                                                                                                                           cAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAACGACTGTGTGCGAGCTGCACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTTGAGCCTCTGGcggagtcgg                                                             >  1:356379/1‑243 (MQ=255)
                                                                                                                                                                               cggAATAAGATGGCAATGGTCAGCGTTTCTGACCCGGGGATGAAAGGGATGGGTGGCATGGCGGCGCGCGTGTATGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAAACAGCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAGTGCAGGAAGAGATCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGAGGTGACGTAACg                                                         <  2:46057/249‑1 (MQ=255)
                                                                                                                                                                                 gAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACggg                                                       >  1:231771/1‑250 (MQ=255)
                                                                                                                                                                                  aaTAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGTCGGCATGGCGGCGCGCGTCTTTGCAGCGCTGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACCGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCCGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTAGTGGAGCCGCTGGCAGGGACGGAACGGCTg                                                     >  2:144553/1‑251 (MQ=255)
                                                                                                                                                                                    taataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCAGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGAGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGAGACGGAACGGATgg                                                    >  1:216803/5‑251 (MQ=255)
                                                                                                                                                                                    taataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGCGGAACGGGTgg                                                    >  1:471764/5‑251 (MQ=255)
                                                                                                                                                                                    aaagaTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTgg                                                    <  2:78371/247‑1 (MQ=255)
                                                                                                                                                                                      ataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGGGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGCAAGAAGGCTTACTGGGGCCGCTGGCGGTgcagggaagggtgggc                                                  >  2:224728/3‑236 (MQ=255)
                                                                                                                                                                                         tGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTGTGCGAGCTGAACGGGCAATGCAGGGAGAGTTCTACCTGGAACTGAAAGAAGGCCTACTGGAGCCTCTGGCGGTggagggacaggcgggcact                                               >  2:321932/1‑233 (MQ=255)
                                                                                                                                                                                         tGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGcgggggcagatcggaagaggccc                                               >  1:228405/1‑229 (MQ=255)
                                                                                                                                                                                         tGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAc                                                               <  2:228405/235‑1 (MQ=255)
                                                                                                                                                                                          ggCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGcccaaaggt                                               <  1:126821/250‑9 (MQ=255)
                                                                                                                                                                                              aTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGGTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTAtctc                                          <  2:231771/251‑1 (MQ=255)
                                                                                                                                                                                                      cGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCCGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCCTCTTCCGAATACAGTATCAGTTTCTGCGGTCCGCAAAGCGACTGTGTGAGAGCTGAACGGGCAATGCAGGCAGAGTTCTACCTGGAACTGGAAGAAGGCTTTCTGGGGCCGATGGCGGTGAGGGGACGGCTGGGCAttttcccgggggtgg                                  >  2:482138/1‑238 (MQ=255)
                                                                                                                                                                                                       gTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGCTGGCCATTATCTCGGTGGGAgg                                 >  1:92439/1‑251 (MQ=255)
                                                                                                                                                                                                       gTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACGGAACGGGTGGCCCTTTTCTCggggggggg                                 >  1:246939/1‑244 (MQ=255)
                                                                                                                                                                                                           ccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGttt                                                                                                                                                                  >  2:145823/6‑118 (MQ=255)
                                                                                                                                                                                                           ccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGttt                                                                                                                                                                  <  1:145823/113‑1 (MQ=255)
                                                                                                                                                                                                               cccGGGGATGAAAGGGATGGTTGGCATGGCGGCGAGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATggg                          >  1:395633/2‑249 (MQ=255)
                                                                                                                                                                                                                ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTAt                        <  1:12345/251‑1 (MQ=255)
                                                                                                                                                                                                                   gggATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTAACAATTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCCGGAAGAGTTCTACCTGGAACTGAAAGAGGGCATACTGGAGCCGCTGGCGGGGACGGAGAAGCTGGCCCTTATCTCgtggaggggggagggtatgtc                     >  2:398521/1‑231 (MQ=255)
                                                                                                                                                                                                                         aaaGGGATGGGTGGCAGGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCtt               <  2:457735/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tggCATGGTGGCGGGCGTGTTGGCCGCGATGTCACGCGCCCGTAGTTCCGTGGTGCTGATTACGCAATCAGCTTCCAAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGATTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGTTTACGGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGAtc     <  2:25293/249‑1 (MQ=255)
                                                                                                                                                                                                                                      gCATGGCGTCGCGCGTGTTTGCAGCGATGTCCCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCg  <  1:357808/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
AGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCG  >  NC_000913/1095‑1575

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