Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,638,227 A→G 22.4% E196G (GAG→GGG)  pitA → phosphate transporter, low‑affinity, tellurite importer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,638,2270AG22.4% 111.4 / 4.2 67E196G (GAG→GGG) pitAphosphate transporter, low‑affinity, tellurite importer
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (11/6);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAG  >  NC_000913/3638000‑3638465
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
ctctCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTTTGATTTTCTTGCTGGCTTGCTCCCGGAGCGGGACCAAAAAACGCCCCCCGTTCCCCCTGACCCCACCCGGGCGGGaaacaa                                                                                                                                                                                                                                       >  2:960190/1‑236 (MQ=255)
   tcATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAaaagaaa                                                                                                                                                                                                                                     <  2:388284/238‑1 (MQ=255)
         cgtgtttggccgctttctcGGGATTGGTTTAACCAATGCTTTTGTGACCGGGACCTCAGTGGTTGATGCACTCAATTTCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCa                                                                                                                                                                                                                              <  1:525071/222‑1 (MQ=255)
             aTTGGCGCGATCATCGGGTTTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAGTATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGCTTTTTTTGCTTCGTCGCTACTGGGTCGGGCCCACAGAAACGCCCCGGTACCCCCCTTACCCAACGGGGCGGGAAAAGAAAAACgggcaaaa                                                                                                                                                                                                                          >  1:388284/1‑233 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGCTGCACTCAATATCCCGACAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTTTGCTTTTCTTTCTGCGTAGCTACCGGAGCGGCCCCCAGAAACGAGCCCGATTCCACCCGACCCCAACGGCCGGTGAAAAGAAAGTAGGCCAGaa                                                                                                                                                                                                                          >  2:920955/1‑239 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAACGAATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCTGCCTGTCGTTTGCTGGCGGTCTGCTTTTCTTGTCTCGCCGTTTCTGTAGCGGCCCAACGAAAACCCCCCGTATCCACCTGACCCCCCCGGGACCTCAAAAAAAAAACCACAAGaa                                                                                                                                                                                                                          >  1:868721/1‑239 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGTTGCGTCGCTACTGGAGCGGCACCAAGAAACCCCCCCGGATCCCCCTGACCCCAGCGGGGCGTGGAAACAAAGACGGCAAGaa                                                                                                                                                                                                                          >  2:610860/1‑239 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGGGCGGCACCCAGAAACGCGCCCGTATCCACCCGACCCCAACGGGGCGTGAAAAGAAAGACgcacaaaa                                                                                                                                                                                                                          >  2:1012590/1‑232 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAACGGGCCCCAGAAACGCGCCCGTATCCACCTTACCCCAGCGGAGCGTGAAAAGAAAACAGGCAAGaa                                                                                                                                                                                                                          >  2:558317/1‑239 (MQ=255)
                            gggTTTGTTTTACCCAGTGGTTTGAGGCCCGGGGCGGCGGTGTGTGATGCCCTCACTATCCCGAAAGTTTTAATTTTTTTCGGTTCTCTTATCGTTTCCCCGATTGTCGGCCTGTGGTTCGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGGGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTg                                                                                                                                                                                                           <  1:1012590/239‑1 (MQ=255)
                                ttGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGAGGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAc                                                                                                                                                                                                        <  2:241893/238‑1 (MQ=255)
                                ttGGTTTAACCAAGGCGTTGATGCCCGGGCCGTCAGTGGGGGATGCACTCAATATCCCGAAATTATTAAGTTTTTTCGTTTCTTTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTCCTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGCCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAc                                                                                                                                                                                                        <  1:719891/238‑1 (MQ=255)
                                 tGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTACTGGAGCGGCACCCAGAAACGCGCCCGTATCCACCTGACCCCAGCGGGGCGTGGAAAGAAAAACGGGAAGAAAAAACCGCCGTTTCGGAccg                                                                                                                                                                                                      >  2:549952/1‑237 (MQ=255)
                                    tttacccgttGCTTTGGTTACGGGGAGGTGAGGGTTGGATGCACTAATTTGCCCGAATGTATTGAGTATTTTCGGTTTTCTGATCGTTTCCCCTATGGTCGGCCTGGTGTTTGCGGGCGGTCTGATTTTCTTGCTGCGTCGCTCCTGGGGCGGCACCAAGAACCGCGCCCGTATCCACCTGCCCCCATCGGAGCGTGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTa                                                                                                                                                                                                   <  1:484127/230‑1 (MQ=255)
                                    tttAACCTGTGTGGTTATGACGGGGTCTTCGGTGGTTGATGCTCTCAATAACCCGAATGTATTAAGTATTTTCGTTTCTCTGATCGTTCCCCCTATTGGCGGCCTGGTGTTTGCTGGCGGTCTGATTTCCTTGCTGCGTCGCTACGGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTa                                                                                                                                                                                                   <  1:460695/239‑1 (MQ=255)
                                    tttAACCAATGGGTTGATGCCCGGGCCGGTAGTGGTGGATGCACTCAATATCCCGAAAGTTTTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTa                                                                                                                                                                                                   <  2:975727/239‑1 (MQ=255)
                                        aCCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGCTGCGTAGCTCCTGGGGAGGCACCCAGAAACGCGCCCGGATCCACCGTGCCCCCTCGGAGCGCGAAAAAAAAAACGGCAAGAAAAAACCGCCCTTTTggcagcgcgttgt                                                                                                                                                                                               >  1:566084/1‑228 (MQ=255)
                                        aCCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGGGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCACCGGGGCGTGAAAAGAAAAACGGCAAGAAAAAACCGCCCTTCTGGACGCGGATTg                                                                                                                                                                                                >  2:719891/1‑238 (MQ=255)
                                             tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTTCTGGCTCGCTACTGGGGCGGCCCCAAGAAACGCGCCCGTAACCACCCGCCCCCAGCGGGGCGGGAAAAGAAAACAGGCAAAAAAAAACCGCCCCTCTgggcacggtttgggggga                                                                                                                                                                                          >  1:611450/1‑223 (MQ=255)
                                             tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTTCCCGTCGCTCCTGGCGCGGCACCCAGAAACGCGCCCGTATCCACCTCACCCCAGCGGAGGCAGAAAAGAAAGACGGCAAGAAAAAACAGCCGCTCTGGGCGCGGATTgggccg                                                                                                                                                                                           >  1:810793/1‑236 (MQ=255)
                                             tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGGATCCACCTGACCCCCGCGGAGCGTGAAAAGAAAAACGGCAAGAAAAAGCCGCCGTTCTGGACGCGGATTGGGCTGa                                                                                                                                                                                          >  2:502833/1‑239 (MQ=255)
                                                   gatgaCCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCCCTCCTGGGGCGGCCCCAAGAAACGCCCCCGTATTCCCCCTACCCCACCGGGACGGGAAAAGAAAAAAGGACAGGAAAAACCACCGCTCCGGACCCGTATTgccataacacct                                                                                                                                                                                     >  1:313096/1‑228 (MQ=255)
                                                     ttgcgggggcttgagggggggATGCCCTCAATATCCGAAAAGTTTTAAGTTTTTTGGGTTCTCTGATCGTTTCCCCTATTTCCGGCTGGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGCCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCTTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTc                                                                                                                                                                                   <  1:46374/219‑1 (MQ=255)
                                                          gggACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCCCCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGGGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGGCGCGTATTGCGCTGATCCTTTCCGCTa                                                                                                                                                                              >  2:943735/1‑238 (MQ=255)
                                                                         gATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTACCGGGCCGGCCCCACGAAAACCACCCGTATCCACCTGACCCCCCCGGGGCGGCAAAACAAAGACGGCACGAAAAAGCCGCCGGTTCGGACGGGCATTTCGCTGATACttttcccccttacccgggcgtgttt                                                                                                                                                              >  1:757999/1‑217 (MQ=255)
                                                                         gATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGGATTGCCCTGATCCTTTCCGCTATCGGCGGGGGGgttt                                                                                                                                                               >  2:262501/1‑238 (MQ=255)
                                                                             cACGCAATACCCGAAAAGTATGAGTTTTTTTCGGTTCTCTGACCGTTTCCCTTATGGTCGGCCTGGTGTTTGCGGGCGGTCTGATTTTCTTGCTGCGGCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGCCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCa                                                                                                                                                           <  2:810793/238‑1 (MQ=255)
                                                                                          gAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGCGGGTTTCgaccgggggcgaccagg                                                                                                                                             >  2:786059/1‑223 (MQ=255)
                                                                                            aagtattaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGGGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACGCGGATTGCGCTGATCCTTTCCGCTATCGGCGGGGGGGTTTGCCACGGCGGGGACCATGGt                                                                                                                                           >  2:178802/1‑239 (MQ=255)
                                                                                                     gttgttgcGGTTCTTGTATCGTTTCCCTTTTTGCCGGCCTGGTGTTTGCGGGCGGTTTAATTTTCTGTCGGGCTCCGTACGGGGGCGGCACAAAGAACCGCGCCCGTATCCCCCTGTCCCCAGCGGAGCGGGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGc                                                                                                                                  <  1:943735/232‑1 (MQ=25)
                                                                                                           ccGGTTTTCTGATCGTTTCCCCTTTTGTCGGCCTGGTTTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGt                                                                                                                            <  2:727954/238‑1 (MQ=255)
                                                                                                                tctcTGCTCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTTTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGCCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGt                                                                                                                       <  2:611450/239‑1 (MQ=255)
                                                                                                                  tctGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTTCGCTATCGGCGTGGCGGTTTTCCACGGGGCGCACCGTGGTCGAGAAGGCAATTGGCTGGttt                                                                                                                     >  2:1000873/1‑238 (MQ=255)
                                                                                                                     gtctttttCCCTTATTGGCGGCTGGGGTTTTGGGGGCGGTCCTTTTTTCTTGGCGCGTCGCTCGGGGGGGGCCACCAAAAAACGCCCCGTTATCCACCTGCCCCGAGCGGAGCGTGAAAAGAAAGACGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATg                                                                                                                   <  1:107938/233‑1 (MQ=25)
                                                                                                                                     ggcGGCCGGGTGTTTGCGGGCGGTTGGATTTTCTGGCTGCGTCGCTACGGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGGGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGt                                                                                                  <  2:168397/237‑1 (MQ=255)
                                                                                                                                     gTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCCCTATCGGCGTGGGGGTTTCGCCCGGCGCGAACGAAGGTGCGAAAGGGATTGGTGTGGTTTTGTTGGGattgatctgtt                                                                                                   >  2:107938/1‑233 (MQ=255)
                                                                                                                                              tttttttgtGGCGGTCTGTTTTTTTTGTTGGCTCGCTACGGGGGCGGCCCCAAGAAACGCCCCCGTATCCCCCTGACCCCAGCGGAGGGTAAAAAGAAAGACGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCcagc                                                                                         <  1:1000873/231‑1 (MQ=37)
                                                                                                                                                         cGGTTTGATTTCCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCgtg                                                                             <  2:755303/238‑1 (MQ=255)
                                                                                                                                                                        gcggcgtcgccgggggggCGGCCCCAAGAAACCCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCa                                                             <  2:536532/222‑1 (MQ=255)
                                                                                                                                                                          tgcGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTTGCCCGGCGCCAACGATGGGCAGCAAGGCATTTGGCTGGTTTTGGTGGGCTTGATTGGTGGCGGCCCACCAGGCTTctgggggaaaagcaaacccccc                                                            >  1:95384/1‑218 (MQ=255)
                                                                                                                                                                          tgcGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCCGAAAGGCa                                                                                                                                 >  1:613669/1‑171 (MQ=255)
                                                                                                                                                                          tgcGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCa                                                                                                                                 <  2:613669/171‑1 (MQ=255)
                                                                                                                                                                             cgcgcgACGGGAGGGGCCACAAAAAACCCCCCCGGACCCACCGGACCCCAGGGGAGCGTGAAAAGAAAGACGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTgg                                                         <  2:424197/236‑1 (MQ=37)
                                                                                                                                                                              tCGCTACTGGAGCTGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACAGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTACCGGCGTGGCGCTTTCGCACGCCCCGCACCATGGTCCGACAAGCATTTGCGCGGGTATGTTGGTATTGATTGGCGTCGGCCCAGCAACCCTCGCTGGTGaaacagatgacaaaggc                                                         >  1:388291/1‑223 (MQ=255)
                                                                                                                                                                              tCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGGGGACATgctggttcttttaa                                                       >  2:121452/1‑226 (MQ=255)
                                                                                                                                                                               ggctcgtGGGGCGGCAACAAAAACCCCGCCCGTATCCCCCTGACCCGAGGGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCt                                                       <  2:190554/232‑1 (MQ=38)
                                                                                                                                                                                      ggAGCGGCACCAAGAAACCCGCCGGTACCCCCCTGACCCCAGCGGGGGGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATc                                               <  2:359585/239‑1 (MQ=255)
                                                                                                                                                                                         gCGGCCCCAAGAAACGCGCCCTGACCCCCCCGCCCCCAGCGGGGGGTGAAAAAAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCac                                             <  1:786059/238‑1 (MQ=38)
                                                                                                                                                                                            gCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATAacccgt                                         >  2:215960/1‑239 (MQ=255)
                                                                                                                                                                                            gCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCAAGGCGCGAACGATGGTCAGAAAAGCCTTGGTCTGGTTATGTTGGtc                                                                                                             >  2:151456/1‑172 (MQ=255)
                                                                                                                                                                                             cACCAAGAAACGCGCCCGTACCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgt                                         <  2:168677/238‑1 (MQ=255)
                                                                                                                                                                                                      acccccccccttccccccTGCCCCCAGGGGGGGGTGAAAAGAAAGGCGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGa                                 <  1:401915/222‑1 (MQ=37)
                                                                                                                                                                                                        cgcgCCCGTATCCACCTTACCCCAGCGTTCCGCGAAACGAACGCCGGCAAGAAAAAGCCGCCGTTCTGGGCGCGAATTGCGCTGATCCTTTCCGCTATCGGCGTAGCGTTTTCCCCCGAGGCGAGCGGTGGCCCGCACGGCAGTGGGCTCGCTTTTTTGGCATTGTTTTGCGTCCCGCCCCCAGGCCTCGCGGGGCACAGTAATTCCCCTGGGTAGGCAACTacccggacccgtcacac                              >  1:159604/1‑234 (MQ=255)
                                                                                                                                                                                                        cgcgCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCCCGGCGCGAACGATGGTCAAAAAGGCCTTGGTTGGGTTATTTTGGGATTGCTTGGCTGCGCGCCCTCAGGGTGTCGTGGTGAACAGAATGCCCCCTGCCACGAAAAAaccccgaccccggggg                               >  1:736203/1‑232 (MQ=255)
                                                                                                                                                                                                         ccccccTTACCCCCCTGCCCCGAGCGGGGCGAGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGcc                             <  2:789286/236‑1 (MQ=37)
                                                                                                                                                                                                          cgcCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGGATTGATTTGCTTGCCGCCAGCAGGCTTCGTGGTGAACATGAATGCCCCTGGCTCGGAAATCACCCTTACCCGGGGTGcc                             >  2:206589/1‑238 (MQ=38)
                                                                                                                                                                                                               ttATCCCCCTGCCCCCAGGGGGGCGGGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATcaac                       <  2:843975/238‑1 (MQ=255)
                                                                                                                                                                                                                   ccACCTGACCCCAGGGGAGCTGGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACg                   <  2:580230/239‑1 (MQ=255)
                                                                                                                                                                                                                     ccctgccccagggggAGGGTGAAAAAAAAGCCGGCAAGAAAAGTCCGCCGTTCTGGACGCGTAGTGCGCTGATCCTTCCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTc                 <  2:479902/226‑1 (MQ=37)
                                                                                                                                                                                                                     aCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGGCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCCCTGGCTACGAAATCACCCGTACCCCGGATGCCATcaacactgt                  >  2:401915/1‑234 (MQ=255)
                                                                                                                                                                                                                       cTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGGTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGCAATCACCCGTACCCGGGGTGCCACCAAAAACGTCgc               >  2:845229/1‑238 (MQ=255)
                                                                                                                                                                                                                         gACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGCTGGTCAAAAAGGCATTGGTCTGGTTATGTTGGGATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGGGAACATTAATTGCCCTGGCTCCGGAAACCCCCGGCCCCCGGATGcccccaaacaagcccaaa             >  1:34042/1‑223 (MQ=25)
                                                                                                                                                                                                                          aCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGGCTGGGTATGTTGGTATTGATTGGCGTCGCGCCCGCAGGCTTCGTGGGGAAACTGAATGCCACTGGCTACGGAATAACCCGTACCCCTGATGCCCTCCACAACGTcagagg            >  2:76309/1‑234 (MQ=37)
                                                                                                                                                                                                                          aCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGAATTGGTCTGGTTATGTTGGGATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGGGAACATGAATGCCACTGGCCACCAAATACCCCGTACCCGTGATGCCACCCACAAACTCGaa              >  2:423306/1‑237 (MQ=38)
                                                                                                                                                                                                                           cccgcgCGGAGCGTGAAAAAAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTACTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGc            <  2:670521/233‑1 (MQ=255)
                                                                                                                                                                                                                           ccccAGCGGAGGGTAAAAAGAAAGAGGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGc            <  2:313096/238‑1 (MQ=255)
                                                                                                                                                                                                                                  ggAGCGTAAAAAAAAAGCCGGCAAAAAAAGGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACttt     <  2:135619/238‑1 (MQ=255)
                                                                                                                                                                                                                                    aGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGGCCGAGCGATTGTCATAACGGCCTTGTTCTGGTTATGTTGGGATTGATTGGCGGTGCGCCATCAGACTTCGTGGGGAACATTGATGCAACCGGGTACGAAATCCCCCCGACCCCCGAAGCCATCAACAACgcccgaccattcttctac  >  1:705272/1‑222 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
CTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAG  >  NC_000913/3638000‑3638465

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: