Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,638,228 A→G 23.5% 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,2280AG23.5% 85.8 / 3.9 55E196G (GAG→GGG) pitAphosphate transporter, low‑affinity; tellurite importer
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base G (9/5);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTT  >  NC_000913/3638001‑3638463
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
ctctCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTTTGATTTTCTTGCTGGCTTGCTCCCGGAGCGGGACCA‑AAAAACGCCCCCCGTTCCCCCTGACCCCACCCGGGCGGGaaacaa                                                                                                                                                                                                                                   >  1:960190/1‑236 (MQ=255)
   tcATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAaaagaaa                                                                                                                                                                                                                                 <  1:388284/238‑1 (MQ=255)
             aTTGGCGCGATCATCGGGTTTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAGTATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGCTTTTTTTGCTTCGTCGCTACTGGGTCGGGCCCACAGAAACGCCCCGGTACCCCCCTTA‑CCCAACGGGGCGGGAAAAGAAAAACgggcaaaa                                                                                                                                                                                                                      >  2:388284/1‑233 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGCTGCACTCAATATCCCGACAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTTTGCTTTTCTTTCTGCGTAGCTACCGGAGCGGCCCCC‑AGAAACGAGCCCGATTCCACCCGACCCCAACGGCCGGTGAAAAGAAAGTAGGCCAGaa                                                                                                                                                                                                                      >  1:920955/1‑239 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAACGAATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCTGCCTGTCGTTTGCTGGCGGTCTGCTTTTCTTGTCTCGCCGTTTCTGTAGCGGCCCAA‑CGAAAACCCCCCGTATCCACCTGACCCCCCCGGGACCTCAAAAAAAAAACCACAAGaa                                                                                                                                                                                                                      >  2:868721/1‑239 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGTTGCGTCGCTACTGGAGCGGCACCA‑AGAAACCCCCCCGGATCCCCCTGACCCCAGCGGGGCGTGGAAACAAAGACGGCAAGaa                                                                                                                                                                                                                      >  1:610860/1‑239 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGGGCGGCACCC‑AGAAACGCGCCCGTATCCACCCGACCCCAACGGGGCGTGAAAAGAAAGACgcacaaaa                                                                                                                                                                                                                      >  1:1012590/1‑232 (MQ=255)
             aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAACGGGCCCC‑AGAAACGCGCCCGTATCCACCTTACCCCAGCGGAGCGTGAAAAGAAAACAGGCAAGaa                                                                                                                                                                                                                      >  1:558317/1‑239 (MQ=255)
                            gggTTTGTTTTACCCAGTGGTTTGAGGCCCGGGGCGGCGGTGTGTGATGCCCTCACTATCCCGAAAGTTTTAATTTTTTTCGGTTCTCTTATCGTTTCCCCGATTGTCGGCCTGTGGTTCGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCCCCTGACCCCAGCGGGGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTg                                                                                                                                                                                                       <  2:1012590/239‑1 (MQ=255)
                                ttGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGAGGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAc                                                                                                                                                                                                    <  1:241893/238‑1 (MQ=255)
                                ttGGTTTAACCAAGGCGTTGATGCCCGGGCCGTCAGTGGGGGATGCACTCAATATCCCGAAATTATTAAGTTTTTTCGTTTCTTTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTCCTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGCCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAc                                                                                                                                                                                                    <  2:719891/238‑1 (MQ=255)
                                 tGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTACTGGAGCGGCACCC‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGGGCGTGGAAAGAAAAACGGGAAGAAAAAACCGCCGTTTCGGAccg                                                                                                                                                                                                  >  1:549952/1‑237 (MQ=255)
                                    tttacccgttGCTTTGGTTACGGGGAGGTGAGGGTTGGATGCACTAATTTGCCCGAATGTATTGAGTATTTTCGGTTTTCTGATCGTTTCCCCTATGGTCGGCCTGGTGTTTGCGGGCGGTCTGATTTTCTTGCTGCGTCGCTCCTGGGGCGGCACCA‑AGAACCGCGCCCGTATCCACCTGCCCCCATCGGAGCGTGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTa                                                                                                                                                                                               <  2:484127/230‑1 (MQ=255)
                                    tttAACCTGTGTGGTTATGACGGGGTCTTCGGTGGTTGATGCTCTCAATAACCCGAATGTATTAAGTATTTTCGTTTCTCTGATCGTTCCCCCTATTGGCGGCCTGGTGTTTGCTGGCGGTCTGATTTCCTTGCTGCGTCGCTACGGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTa                                                                                                                                                                                               <  2:460695/239‑1 (MQ=255)
                                    tttAACCAATGGGTTGATGCCCGGGCCGGTAGTGGTGGATGCACTCAATATCCCGAAAGTTTTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTa                                                                                                                                                                                               <  1:975727/239‑1 (MQ=255)
                                        aCCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGCTGCGTAGCTCCTGGGGAGGCACCC‑AGAAACGCGCCCGGATCCACCGTGCCCCCTCGGAGCGCGAAAAAAAAAACGGCAAGAAAAAACCGCCCTTTTggcagcgcgttgt                                                                                                                                                                                           >  2:566084/1‑228 (MQ=255)
                                        aCCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGGGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCACCGGGGCGTGAAAAGAAAAACGGCAAGAAAAAACCGCCCTTCTGGACGCGGATTg                                                                                                                                                                                            >  1:719891/1‑238 (MQ=255)
                                             tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTTCCCGTCGCTCCTGGCGCGGCACCC‑AGAAACGCGCCCGTATCCACCTCACCCCAGCGGAGGCAGAAAAGAAAGACGGCAAGAAAAAACAGCCGCTCTGGGCGCGGATTgggccg                                                                                                                                                                                       >  2:810793/1‑236 (MQ=255)
                                             tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGGATCCACCTGACCCCCGCGGAGCGTGAAAAGAAAAACGGCAAGAAAAAGCCGCCGTTCTGGACGCGGATTGGGCTGa                                                                                                                                                                                      >  1:502833/1‑239 (MQ=255)
                                                   gatgaCCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCCCTCCTGGGGCGGCCCCA‑AGAAACGCCCCCGTATTCCCCCTACCCCACCGGGACGGGAAAAGAAAAAAGGACAGGAAAAACCACCGCTCCGGACCCGTATTgccataacacct                                                                                                                                                                                 >  2:313096/1‑228 (MQ=255)
                                                          gggACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCCCCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGGGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGGCGCGTATTGCGCTGATCCTTTCCGCTa                                                                                                                                                                          >  1:943735/1‑238 (MQ=255)
                                                                         gATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGGATTGCCCTGATCCTTTCCGCTATCGGCGGGGGGgttt                                                                                                                                                           >  1:262501/1‑238 (MQ=255)
                                                                             cACGCAATACCCGAAAAGTATGAGTTTTTTTCGGTTCTCTGACCGTTTCCCTTATGGTCGGCCTGGTGTTTGCGGGCGGTCTGATTTTCTTGCTGCGGCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGCCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCa                                                                                                                                                       <  1:810793/238‑1 (MQ=255)
                                                                                            aagtattaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGGGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACGCGGATTGCGCTGATCCTTTCCGCTATCGGCGGGGGGGTTTGCCACGGCGGGGACCATGGt                                                                                                                                       >  1:178802/1‑239 (MQ=255)
                                                                                                     gttgttgcGGTTCTTGTATCGTTTCCCTTTTTGCCGGCCTGGTGTTTGCGGGCGGTTTAATTTTCTGTCGGGCTCCGTACGGGGGCGGCACAA‑AGAACCGCGCCCGTATCCCCCTGTCCCCAGCGGAGCGGGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGc                                                                                                                              <  2:943735/232‑1 (MQ=25)
                                                                                                           ccGGTTTTCTGATCGTTTCCCCTTTTGTCGGCCTGGTTTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGt                                                                                                                        <  1:727954/238‑1 (MQ=255)
                                                                                                                tctcTGCTCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTTTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCCCCTGCCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGt                                                                                                                   <  1:611450/239‑1 (MQ=255)
                                                                                                                  tctGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTTCGCTATCGGCGTGGCGGTTTTCCACGGGGCGCACCGTGGTCGAGAAGGCAATTGGCTGGttt                                                                                                                 >  1:1000873/1‑238 (MQ=255)
                                                                                                                     gtctttttCCCTTATTGGCGGCTGGGGTTTTGGGGGCGGTCCTTTTTTCTTGGCGCGTCGCTCGGGGGGGGCCACCA‑AAAAACGCCCCGTTATCCACCTGCCCCGAGCGGAGCGTGAAAAGAAAGACGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATg                                                                                                               <  2:107938/233‑1 (MQ=25)
                                                                                                                                     ggcGGCCGGGTGTTTGCGGGCGGTTGGATTTTCTGGCTGCGTCGCTACGGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGGGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGt                                                                                              <  1:168397/237‑1 (MQ=255)
                                                                                                                                     gTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCCCTATCGGCGTGGGGGTTTCGCCCGGCGCGAACGAAGGTGCGAAAGGGATTGGTGTGGTTTTGTTGGGattgatctgtt                                                                                               >  1:107938/1‑233 (MQ=255)
                                                                                                                                              tttttttgtGGCGGTCTGTTTTTTTTGTTGGCTCGCTACGGGGGCGGCCCCA‑AGAAACGCCCCCGTATCCCCCTGACCCCAGCGGAGGGTAAAAAGAAAGACGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCcagc                                                                                     <  2:1000873/231‑1 (MQ=37)
                                                                                                                                                         cGGTTTGATTTCCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCgtg                                                                          <  1:755303/238‑1 (MQ=255)
                                                                                                                                                                          tgcGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCCGAAAGGCa                                                                                                                             >  2:613669/1‑171 (MQ=255)
                                                                                                                                                                          tgcGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCa                                                                                                                             <  1:613669/171‑1 (MQ=255)
                                                                                                                                                                             cgcgcgACGGGAGGGGCCACA‑AAAAACCCCCCCGGACCCACCGGACCCCAGGGGAGCGTGAAAAGAAAGACGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTgg                                                      <  1:424197/236‑1 (MQ=37)
                                                                                                                                                                               ggctcgtGGGGCGGCAACA‑AAAACCCCGCCCGTATCCCCCTGACCCGAGGGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCt                                                    <  1:190554/232‑1 (MQ=38)
                                                                                                                                                                                      ggAGCGGCACCA‑AGAAACCCGCCGGTACCCCCCTGACCCCAGCGGGGGGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATc                                            <  1:359585/239‑1 (MQ=255)
                                                                                                                                                                                         gCGGCCCCA‑AGAAACGCGCCCTGACCCCCCCGCCCCCAGCGGGGGGTGAAAAAAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCac                                          <  2:786059/238‑1 (MQ=38)
                                                                                                                                                                                            gCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATAacccgt                                      >  1:215960/1‑239 (MQ=255)
                                                                                                                                                                                            gCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCAAGGCGCGAACGATGGTCAGAAAAGCCTTGGTCTGGTTATGTTGGtc                                                                                                         >  1:151456/1‑172 (MQ=255)
                                                                                                                                                                                             cACCA‑AGAAACGCGCCCGTACCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgt                                      <  1:168677/238‑1 (MQ=255)
                                                                                                                                                                                                        cgcgCCCGTATCCACCTTACCCCAGCGTTCCGCGAAACGAACGCCGGCAAGAAAAAGCCGCCGTTCTGGGCGCGAATTGCGCTGATCCTTTCCGCTATCGGCGTAGCGTTTTCCCCCGAGGCGAGCGGTGGCCCGCACGGCAGTGGGCTCGCTTTTTTGGCATTGTTTTGCGTCCCGCCCCCAGGCC‑TCGCGGGGCACAGTAATTCCCCTGGGTAGGCAACTacccggacccgtcacac                           >  2:159604/1‑234 (MQ=255)
                                                                                                                                                                                                        cgcgCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCCCGGCGCGAACGATGGTCAAAAAGGCCTTGGTTGGGTTATTTTGGGATTGCTTGGCTGCGCGCCCTCAGGGTGTCGTGGTGAACA‑GAATGCCCCCTGCCACGAAAAAaccccgaccccggggg                            >  2:736203/1‑232 (MQ=255)
                                                                                                                                                                                                         ccccccTTACCCCCCTGCCCCGAGCGGGGCGAGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGcc                          <  1:789286/236‑1 (MQ=37)
                                                                                                                                                                                                          cgcCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGGATTGATTTGCTTGCCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCCCTGGCTCGGAAATCACCCTTACCCGGGGTGcc                          >  1:206589/1‑238 (MQ=38)
                                                                                                                                                                                                               ttATCCCCCTGCCCCCAGGGGGGCGGGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATcaac                    <  1:843975/238‑1 (MQ=255)
                                                                                                                                                                                                                   ccACCTGACCCCAGGGGAGCTGGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACg                <  1:580230/239‑1 (MQ=255)
                                                                                                                                                                                                                     aCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGGCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCCCTGGCTACGAAATCACCCGTACCCCGGATGCCATcaacactgt               >  1:401915/1‑234 (MQ=255)
                                                                                                                                                                                                                       cTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGGTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGCAATCACCCGTACCCGGGGTGCCACCAAAAACGTCgc            >  1:845229/1‑238 (MQ=255)
                                                                                                                                                                                                                          aCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGGCTGGGTATGTTGGTATTGATTGGCGTCGCGCCCGCAGGCT‑TCGTGGGGAAACTGAATGCCACTGGCTACGGAATAACCCGTACCCCTGATGCCCTCCACAACGTcagagg         >  1:76309/1‑234 (MQ=37)
                                                                                                                                                                                                                          aCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGAATTGGTCTGGTTATGTTGGGATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGGGAACATGAATGCCACTGGCCACCAAATACCCCGTACCCGTGATGCCACCCACAAACTCGaa           >  1:423306/1‑237 (MQ=38)
                                                                                                                                                                                                                           cccgcgCGGAGCGTGAAAAAAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTACTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGc         <  1:670521/233‑1 (MQ=255)
                                                                                                                                                                                                                           ccccAGCGGAGGGTAAAAAGAAAGAGGGAAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGc         <  1:313096/238‑1 (MQ=255)
                                                                                                                                                                                                                                  ggAGCGTAAAAAAAAAGCCGGCAAAAAAAGGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACttt  <  1:135619/238‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
CTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCA‑AGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCT‑TCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTT  >  NC_000913/3638001‑3638463

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: ‑