Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,782,913 A→C 22.0% V250G (GTG→GGG)  gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,782,9130AC22.0% 96.7 / 10.0 73V250G (GTG→GGG) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base C (3/14);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 2.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGC  >  NC_000913/3782678‑3783149
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
tAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCa                                                                                                                                                                                                                                              <  2:1162383/239‑1 (MQ=255)
        tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCCCCACCGGGTGTCGGGTCCGGGAAGCACCAGaaccccc                                                                                                                                                                                                                                      >  1:1037945/1‑233 (MQ=255)
        tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGCCCGGTTGTCGGTTCCGGGAAGCACCACATCGccc                                                                                                                                                                                                                                     >  1:505674/1‑238 (MQ=255)
        tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGCCCGGGTTTTGGTACCGGTAAACACCCGATTGccc                                                                                                                                                                                                                                     >  1:616706/1‑238 (MQ=255)
         cAATGCTGCGTGGCGGTCTTTTTTCCCGAAATATAAACCTTGAAAATTTCCCTGGTTTATTGGCATTCCAACGCCGAGCCGTGCCCCCGTTTCGCGGCTTTCTTTCGTATTGCGGTACCCTCCCCCCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACACCCAGGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGcc                                                                                                                                                                                                                                      <  2:616706/237‑1 (MQ=255)
         cAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCa                                                                                                                                                                                                                                     <  1:812018/238‑1 (MQ=255)
            tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCACTGGTTGTCGGGACAGGTAAGAACCCAAACGCCaaac                                                                                                                                                                                                                                 >  1:484379/1‑236 (MQ=255)
                   ccggccgtgTTTTTCCCGCAATATAATCCTTGAAAAATTTCCTGTGTTTTTGGCATTTCAGCGCCGAACCGATGCCCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGGGCATCATGCCAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                         <  2:944065/231‑1 (MQ=255)
                        cgcgTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTTAGCGCCCGGGTTTTGGGGACGGGAAAGCCCCAAAAGCCCAAACCCCCCCATGccc                                                                                                                                                                                                                     >  2:147341/1‑239 (MQ=255)
                           gTTTTTCCCGAAATAAAATCCTTAAAAATTTTCCTCGGTTATTGGCATTTCAACTCCGAAGCGATGCCCCAGTTCGCGGCCTTCTTTCGTATGGCGGTAGCCTTCCCCCCCCGGACCAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGCCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATaa                                                                                                                                                                                                                 <  2:1037945/239‑1 (MQ=255)
                                       atatAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGGAAGAACCAGATAGCCCAGCCCCCCCCCGCCCCTACAAGTGTGAGGGt                                                                                                                                                                                                     >  2:155776/1‑239 (MQ=255)
                                        tataATACTTGAAAATTTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                         <  1:1051276/218‑1 (MQ=255)
                                        tagaataCTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                         <  1:536292/215‑1 (MQ=255)
                                         ataatcCTTGATAATTTCCCCCGTTTTTTGACATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTACCCTTCCCCCCCCTGCCCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCg                                                                                                                                                                                                  <  2:484379/239‑1 (MQ=255)
                                            acccTTGATAATTTCCCATGGTTTTTGGCTTTACAACGCCGAAGCGATGCCCCTGTTCGCGGACTTCTTTCGTTTTGCGGTAGCCTTCCCCCACCGGCCCAATCTTCCCCGGCGCGTTTTGTACCTCCATCCCCTGCCCGGGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAAACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCa                                                                                                                                                                                               <  2:669049/236‑1 (MQ=255)
                                              tcTTGATAATTTCCCTCGGTTGTTGCCATTTCACCCCCGAAGCGATCCCCCAGTTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAATCTTCCCCTGCGCGCTTTGTACACCCATGCCCTGCCCGACCATCATGCCAAACCGGCGGTTCCGCGCTTGTTTGTCGGTACAGGTAAGCCCCAGATCCCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                             <  2:176200/238‑1 (MQ=255)
                                              aCTAGATAATTTGCCTCGGTTTATGGCGTTTCACCGCAGAACCGATGCGCCAGTTGGCGGCCTTCTTTCGTATTGCGGTACCCTCCCCCCCCCTGCCCAATCGTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGCCCAAAACGGCGGTTACGCAACTGGTTGTCGGTACGGGTAAGCACCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                             <  2:176547/239‑1 (MQ=255)
                                                ttGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCACTGGGTGTCGGTACAGGGAAGCACCCGATCGCCAAGACCCCCCATTCCCCCAAAGGGGGCAGGGTCGGCAccc                                                                                                                                                                                            >  1:837181/1‑238 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGCCTGGTTGTCGGGTAAGGGTAGAACCCAAACCCCACACCCCCCCCATCCCCCAAAAGGGGCAGGGTCGGGAccccc                                                                                                                                                                                          >  1:1022929/1‑236 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGCCTGGTTGTCGGGACAGGGAAGCCCCCAAACGCCAAACCCCGCCATGCCCCTACAAGTGGCAGGGGCGGCACCCac                                                                                                                                                                                          >  1:425071/1‑237 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTTTCGGTACAGGGAAGCACCAGATCGCCCAACCCCGCCCTGCCCCTAAAGGGGGACGGGGCGGGAcccca                                                                                                                                                                                          >  1:176547/1‑236 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAACCCCGCCCATCCCCTAAAGGTGGCAGGGTCGGCACCCAgc                                                                                                                                                                                         >  1:858266/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCTGGTGTCGGGACAGGTAAGCACCAGAAAGCCAAACCCCCCCATGCCCCTAAAGGGGGCCAGGGCGGGACCCAgc                                                                                                                                                                                         >  1:632862/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCGGTTGTCGGTACAGGGAAGCACCAAATCGCCAAGCCCCGCCACGCCCCCAAAAGGGGCAGGGTCGGCACCCAgc                                                                                                                                                                                         >  1:318431/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTCGGGTCAGGTAAGGACCAGAACTCCCAACCCCGCCATGCCCCCCAAGAAGGCAAGGTCGGccgcccc                                                                                                                                                                                          >  2:468665/1‑232 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTCGGGACAGGGAAACCACAAAACGCCAAACCCCCCCATGCCCATAAAAGTGGCAGGGTCGGCACCCCgc                                                                                                                                                                                         >  1:283221/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGACTTGTTTTCGGGACAAGGAAGCACCCGAATGCCAAACCCCCCCCAGGCCACTAAAGTGGCAGGGGCGGCACCCac                                                                                                                                                                                          >  1:669049/1‑237 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTACGCGACTGGGTGTCGGTACAGGTAAGAACCAAAACGCCAAGACCCGCCATGGGCCTAAAAGTGGCAgggggcgggccccc                                                                                                                                                                                           >  2:365181/1‑227 (MQ=255)
                                                   aaaaaTTTCCGGGGTTATTGCCATTCAAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCAGGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAgcg                                                                                                                                                                                        <  2:425071/237‑1 (MQ=255)
                                                         ttCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCg                                                                                                                                                                                                                             <  1:802272/197‑1 (MQ=255)
                                                                  tATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGAcgcg                                                                                                                                                                         <  1:365181/239‑1 (MQ=255)
                                                                           ttCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGCCTGGTTGTCGGGGCAGGGAAGGACCAGAAGGCCCAACCCCCCCCAGCCCCCTAAAGGGCGGAGGGTGGGGCCCCCAGCGCCCCCCAAAGCcgcacctact                                                                                                                                                                   >  2:1133558/1‑231 (MQ=255)
                                                                            tCACCGCCGAAGCGATCCGCAATTCCCCGGCTTTCTGCGGTATTGCGGTAGCCTTCCCCCACCGGACCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCGGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTAAAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                               <  2:505674/239‑1 (MQ=255)
                                                                                       gCGATGCCCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACgg                                                                                                                                                     <  1:159351/238‑1 (MQ=255)
                                                                                                     tCGGGGACTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcg                                                                                                                                           <  1:1241212/235‑1 (MQ=255)
                                                                                                     tCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcgc                                                                                                                                          <  1:434305/236‑1 (MQ=255)
                                                                                                         ggATTTCGTGCGTATTGGGGTAGCCTTCCACCCCCTGACTAATCTTCCCCTGCGCGTTTTGAACTTCCAGGCCCGGACCGAGCATCAGGCCAAAACGGCGGTTACGCGACGGGTTGTCGGTACAGGTAAGCACCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTAcg                                                                                                                                    <  2:862122/238‑1 (MQ=255)
                                                                                                                    gTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCAAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc                                                                                                                        <  1:468665/239‑1 (MQ=255)
                                                                                                                     tATTGCGGTAGCCTTCCCCCCCCTGACCAATCTTCCCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATaaa                                                                                                                                                                                                                <  2:488724/151‑1 (MQ=255)
                                                                                                                     tATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTATCCCCCGGGTGTCGGTACAGGTAAGCACCAAAACGCCAAACCCCCCCCTGCCCCTCAAAGGGGCCGGGGCGGGACCCAGCGCCGCACCCACACGCGACCTTTTATCCCGCCCACCGGCGGTCAAGGCCGATCgcgccgtcgtcccc                                                                                                                       >  2:1010800/1‑229 (MQ=255)
                                                                                                                           ggTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                             <  1:330607/163‑1 (MQ=255)
                                                                                                                               gCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTTTCGGGACAGGGAAGCACCAGAACGCCAAACCCCCCCCTGCCCCTAAAAGGGGCAggg                                                                                                                                                                                                      >  1:488724/1‑151 (MQ=255)
                                                                                                                                ccTTCCACCACCGGACCAATCTTCCCCGGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATg                                                                                                                                                                                                                        <  2:259762/132‑1 (MQ=255)
                                                                                                                                    ccacccCCTGCCCAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGGCGGTACAGGTAAGCACCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcg                                                                                                                                           <  1:438757/204‑1 (MQ=255)
                                                                                                                                         cccTGACCAATCTGCCCCTGCGCGCTTTGTACATCAAGGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACGGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa                                                                                                      <  1:529693/235‑1 (MQ=255)
                                                                                                                                          ccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGGAAGCACCAGATCGCCAAGCCCCCCCCTGCCCCTAAAgg                                                                                                                                                                                                              >  1:259762/1‑132 (MQ=255)
                                                                                                                                            tGACCAATCTTCCCCTGCGCCTTTTGTCCATCCATGCCCTGACCGAGCATCATCCCAAACCGGCGTTTCCGCGACTTGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCTCCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATccc                                                                                                <  1:995961/239‑1 (MQ=255)
                                                                                                                                              cccACTCTCCACCTGCGCGCTTTGTCCATCCATGCCCTAACCGGGCTGCATGCCAAAGCGGCGGTCCCGCGACTGTTTTCCGGTACGGGTAAGCCCCAGACCCCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCCCCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCg                                                                                              <  2:318431/238‑1 (MQ=255)
                                                                                                                                              aCCAATCTGCCCCGGGCCGTTTTGTACATCCACTCCCTACCCGACCATCATGCAAAAGCGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCg                                                                                              <  2:285254/239‑1 (MQ=255)
                                                                                                                                                 aaTCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGACCTGTTTTCGGGACAGGTAAACACCAAAACGCCAAGCCCCGCCCTGCCCCTAAAAGGGGCACGGGCGGGCCCCCGCGCCCCCCCCACACGCGGAACTTCAGCCAGACCACGCGCGCTCAACCCCCGACGCGCACTCGCACCAAAACACAAACCGTCGGGCagtcccaca                                                                                            >  2:687014/1‑230 (MQ=255)
                                                                                                                                                   ttttctCCTGCGCGTTTTGTCCATCCATGCCCTGCCCGACCATCATGCAAAACCGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc                                                                                                                        <  1:824167/206‑1 (MQ=255)
                                                                                                                                                   tcttctCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGGCAAAAGCGGGGGTTACGCGACTGGTTGTGGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                               <  1:588800/169‑1 (MQ=255)
                                                                                                                                                                         aTCCATGCCCTGACCGAGCATCAGGCCAAAACGGCGGTTACGCGACGGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcg                                                                   <  1:1133558/239‑1 (MQ=255)
                                                                                                                                                                           ccATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg                                                                 <  1:888145/239‑1 (MQ=255)
                                                                                                                                                                           ccATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGa                                                                                                                                                                             <  1:989866/131‑1 (MQ=255)
                                                                                                                                                                                 ccTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCa                                                           <  1:147341/239‑1 (MQ=255)
                                                                                                                                                                                     aCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAGCACCAAAACGCCCAGCCCCCCCATGCCCCTAAAGGGGGCCGGGGCGGCCCCCCGCCCCCCCCCCAAACGCGGCAttt                                                                                                                                                                   >  2:989866/1‑131 (MQ=255)
                                                                                                                                                                                     aCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGTAAACACCAAAACGCCCAACCCCCCCCCGCCCCTAAAAGGGGGAGGGGCGGCACCCCGCGCCCCCCCAAGAACCCACCTTTTAGCCCCCCCCCCGGGGAGCAGCGCCGGACGCGCCTTCGCACCCAAACCGATACCCCCGGACCTCCCCGCCCCCCTGGCCACAAAGTATTTTATCGCgccgccccgc                                                        >  2:886736/1‑234 (MQ=255)
                                                                                                                                                                                      cggAGCATAACGCCAAACCGGCGGTTCCGCGCCTGTTGGCGGTTACAGGTAACCCCCAGCCCGCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                      <  2:858266/237‑1 (MQ=255)
                                                                                                                                                                                       cGAGCATCATGCCAAAACGGGGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                      <  1:1150249/238‑1 (MQ=255)
                                                                                                                                                                                         aGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTGTCGGGACAGGTAAGCACCAGAACGCCAAGCCCCCCCCTGCCCCTAAAGGGGGCAGGGTCGGCACCCAACGCCCCCCCCAGACGCGACATTTCAGCCAGCCCCCGGGTGATCAGCGCCCTACCCGCCCTCGCACCCAAACAGCCAACCTCACAAATCCCCTCCACAATAGCCATAACGGATTTTACAGCGCCGCCaccagtca                                                    >  2:847472/1‑231 (MQ=255)
                                                                                                                                                                                          gCATCATGCAAAAACGGCGGTTACGGGACTGTTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                      <  1:1176027/235‑1 (MQ=255)
                                                                                                                                                                                           catcatGCCAAAGCGGCGGTTCCGCGACTGGTTGCCGGTACAGGAAAGCCCCAGACCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACg                                                  <  1:513776/238‑1 (MQ=255)
                                                                                                                                                                                               atGCCAAACGGGGGTTTCCGCGCCTGGTTGGCGGTACAGGTAAGCCCCAGATCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAt                                             <  1:155776/239‑1 (MQ=255)
                                                                                                                                                                                                         ggccttcccgcgTCTGGTTGTCGGAACAGGAAAGCACCAGACCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCCCCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGa                                   <  2:283221/231‑1 (MQ=255)
                                                                                                                                                                                                                          gTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGa                   >  1:5143/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGcc         >  1:1141603/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGc          >  1:1228092/1‑238 (MQ=255)
                                                                                                                                                                                                                                      gggTAAGCCCCAGACCCCCCAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg        <  1:1075470/238‑1 (MQ=255)
                                                                                                                                                                                                                                      aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCgccgcc                                                            >  2:1054638/1‑187 (MQ=255)
                                                                                                                                                                                                                                      aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCgccgcc                                                            <  1:1054638/187‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGCCATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGt    >  1:542455/1‑238 (MQ=255)
                                                                                                                                                                                                                                             cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgc  >  2:486184/1‑238 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
TAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGC  >  NC_000913/3782678‑3783149

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