Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,950,171 A→G 20.5% E311G (GAG→GGG)  argA → amino acid N‑acetyltransferase and inactive acetylglutamate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,950,1710AG20.5% 124.8 / 5.5 56E311G (GAG→GGG) argAamino acid N‑acetyltransferase and inactive acetylglutamate kinase
Reads supporting (aligned to +/- strand):  ref base A (3/40);  new base G (9/3);  total (13/43)
Fisher's exact test for biased strand distribution p-value = 6.33e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGG  >  NC_000913/2949938‑2950404
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
ggTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTTTGGGGTTgt                                                                                                                                                                                                                                            >  2:96078/1‑238 (MQ=255)
  aaaaaGCCCCGGAAGGGGAAGGCGATTACAACTCCGGTACGGTGCGCTTTGTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGGTGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGGGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGAt                                                                                                                                                                                                                                           <  2:832942/235‑1 (MQ=255)
     aaGCCCAGGAAGAGAAAGGCGATTACAACCCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGc                                                                                                                                                                                                                                       <  2:364835/239‑1 (MQ=255)
      atcccGGGTAGAGAAAGGCGATTACCACGCCGGTATGGTGCGCTTTTTGCGTGGCCCAGTGAAAGCCTGCCGCAGCGGCGTGCGGCGATGTTATTTAATCAGTTACCAGGAAGATGGCGTGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGcc                                                                                                                                                                                                                                      <  1:316794/237‑1 (MQ=255)
      aGCCCGGGAAGAGAAAGGCGTGTCCACCGCCGGTCCGGGGCGCTGTGTGGGTCGCGAAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTTCGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGATTTGATTCGcc                                                                                                                                                                                                                                      <  1:1156613/239‑1 (MQ=255)
       gCCCAGGAAGAGAAAGGGGTTTACCCCCCCGGTTCGGTGCGTTTTTTGCGTGGCCCAGTGAAAGCCTGGCGCAGCGGCGTGCGTCGCTGTCATTTAAGCAGGTAGCAGGAAGATGGCGCGCTGTGGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGcc                                                                                                                                                                                                                                      <  1:498364/238‑1 (MQ=255)
            ggAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGCAGATGGCGCGCTGTTGCAGAAGTTGTTCTCACCCACCGGGATCGGTACGCCGTTTGGGATGCAAAACGCCAAGACAATCCGCGCGCCAACAACCAAAGAAATTGGCGGTGTTTTGGGGtggtgttagcccccgg                                                                                                                                                                                                                                >  1:1139099/1‑224 (MQ=255)
                 agaAAGGGCGTTGCAACCCCGGTACGGTGGTCTTTTTGCGGTGGGGAGGGAAAGCCTGGCGCAGGGGCGTGGGGCGCTGGCATTTAATCAGTTGTCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGagca                                                                                                                                                                                                                            <  1:1261064/238‑1 (MQ=255)
                   aaaGGCGATTACAACTCCGGTACGGTTCCCTTTCTGGGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAgcagca                                                                                                                                                                                                                         <  2:865964/239‑1 (MQ=255)
                   aaaGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGCTATTGGCGGTATTTTGGGGTTTGTTCGCCCCCTGTagacgcc                                                                                                                                                                                                                         >  2:368809/1‑234 (MQ=255)
                    aaGGGGATTCCAACTCCGGTACGGTGCGGTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAgcagca                                                                                                                                                                                                                         <  2:1136547/238‑1 (MQ=255)
                    aaGGCGATTCAACCTCCGTTAGGGGGCGTTTTTTGGGTGGCGCGGGGAAAGCCTGCCGCACGGGGGTGCGTCCCTGTCATTTAAGCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTAGCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAgcagca                                                                                                                                                                                                                         <  1:546421/238‑1 (MQ=255)
                    aaGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAgcagca                                                                                                                                                                                                                         <  2:878838/238‑1 (MQ=255)
                        cggttgcagttCCGGTACGTTGCGTTTTTTGCGGGGCGCGGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGGCATTTAATCAGTTATCAGGAATATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAgg                                                                                                                                                                                                                      <  1:96078/227‑1 (MQ=255)
                          aTTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCTGGGGTTTGGTTGCCCACCTGGGCAGCAAggggtt                                                                                                                                                                                                                  >  2:643421/1‑235 (MQ=255)
                                            tGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGGGTTATTTGCCCCACGGCGGAGCAAGGTATTCTGGTATGCCCtttt                                                                                                                                                                                                   >  2:1042750/1‑234 (MQ=255)
                                              cgcTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCTTGAGGTGATTCTCCCCCCGGGGCAGCAAGGGAGTCtttgaagcccgttccccc                                                                                                                                                                                               >  2:270929/1‑221 (MQ=255)
                                                    ttGCGTGGCGCAGTGAAAGCCTGCCGCAGGGGCGTGCGTCCGTGTCATTTAAGCAGTTATCAGGAAGGTGGCGCGCGGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCt                                                                                                                                                                                        <  2:692531/239‑1 (MQ=255)
                                                     tGCGGGGCGCAGTGAAAGCCGGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTg                                                                                                                                                                                       <  2:584268/239‑1 (MQ=255)
                                                       cGTGGCGTAGTGAAACCCTG‑‑GCAGCGCGGTGCGTCGCTGTCATTTAAGTGGTTAGCTGTAATATGGCGCGCTGTTGCAAGATTTGCTCCCACGCGACGGTATCGGTACACAGATTGTGATGGAAAGCGCCGATCAGATTCGTCCCGCAACAATCAATGATATTGGCGGTATTCTCGAGTTGATTCGCGCTCTGGAGCAGCAAGGTACTCTGGTACGCCGTTCTCGCGAGCAGCTGgaga                                                                                                                                                                                   <  1:19445/239‑1 (MQ=255)
                                                                         tgcgggtgtggcgtgcgggtgtGTCATTTAATCTATGTGCAGGAAGGTGGTGCGCTGTTGGAAGAGTTGTTCTCCCGGGACGGTATCGGTACGCAGATTGTGATGGAAATCGGCGGGCAGTTTCGGCGGCGAACAATCATCGGTATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGGTCTCGCGAGCAGCTGGAGATGGAAATCGACAAAtt                                                                                                                                                                   <  1:970020/218‑1 (MQ=255)
                                                                         tGCCGCAGGGGCGTGCGTCGCTGGCATTTAATCAGTTAACAGGAAGATGGCGCGCTGTTGCAAGAGTTGGTCTCACGCGACGGTAGCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAAtt                                                                                                                                                                   <  2:276355/239‑1 (MQ=255)
                                                                         tGCCGCAGCGGGCTGCGTCGCTGTCATTTAATCAGTTACCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAAtt                                                                                                                                                                   <  2:1079236/239‑1 (MQ=255)
                                                                           ccGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGCATTCTTGGGTTGATTCTGCCACCGGAGCAGCAAGGTATTCtt                                                                                                                                                                                                               >  2:42418/1‑192 (MQ=255)
                                                                             gCAGTGGCATGCGTCGCTGGCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGGGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCAcc                                                                                                                                                               <  2:915692/239‑1 (MQ=255)
                                                                                          tCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAgcgc                                                                                                                                                  <  2:388219/238‑1 (MQ=255)
                                                                                                            ttATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTTGGGTTGATTTGCCCACCTGCGCCGCCAGGTATTCTGGTACGCCGTTCTCCCGCGCAGCTGGGGATTGGAATCGACAAATTTACCCTTATTTAGCGCGATAAAACGACTATTTcc                                                                                                                               >  2:1018675/1‑238 (MQ=255)
                                                                                                                 aGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGcgcc                                                                                                                        <  2:493471/239‑1 (MQ=255)
                                                                                                                   gAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAaca                                                                                                                                           <  2:895469/220‑1 (MQ=255)
                                                                                                                    agtgtGGCGCTCTGTTGTAAGCGTTTTTCTCACGGGACGGGCTCGGTAGGCTGATTGTGATGGAAAGCGCCGTTCAGTTGCGTCGAGCAATAATCAACGGTAGTGGCCGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACTCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCTCCACTCTTCAGCGCGATAACACGACTATTGCCTGcgcc                                                                                                                        <  1:623348/232‑1 (MQ=255)
                                                                                                                      gATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCTTGCGTTGATTCGCCCACTCGCGCAGCAACGGATTCTGCGACGCCGCTTTCGTGGCGCGCTGGGGGTGGGAAAACACAAAATCACCATTATTCTCCGCGGTAAcaccacaaatatccgcggcgcgg                                                                                                                      >  1:296803/1‑221 (MQ=255)
                                                                                                                           cgcgcTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGGTTGATTTGCCCCCCGGAGCAGCAAGGGATTCTTGTACGCCGTTCTCGCGGGCAGCCGGAGATGGAAATTGACAAATTCACCCTTTTTCAGCGCGATAAAACCGCCATTGACTggcacaggccccta                                                                                                               >  2:368667/1‑225 (MQ=255)
                                                                                                                             cgcTGTTGCAAGGTTTGTTTTCACGGGCCGGTATCGGTTAGCGAGTTGGGAAGGAAAGGCCCGGGCAGGTTCGTCGCGGAACAAGCAACGATGATGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGGAGCAATGTATGCTGGTTCGCCGTTCGCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGACGCGCTCTATcc                                                                                                             <  1:550897/238‑1 (MQ=255)
                                                                                                                                   tGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCg                                                                                                    <  2:1105660/239‑1 (MQ=255)
                                                                                                                                          gttgttCTCACGCGACGGTATCGGTACGCAGATTGGGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGgaagag                                                                                              <  2:596445/238‑1 (MQ=255)
                                                                                                                                          gttggtCCCCCGCGCCGGTATCGGTACGCAGATTGTGAAGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGgaagag                                                                                              <  2:762097/238‑1 (MQ=255)
                                                                                                                                           ttgttCTCACGCGGCGGTATCGGTACGCAGATTGGGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGgaagaga                                                                                             <  2:564500/238‑1 (MQ=255)
                                                                                                                                              ttCTCACACGCCGGTTTCGGTACGCAGATTGTGTTGGAACTCGCCGATCAGTTTCGTCGTGCAACAAACAACGGTATGTGCGGTAGTCTGGAGTGGAGTAGCTGACTGGAGCAGAAAGGTATTTTGGTTCGCCGTTCTCGCGAGCAGCTTGAGATTGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGGTCCCGGAAGAGAAGAt                                                                                        <  1:592027/239‑1 (MQ=255)
                                                                                                                                                tctccgggtcggtttgtgttccgaGATTGTGAGGGAAAGGCCCGGGCAAGTTCGTGCCGCAAAAATCCACGATATTGGCGGTATTCTGGAGTTTATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTg                                                                                      <  1:643421/216‑1 (MQ=255)
                                                                                                                                                ctcACGCGCCGGTTTCGGTACGTAGATTTGGATGGAAAGCGCCGAGCAGGTTCGTCGCGGAACAATCAACGATATTGGCCGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTg                                                                                      <  2:554754/239‑1 (MQ=255)
                                                                                                                                                                 gTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGcgccg                                                                                                                       <  2:1000315/192‑1 (MQ=255)
                                                                                                                                                                    cGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTGTTGAGGTTGTTTGGCCACCGGGGCAGCAAGGGATTCTGTGGTACCGTTCTCGCCGGCAGCCGGAGACGGATATTGCACAATTCACCATCATTTTAGCCGCAAACAGGACCATTTCCTTGCGCGCCGCCTATCCGCTTCCCGGAGAGAAGACTGGGGGAAATTCCCTGTggggg                                                                   >  1:982817/1‑238 (MQ=255)
                                                                                                                                                                         attggtaTGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTg                                                                                      <  2:242047/208‑1 (MQ=255)
                                                                                                                                                                                   aaaGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGTTTTTGGGGGTGATTCGCCCACCGGGGCAGCAAGGGATTTTGGTACACCGTTCTCGCGAGCAGCCGTGGATGGGAAATGACCAATTTACCACTATTTATCGCGGTAAAAAGGACATTTCCTGCGCCCGCCTCCATCCTGTGCCCGAAGAAAAACATTTGGGAAAAGCCCTGGGGGGGGGTTCCCCCgagaa                                                     >  1:452476/1‑234 (MQ=255)
                                                                                                                                                                                     agggcCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCg                                                 <  2:1053250/236‑1 (MQ=255)
                                                                                                                                                                                      gcgcCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGATTTGATTCGCACACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATTGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCg                                                 <  2:1134883/238‑1 (MQ=255)
                                                                                                                                                                                          cgcgCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGt                                             <  2:1006467/235‑1 (MQ=255)
                                                                                                                                                                                                  ttCGTCGCGAAACAATCAACGATATTGGCGGTATTCTGGGGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAg                                     <  2:550765/238‑1 (MQ=255)
                                                                                                                                                                                                      tcgCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAac                                                                                                                                            <  2:569488/136‑1 (MQ=255)
                                                                                                                                                                                                       cgcgCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGa                               <  2:296803/239‑1 (MQ=255)
                                                                                                                                                                                                                 aTCAACGATATTGGCGGTATTCTGGAGTTGAGTCGGCCACCGGAGCAGCAAGGGATTCCGGGACCCCCTTCTCGCCGGCCGCTGGGGGGGGGAATCAGAAAATTTCCCACTATTCAGCGCGGGCACACGACATTTTCCCTGCCCGCCCCCCACTCGATCCCGGCAGGCAGAGTTGGTGGAAAGGCCCGGCTGGCCGTTTACCCGCATTTACGCAGTTCCTCAAgcgggaagagtcg                         >  1:817391/1‑227 (MQ=255)
                                                                                                                                                                                                                  tCAACGATTTTGGCGGTAGTCTGGAGTTGAGTCGCCCACTGGAACAGCAAGGTAGTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGTGTGAAGTTctgctg                     <  2:958346/239‑1 (MQ=255)
                                                                                                                                                                                                                         tatTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                                                                                                                                                                                    <  1:562753/47‑1 (MQ=255)
                                                                                                                                                                                                                                gggaTTCTGGGGTTGATTCGCCCACTGGAGCAGCAAGGTATTGTGGTACGCCGTTCTCGCGAGCAACTGGAGATGGAAATCGACAAATTCCCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGc       <  2:791013/235‑1 (MQ=255)
                                                                                                                                                                                                                                     tCTGGGGTTGATTCCCCCAGTGGAGCAGCAAGGTATTCCGTTACGCCGTTCTCGCAAGCTGCTGGTGATGGAAAACCACAAATTCACCATTATGCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAgg  <  2:1253614/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAgg  <  2:753477/239‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
GGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGG  >  NC_000913/2949938‑2950404

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