Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,782,912 A→C 19.9% 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,9120AC19.9% 112.3 / 5.8 77V250G (GTG→GGG) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base C (6/11);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  NC_000913/3782678‑3783149
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
aGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGTGACTGGTTGTCGGTACAGGTAAGCACCAg                                                                                                                                                                                                                                                 <  2:329126/239‑1 (MQ=255)
       tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGGCTTGGTGTCGGGTAAGGTAAACACCAAATCGCCa                                                                                                                                                                                                                                          >  2:857832/1‑239 (MQ=255)
       tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCCTCACTGGTTGTCGGGACAGCGAAGCACCAAAACGCCa                                                                                                                                                                                                                                          >  1:917818/1‑239 (MQ=255)
       tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGGTGTCGGTGAAGGGAAGGACCCAAACGCCa                                                                                                                                                                                                                                          >  2:574333/1‑239 (MQ=255)
           tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTTTCCGTACCGGGGAGCACCaacacgcccagc                                                                                                                                                                                                                                        >  2:433040/1‑227 (MQ=255)
                 cTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCa                                                                                                                                                                                                                                <  2:1125230/239‑1 (MQ=255)
                  tcgcggcctttttttCCGAATTATAATCCTAGTTAAATTCCCTCTGTTGTTGGCATTTCAACGCGAACGCGATGCGCTCGTTGGCGGACTCCTTCCGTATCGCGCTAGCTTCCCACCACCAGTCCAATCTTCTCCTGCCCCCTTTGTGCATCAAGGCCAGGACCAACCGTCAGCCCAAAACGGCGGTTACGCGACTGGTTGTCGGACAAGGTAAGCACCAGACCGCCAAGCCCCGCCAt                                                                                                                                                                                                                               <  1:699977/230‑1 (MQ=255)
                  ccgcgcgcgTTTTTTCCTCAATATAATACTGAAGAAATTTCCTCGTTTATTGGCATTTCACCGCCGAAGCGATGCGCAGGTCCGCGGACTTCTTTCTTATTGCGTTGGCCTCCCACCCCCTGACCAATCTTCTCCGGCACGTTTGTTCCATCCATGCCCTGACCGACCTTCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGTTACAGGTAACCCCCAGACCGCCAAGCCCCGCCAt                                                                                                                                                                                                                               <  1:121969/238‑1 (MQ=255)
                     cgcgcgTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCGGGTTTCGGTAAAGGTCAGCaaccttcccttataccaatccgca                                                                                                                                                                                                                               >  1:1081090/1‑214 (MQ=255)
                          gTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTGTCGGTACGGGGTAGCACCAGATCCCCAAGCCCCCCCCTGCCCCTaa                                                                                                                                                                                                                      >  2:460549/1‑239 (MQ=255)
                          gTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTCCACCACCTGCCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGTTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATaa                                                                                                                                                                                                                      <  2:1029668/238‑1 (MQ=255)
                                 ccGCAATATAATACTTGTTAAATTCCCTGGGTTATTGGATTTCCAACGCCGAAGCGATGCGCAAGTTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGCCAAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGACCGCAAAGCCCCGCCATGCCCATAAAGGTGGc                                                                                                                                                                                                             <  1:940548/239‑1 (MQ=255)
                                      atatAATCCTTGTTAAATTTCCTCGTTTTTGGGCATTCCAACGCCGAAGCGATGCGCCGGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGACCGCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGt                                                                                                                                                                                                        <  1:471277/239‑1 (MQ=255)
                                             gctggATAAATTCCCGCGGTTTTTGGCATTCCAACGCCGAACCGATGCGCCAGTTCGCGGCCTTCTTTCTTTTTGCGGTACCCTTCCTCCCCTTGCCCAATTGTCTCCTGCGGGCTTTGTACAGCCATGCCGTGACCGAGCATGTTGCCAAAGCGGGGGTTACGCGACTGGTTGTCGGTCCAGGAAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                                <  1:969787/235‑1 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGCCTCGTTGTCGGGACCGGGAAGCCCCAAAACACCAAACCCCCCCCCTCCCACAAAGGTGGGCAGGGGCGCAcccca                                                                                                                                                                                             >  1:143559/1‑236 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAAATCGCCAAGCCCCCCCATGCCCATAAAAGGGGCAGGGGCGGCACCCCgc                                                                                                                                                                                           >  2:981734/1‑239 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGGAAGCACCAAATCGCCAAGCCCCCCCCTACCCCTAAAAGTGGCAGGGGCGGCACCCAgc                                                                                                                                                                                           >  2:172493/1‑239 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCACTGGGTGTCGGGACAGGGAAACACCAAAACGCCACACCCCCCCCTGCCCATACAGGCGGGAGGGTGGGCAccccg                                                                                                                                                                                             >  1:1190638/1‑236 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTCGGGACAGGTAATGACCAAATCCCCCACACACCCCCAGACCACAAAGCAGGCAGGGTCGGGAccacac                                                                                                                                                                                            >  1:570572/1‑235 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTCGGGACAGGGAAGCACCAAAGCGCCCAACCCCCCCATGCCCCTAAAAGGGGGAGGGGCGGCAccccg                                                                                                                                                                                             >  2:703181/1‑236 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTCACGCGACTTGTGGTAGGTGCACGGCAACACCACAATCCCCAAACCCCCCCAGCACCTAAAACGGGCAGGGTGGCGCAcccca                                                                                                                                                                                             >  1:842522/1‑237 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGCGACCCCAATGGGTGTCCGGTAGGAGAAGCACCCCATCGCCCAACCCCGCCCCACCCATAAAAGACGCCGGggcggccgccac                                                                                                                                                                                             >  1:763724/1‑232 (MQ=255)
                                                      aTTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTTCGCGACTGGTTGTCGGTACAGGTAAGACCCACATCGCCCAACCCCGCCCTGCCCCTAAAAGTGGGCGGGGCGGCACCCCAcgccg                                                                                                                                                                                       >  2:1177600/1‑238 (MQ=255)
                                                         tCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAGCACCAAATCGCCAAGCCCCCCCATGCCCCTAAAGGTGGCAGGGGCGGGACCCACCGCCGCAc                                                                                                                                                                                    >  2:831939/1‑238 (MQ=255)
                                                                ttATTGGCATTTCAACCCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGa                                                                                                                                                                               <  2:1155246/236‑1 (MQ=255)
                                                                ttATTGGCATTGCAATCCCGAGGCTATCCGCCGGTGCGCGGACTTCGTTCCTATTGCGGTCGCCTCCCACCACCTGACCAATCTTCCCCGGCGCGCTTTGTACACCCTGGCCCGGACCGAGCACCATGCCAAACCGGCGGTTACGCGACTGTTTTTCGTTACAGGTAAGCCCCAGACCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGAcgc                                                                                                                                                                            <  1:685783/239‑1 (MQ=255)
                                                                       cattcaaaCGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATGGCGGTAGCCTTCCACCCCCGGACCAATCTTCTCCTGCGCGCTTTGTACACCCATGCCCTGACCGACCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACAtt                                                                                                                                                                      <  1:981734/232‑1 (MQ=255)
                                                                           tCAACGCCGAAGCGATGCGCCAGTTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAATCTTCCCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                                 <  2:906319/239‑1 (MQ=255)
                                                                           tCAACGCCGAAGCGATGCGCCAGTGCGCGGACTTCTTTCGTATTGCGGTAGCCTCCCACCCCCTGCCCAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTGGGTACAGGTAAGCACCAGAGCGCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                                 <  2:260333/239‑1 (MQ=255)
                                                                           tCAACCCCGAGGCGATGCGCCTGTTCGCGGCCTTCTTTCGTATTGCGGTACCCTTCCCCCCCCTGCCCAATCTTCCCCTGCGCGCTTTGTACATCCATGCCCTGACGGAGCATCAGGCCAAAAGGGCGGTTCCGCGACTGGTGGTCGGTACAGGTAAGCCCCAGACCGCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                                 <  1:845969/239‑1 (MQ=255)
                                                                           gcAACGCCGAAGCTATGCGCCCTTTCGCGGACTTCTTTCGTATTGCGTTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                                 <  2:657570/238‑1 (MQ=255)
                                                                              aCGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACACCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagccag                                                                                                                                                              <  2:917818/239‑1 (MQ=255)
                                                                                 ccGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGCCTGGGTGTCGGTCCAGGGAAGCACCAGAAGCCCAAGCCCCGCCATGCCCCTAAAGGTGGCAGGGGCGGCACCCACGGCCGCCCCCAGACGGCACACTTccacccacc                                                                                                                                                            >  1:568568/1‑230 (MQ=255)
                                                                                  cGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGATACCCCGCCGGCTTGCGGGACAGGGGTACCACAAAAAGCCAAACCCCGCCACGCCCAAAAAAGGAGCAGGGGGGGAGCCCCGCGCGCCCCCCAGAcgcgcaacctcacacaccac                                                                                                                                                           >  1:527562/1‑222 (MQ=255)
                                                                                    aaGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAAATCGCCAAACCCCCCCCTGCCCCTAAAAGGGGCAGGGGCGGCACCCCGCGCCGCACCCAGACGCGACCTTTCACCCAAcccccc                                                                                                                                                        >  2:1153798/1‑236 (MQ=255)
                                                                                    aaGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGAATGGGTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGGGGCAGGGGCGGCACCCAGCGCCGCACCAAGAAGCGGCAATTCAGCCAACCCACg                                                                                                                                                        >  2:959680/1‑239 (MQ=255)
                                                                                               cAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGCCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcg                                                                                                                                             <  2:884660/239‑1 (MQ=255)
                                                                                                                        cGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGCCTGGGTGTCGGTACAGGTAAGCACCACAAAGCCCAACCCCGCCATGCCCCTAAAGGGGGCAGGGGCGGGACCCAGCGCCGCCCCAACAAGCGCCAATTTCGCCCGCCCACCGGGGATCAGACCCGCACGCGCACTCGCACCAAAAc                                                                                                                     >  1:310229/1‑238 (MQ=255)
                                                                                                                            acccTCCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATa                                                                                                               <  2:310229/237‑1 (MQ=255)
                                                                                                                                   caccacCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTCCCCCCCTGGTTGTCGGTAACGGTAAGCCCCCGATCGCCAAGCCCCCCCCTGCCCCTAAAGGTGGCAGGGGCGGCACCCACCGCCGCCCCCCGCCCCGACACTTCAGCCAGCCCCCGGGCGGTCTGCGCCCCCCCCCCATTCGCACCACAACCGATACCGTCg                                                                                                         >  1:593136/1‑238 (MQ=255)
                                                                                                                                    accaccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGCCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa                                                                                                       <  2:171357/239‑1 (MQ=255)
                                                                                                                                     cccccTGCCCAATCTTCCCCTGCCCGCTTTGTCCACCCATGCCCTCACCGAGCACCACCCCAAACCGGCGGTCACGCGCTTGGTTGTCGTTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGAc                                                                                                      <  1:959680/236‑1 (MQ=255)
                                                                                                                                     ccacgttcccAATTTTTCCCGGCGCGTTTGTTAAATCCAGGCCCCGCCCGAGCATAATGCCAAAACGGCGTTTCCGCGACTGTTTGTCGGTACAGGTAAGCACCAAACCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGAc                                                                                                      <  1:901853/231‑1 (MQ=255)
                                                                                                                                          ggcccATTCTTCCCCTGCGCCCTTTGTACACCCATGCCTGGACCGACCACCAGCCAAAAACGGCGGTTCCGCGCCGGGTTGCGGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATccc                                                                                                 <  1:318426/236‑1 (MQ=255)
                                                                                                                                           gccagaTCTTCTCCGGCGCGTTTTGTACACCCATGCCCTGACCGAGCATCAGGCAAAAACGGCGGTTCCGCGACGGGTTGCCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATccc                                                                                                 <  2:1106378/233‑1 (MQ=255)
                                                                                                                                                 tcttttCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCaaa                                                                                                                       <  2:747130/211‑1 (MQ=255)
                                                                                                                                                 tcttctCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCCACTGGTTGTCGGGACCGGTAAACACCAGCTCGCCAAGCCCCCCCATTCCCCTTAAGCCGGCAACGTGGGCACCCCGACGCGCAACCACTACCCGCCATTCCACTACCCCACGCGCGGTCAGCCCCCCACACGCATCCGCCCCCAAACCGCAACACCCGGACATCACCccacca                                                                                          >  1:1188485/1‑239 (MQ=255)
                                                                                                                                                 tcttctCCTGCGCGCTTTGTACACCCAGGCCCTGCCCGAGCATCATGCCAAAACGGCGTTTACGCGCCTGTTTGCCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGa                                                                                                                  <  2:944407/216‑1 (MQ=255)
                                                                                                                                                      tcCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAAGTAAAACCCCAAACGCCAAGCCCCCCCCCGCCCCATAAAGGGGCAGGGGCGGCCCCCCCCGCCCCACCCAGACGCGCCATTTCAGCCAGCCCCCGGGGGAGCAGCGCCCCACCCGCACTCGCCCCCACACAGATAcccccggccacaccccccccaccac                                                                                      >  1:1044119/1‑215 (MQ=255)
                                                                                                                                                        cTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGCCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCaa                                                                                   <  2:568568/239‑1 (MQ=255)
                                                                                                                                                          gcgcgcTTTGTCCACCCTTGCCCTGACCGAGCATCATGCCAAAACGGCGTTTACGCGTCTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATa                                                                                 <  2:963057/239‑1 (MQ=255)
                                                                                                                                                           cgcgcTTTGTACATCCATGCCCTGACCGAGCATGATGCCAAAACGGCGGTTACGCGACCGGGTGTCGGTACAGGCAAGCACCAGAACGCCCAACCCCGCCCTCCCCCCAAAAGGGGCAGGGGCGGCACCCACCGCCCCACCAACAAGCCACATTTTCTCCAACCCACGGGGGATCAGCGCCCCACGCCCACTCGCCCCCCAACCACGACCGCCGCACATCCCCCCACCACTAGCacaa                                                                                 >  1:1018894/1‑235 (MQ=255)
                                                                                                                                                           cgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTTTCGGTACAGGTAAACACCAAAAAGCCCAACCCCCCCACGCCCCTAAAGGTGGCAGGGGCCGCACCCAGCGCCGCCCCCACAACCGCCACTTTAGCCACCCCACCGGGGGGCAGCGCCGCACGCGCGCTCGCACCCACACCCATACCGCCCGCCATCCCCGCCCCCATGGcgctc                                                                                 >  1:1106378/1‑234 (MQ=255)
                                                                                                                                                           cgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTTGGGAAAGGGAAGCACCAAATCGCCCAGCCCCCCCCCGCCCCTAAAAGGGGCGGGGGCGGCACCCCGCGCCCCCCCACGACGCGCCATTTTAGCCAGCCCCCGGGGGATCAGCGCCGTACGCGCCGTCGCCCCCAAACCCCTAACGCCGGGAATCCCCccaccaccgcccata                                                                                 >  1:995549/1‑229 (MQ=255)
                                                                                                                                                                       aTCCATGCCCTGCCCGAGCATCATGCAAAAGCGGCGTTTACTCGACTGGTTGTCGGTACGGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCGCcg                                                                    <  2:520/239‑1 (MQ=255)
                                                                                                                                                                         ccATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                                 <  2:167100/145‑1 (MQ=255)
                                                                                                                                                                         ccATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg                                                                  <  2:905556/239‑1 (MQ=255)
                                                                                                                                                                         ccATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg                                                                  <  2:6381/239‑1 (MQ=255)
                                                                                                                                                                            tGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAAGGAACCACCAGAACCCCCAACCCCCCCATGCCCCTAAAAGGGGGCGGGGCGGCACCCACCGCCCCCCCCCGACGCGACAATTATACCACCCCACGCGGGGGCACGGCCCGCAGCGCCCTTGCCCCAAAACCAATAACCCACGGCATCCCCGCCCCACAGCGAAGAACGCTTTTCACCCcgcc                                                                >  1:815248/1‑238 (MQ=255)
                                                                                                                                                                              cccTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAACCACCAGAGCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCgccgcc                                                             <  2:914199/239‑1 (MQ=255)
                                                                                                                                                                              cccTGACCGAGCATCATGCCAAAACGGCGGCTTAGCCGCTGGTTGTCGGGAAAGGTAAACCCCAAAACGCCAAACCCCGCCCTTCCCCTAACAGGGGGCGGGGCCGGCCCCAGCGCCGCCACCAGACGCGACCTTTCATCCAACCCCCAGGGGAACACCGCCAGAACCGCCATCCCACCAAAACACACACCGGCCGCCAAACCCGCACCACCGTGGATAACGTATTTCACCCcgccgc                                                              >  1:157676/1‑238 (MQ=255)
                                                                                                                                                                               ccTGCCCGAGCATCAGGCCAAACCGGCGGTTGCGCGACTGTTTGGCGGTACGGGTAAGCACAAGCCCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCa                                                            <  1:851479/239‑1 (MQ=255)
                                                                                                                                                                                 tGACCGAGCATAATGCCAAACCGGCGTTGCCGCGACTGGTTGTCGGTACAGGAAAGCCCCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCCCCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATCCCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAg                                                          <  2:842522/239‑1 (MQ=255)
                                                                                                                                                                                   aCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTTTCGGGACAGGGAAGCACCAAAACGCCCAACCCCCCCACGCCCCTAAAAGGGGGCGGGTCGGCCCCCCGCGCCGCACCCAGACGCGACATTTCAGCCCGCCCACGGGGGATCAGCGCCGGTCGCGCATTCGCCCCCAAACCGGTACCGCCGGGCATCCCCCCCCCCAAGGCCATAAAGTTTTTCTCCCCgccgcccaccg                                                        >  1:572207/1‑234 (MQ=255)
                                                                                                                                                                                    gcGAGCACCAGGCCAAACCGGCGGTCACGCGGTTGGTTGTGGGTCCGGGAAGGCACAAGACCGCCAGCCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                       <  1:831939/238‑1 (MQ=255)
                                                                                                                                                                                      gacaaccatGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCa                                                     <  2:1190638/233‑1 (MQ=255)
                                                                                                                                                                                      gAGCATCATGCCAAAACGGCGTTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCa                                                     <  2:143559/239‑1 (MQ=255)
                                                                                                                                                                                         catcatGCCAAACCGGCGGTTACGCGACTGGTTGGCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa                                                                                                       <  2:331642/186‑1 (MQ=255)
                                                                                                                                                                                         catcatGCCAAAACGGCGGTTACGCGGCTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGc                                                  <  2:681111/239‑1 (MQ=255)
                                                                                                                                                                                         caccaggCCAAAACGGGGGTCCCGCGATGGGTTGTCGGTACGGGTAAGCCCCAGATCGCCAGCCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGc                                                  <  1:172493/233‑1 (MQ=255)
                                                                                                                                                                                                ccAAAACGGGGGTTACGCGACTGGTGGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATg                                             <  2:1044119/237‑1 (MQ=255)
                                                                                                                                                                                                         cggTTACGCGGCTGGTTGTCGGTACAAGGAAGCACCAAAACGCCAAGCCCCCCCATGCCCCTAAAAGGGGGAGGGTCGGCCCCCAGCGCCCCCCCAAGACGCGACATTTCAGCCAGCCCACGGGGGATCAACGCCGTACGCCCATTCGCCCCAAAACCGATACCGTCGGGCAACCCCCCCCCAATGGCAATAACGTTTTTCACCCCCCCGCCAAGCCGCCCGCCACAGCAAACCGGAtt                                  >  1:824376/1‑239 (MQ=255)
                                                                                                                                                                                                         cggTTACGCGACTGGTGGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGAt                                   <  2:824376/238‑1 (MQ=255)
                                                                                                                                                                                                                          tgccGGTAAAGATAAGCCCCAGCCCGCCCGCCCCCGCCATGCCCATAAGGGGGGCAGGGTGTGCGCCCAGCGCCGCACCAACACGCGACATGTCAGCCAGCCCACGGGTGATCAGCTCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGaa                   <  1:1043258/235‑1 (MQ=255)
                                                                                                                                                                                                                               gTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtttt             >  2:563985/1‑239 (MQ=255)
                                                                                                                                                                                                                               gTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtttt             >  2:1077633/1‑238 (MQ=255)
                                                                                                                                                                                                                               gTACAGGTAAGCACCAGATCGCCAAGCACCGCCATGCCCATAAAGGTGGCCGGGTCGGCACCCAGCGCCGCACCTAGACGCGCCCTTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGCCGGACATCCCCGCACCAACGGCAATAACGTTTTTCAACGCGCCGCCAAGCTGCAGCCCAATGAAATCCGCATTGCTgccaaaccggaacgtttt             >  1:707642/1‑222 (MQ=255)
                                                                                                                                                                                                                                          cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca  >  2:613472/1‑239 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
AGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  NC_000913/3782678‑3783149

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: