Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,509 A→G A391A (GCA→GCG 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,5090AG93.3% 166.3 / 0.3 60A391A (GCA→GCGthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (32/24);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTG  >  NC_000913/1270‑1750
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
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)
          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)
                                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)
                                                                                                                                         taTCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCATTATCTCGGTGGTAGGTGCTGGTATGCGAACCTTGCGTGGGGTCTCGGCGAAATTCTTTGCCGCGGTTGCCCGCGCCCATATCAACATTGTCGCCATTTCTCAGGGGTCTTTTGCACGCGCCATCTCTg                                                                                               >  2:395493/2‑251 (MQ=255)
                                                                                                                                                 tCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCCCAATCTCTGTCGGGGt                                                                                        >  2:397278/1‑250 (MQ=255)
                                                                                                                                                      gTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATaa                                                                                  <  1:184590/251‑1 (MQ=255)
                                                                                                                                                            caAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGAACGGATGGCCATTATCTCGGGGGTAGGGGATGGTAAGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGGTGGCCCGCGCCAATATCAACATTGCCGCCATTGCACCGGGATCTTCTGAACGCTCATTCTCTGTCGTGGTAAATAACcatg                                                                             >  1:363077/1‑250 (MQ=255)
                                                                                                                                                             aaaGCGACTGTGTGCGAGCTGAACGGGAAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGGGGGAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCAGCGATGGCCCGCGCCAAAATCAAACTTGTCGACATTGCTCAGGGATCTTCTTAACGCCCAATATCTGTCGAGGAAAATAAAgatgat                                                                           >  1:156817/1‑251 (MQ=255)
                                                                                                                                                               aGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAACGCAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGCTGGCCATTATCTCGGTGGTAGGCGAAGGTATGCGAAACTTGCGTGGGATCTCGGCGAAGTTCTTTGCCGCACTGGCCCGCGCCCATATCAAAATTGTCGCCCTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGGAAATAACGATGATGc                                                                         >  2:460917/1‑251 (MQ=255)
                                                                                                                                                                   aCTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGAcc                                                                     <  1:61588/251‑1 (MQ=255)
                                                                                                                                                                           cgcgCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTATCTCGGTGGTCGGTGGTGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTTTTCACGCTCAATCCCTGTCGTGGTAAATAACGGTGATGCGACCACCGGCgg                                                             >  2:289490/4‑250 (MQ=255)
                                                                                                                                                                            gAGCTGAACGGGCAATGCAGGAAGAGTTCTCCCTGGACCTGAAAGAAGGCTTACTGGAGCCGCTGGCGGAGACGGCACGGCTGGCCATTATCTCGGGGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGACGCGCTGGCCCACGCCACTATCAACATTTTACCCATTACTCAGGGAACTTCCTGAATGTCAATCTCTGTCATGGTAAATAACgatattgccaccaatagcagg                                                            >  2:309720/1‑233 (MQ=255)
                                                                                                                                                                            gAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAAGGAACGGCTGGCAATTATCTCGGTGGGAGGGGATGGTATGGGCACCTTGCGTGGGATCTCGGCGAAATTATTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCCTTGCTCAGGGGTCTTCTGAACGCTCAATCTCTGTCAGGGGAAATAACGATGTGGCGACCACTGGCgcg                                                            >  1:451336/1‑249 (MQ=255)
                                                                                                                                                                                     gggCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGTCGGAACGGCTGGCCATTATCTCGGTGGTAGGAGTTGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAACTTCTGAACGCCCAATCTCTGTCGTGGTAAATGACGATGGTGCAAACACTGGCGTGCGCGTTACt                                                   >  2:119246/1‑251 (MQ=255)
                                                                                                                                                                                       gCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGGGTTACtc                                                  >  1:98430/1‑250 (MQ=255)
                                                                                                                                                                                              gggaagaggctggccTGGAACTGAAAGAAGGCTTACTGGTGCCGCTGGCGGTGACGGAACGGTTGGCCATTATCACGGTGGTAGGTGGTGGTATGCGCACCTTGCGTGGGATGTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAACTTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAACCTCTGCCGTGGTAAATAACGAGGATGCGACCACTGGCGTGCGCGTTACTCATCAGATg                                          <  2:143812/238‑1 (MQ=255)
                                                                                                                                                                                              aGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGCACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATg                                          >  1:238868/1‑251 (MQ=255)
                                                                                                                                                                                               ggAAGAGTTCTACCTGGAACTGAAAGCAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCa                                                                                                      >  1:364290/1‑190 (MQ=255)
                                                                                                                                                                                               ggAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGGAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGc                                         >  1:219127/1‑251 (MQ=255)
                                                                                                                                                                                               ggAAGAGGTCTCCCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCTGTCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCa                                                                                                      <  2:364290/190‑1 (MQ=255)
                                                                                                                                                                                                gAAGAGTTTTACCTGGAACGGAAAGAAGGCTTACTGGAGGCGCGGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCt                                        <  2:151720/251‑1 (MQ=255)
                                                                                                                                                                                                gAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTg                                                                                                                                                                     <  1:18876/127‑1 (MQ=255)
                                                                                                                                                                                                               gAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa                         <  1:224728/251‑1 (MQ=255)
                                                                                                                                                                                                                   tgatagcagGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGCTTGTATGCGCACCTTGCGTGGGTTCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCCACATTGTCGCCATTGCTCAGGGCTCTTCTGAACGCTCATTCTCTGTCGTGTTAAATACCGCTGATGCTACCACTGGCGTGCGCGTTTCTCTTCCGATGCTTTTCAATAcctcttcggtt                     >  2:192260/8‑242 (MQ=255)
                                                                                                                                                                                                                                      gCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTg  >  1:425470/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
CTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTG  >  NC_000913/1270‑1750

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