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,4070CT96.4% 157.5 / ‑2.3 58S357S (AGC→AGTthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/25);  total (31/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 = 6.82e-01

CCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAAT  >  NC_000913/1170‑1652
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
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                                                                                                                                                                                                                                                <  1:145823/113‑1 (MQ=255)
                                                                                                                                ccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGttt                                                                                                                                                                                                                                                >  2:145823/6‑118 (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)
                                                                                                                                                                               gCAGCGATGTCACGCGCCCGTATTGCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGcgct                                                           <  1:177271/251‑4 (MQ=255)
                                                                                                                                                                                             cgcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCCGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAata                                             <  1:405549/251‑1 (MQ=255)
                                                                                                                                                                                              gcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAatat                                            >  1:375113/1‑251 (MQ=255)
                                                                                                                                                                                                  gTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCGCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAc                                        <  2:92439/251‑1 (MQ=255)
                                                                                                                                                                                                    aTTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACa                                       <  1:121379/250‑1 (MQ=255)
                                                                                                                                                                                                     tttCCGTGGTGCTGATTACGCAATCATCTTCCGAATCCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGGGGTAGGGGGTGGGATGCGCACCTTGCGAGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCCATATCAACAtt                                     >  2:61588/1‑251 (MQ=255)
                                                                                                                                                                                                                cTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTATCTCGGTGGGAGGGGATGGTAATAGCAACTGTGGTGGGATATCGGCGCAATACTTTTACGCCGTGGCCCGCGCCCATATCCACATTGTTGccctagt                          >  2:21150/1‑246 (MQ=255)
                                                                                                                                                                                                                cTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGGGGCAATATCAACATTTTCGCCATTGc                          >  1:81001/1‑251 (MQ=255)
                                                                                                                                                                                                                  gATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTAGTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGc                          <  1:144553/249‑1 (MQ=255)
                                                                                                                                                                                                                   gtttcGCAATCATCTTTCGTATACTGGATCAGTTTCGGCGTTCCGCATAGCGATGGGGTGCGAGCTAAACGGGCAATGCAGGAAGATTTCGACGTGGAAGTGACAGAAGGCTGACTGGAGGCGTAGGCGGTGACGGAACGGCTGGCCATGATCGCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAGTATCAACATTGTCGCCATTGCTCa                       <  2:239051/247‑1 (MQ=255)
                                                                                                                                                                                                                        gCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAt                  <  1:49764/251‑1 (MQ=255)
                                                                                                                                                                                                                           atcatcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAtctt               >  1:146347/1‑251 (MQ=255)
                                                                                                                                                                                                                            tcatcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTTCTCAGGGGtcttc              >  1:481598/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ttCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAgggg                                                                                                             >  1:428648/1‑150 (MQ=255)
                                                                                                                                                                                                                                 ttCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTg                                                                                                             <  2:428648/152‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGCTGGCCATTATCTCGGTGGTAGGCGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTTCCGCGCTGGCCCGCGCCAATATCAAAATTGTCGCCATTGCTCAGGGATCTTCTGaaagc       >  1:215168/1‑248 (MQ=255)
                                                                                                                                                                                                                                       aTACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGa                                                                                                            <  1:333147/147‑1 (MQ=255)
                                                                                                                                                                                                                                        tACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTTAACGGGCAATGCAGGCAGAGTTCTACCTGGAACTGAAAGAATGCTTACTGGAGCCCCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGCGGCAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCTCCAATATCAACATTGCCGCCATTGCTCAGGGATCTTCTGAACGCTCAAt  >  1:143812/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
CCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAAT  >  NC_000913/1170‑1652

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