Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,638,227 A→G 22.0% 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.0% 81.8 / 3.0 56E196G (GAG→GGG) pitAphosphate transporter, low‑affinity, tellurite importer
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (5/9);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGT  >  NC_000913/3638013‑3638449
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTTCCCCGTTCCCACTTGATCGGCACCAAGAAAAGCGCCCCTACCCCCCTGACCCCCACGGGGCGGGAaaaggaa                                                                                                                                                                                                                        >  1:463877/1‑228 (MQ=255)
  tGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCGATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGGGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTCCTGGGGGGGGCACAAGAAACGCGCCCGTATCCACCCGACCCCCCCGGGGGGGGAAaaaaaaga                                                                                                                                                                                                                      >  1:619236/1‑228 (MQ=255)
                        taaacaaagggtgtgatgaCGGGGCGGGTAGTGGGGGATCCCCTCAAATCCCCTGAAGTATTAAGTTTTTCCGGTTCTTGGATCTTTTCCCCGTGTGGCGGCCGGGGGTTTGCTGGGGTGCGAATTTTCTTGCTGCGTCGCTACGGGGGCGGCACCAAGAAGCGCGCCCGTACCCACCGGACCCAAGGGGGGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTc                                                                                                                                                                                               <  1:397977/215‑1 (MQ=255)
                              aatggtttaatggCCGGGACGGTGGGGGTGGATGCCCTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAc                                                                                                                                                                                          <  2:338532/218‑1 (MQ=255)
                              aagggtgtgatgaCCGGGACGGCCGGGGTGGAAGCACTCAATTTCCCGAAAGTATTAAGTATTTCCGGTTCCCTGATCGTTTCCCCTATTGTCGGCCGGTTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAc                                                                                                                                                                                          <  2:110958/220‑1 (MQ=255)
                                tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGGAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGCTTCGTCGCTACTGGGGCGGCCCCAAGAAACGGGCCCCTATTCACCTGACCCCACCGGAGCGTGCAAAGAAAGACGGCAAGAAAAAACCGCCCTTCtcggccc                                                                                                                                                                                        >  2:34107/1‑221 (MQ=255)
                                tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGGTTTTCTTGCTGCCTCGCTCCTTGCGCGGCACCAAGAAACTCGCCCCGTTCCCCCCGCCCCCAACGGAGCGGCGAAAAAAAAACGGGAAGAAAAAaccacccctccgggccgg                                                                                                                                                                                       >  1:601671/1‑216 (MQ=255)
                                tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCCGCGGGGGGTGAAAAGAAAGACGGCAAGACAAAGCCTCCCtt                                                                                                                                                                                                >  2:397977/1‑219 (MQ=255)
                                     tgatgaCCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTTCTGCGTCGCTACTTGGACGGCACCAAAAAAAGCGCCCTTTTTCACCTGCCCCCCGCGGCGCGTGAAAAACAAAACGGACAGAAAAAACCGCCGGTATggggccgaattt                                                                                                                                                                                  >  1:502264/1‑218 (MQ=255)
                                            cGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCCCTACTGGAGCGGCACCAAGAAACGCGCCCGGATCCACCCGCACCCAGCGGAGCGTTGAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACCCGCAtttggcggat                                                                                                                                                                           >  2:206869/1‑220 (MQ=255)
                                                 cGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTTTTGCCGCGTTGCTACTGGGGGGGGACCAAGGAAAGCGCCCGGTTTCCCCCTGCACCAACGGAGCGTGCAAAAAAAAACGGCCAGAAAAACCCCGCCGACTGGCCGCGGATTGCGCTGAGCCtgt                                                                                                                                                                      >  1:473422/1‑226 (MQ=255)
                                                  ttCAGTGGTATATGTACTCCATTCCACGAAAGTATTAGTTATTTGCGTTTCTCTGATCGTCTCCCCTGTTGTCGGCCTGTTTTTTGCTGGCGGTCTGATTTTCTTGCTGCGGCCCATCTGGGCCGGCACCAAGAAACGCGCCGGTTTCCACCTGACCCCAGCGGGGCGTGAAAAAAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCttgc                                                                                                                                                                     <  1:256084/227‑3 (MQ=255)
                                                        agggccagCACTCAATACCCGGAAAGTATTAAGGATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCGTGGGGTTTGCTGGCGGTCGGATTTTCTTGCTGCGTCGCTACTGGAGCGGACCCAAGAAACGCGCCCGTACCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTa                                                                                                                                                                <  2:60072/220‑1 (MQ=255)
                                                                       tatCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGGATCCACCTGACCCCAGCGGAGCGGGAAAAGAAAGAAGGCAAGAAAAAACCGCCGTTCTGGACGCGACATGCGCTGATCCTTTCGGCTATCGGCGTGGGGgtttt                                                                                                                                              >  2:578429/1‑227 (MQ=255)
                                                                          cccGAAAGTTTTAAGGTTTTCGGGTTCTCTGCTCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCGTGCTGCGTCGCTACGGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCa                                                                                                                                           <  2:247596/228‑1 (MQ=255)
                                                                           cgggatgttttgatttttttcggttttcggACGGTTTCCCCTGTGTCGGGCCGGGTGTTTGCTGGCGGTCTTATTTGTTTGATGCGTCGCTGCTGGGGCGGCACCAAAAAACGCGCCCGTATCCACCTGACCCGAGCGGAGTGTAAAAAGAAAAACGGCAAGAAAACGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCAc                                                                                                                                          <  1:367613/212‑1 (MQ=255)
                                                                            cGAAAGTATTAATTATTTTCGGTTCTCTGATCGTTTCCCCTATTGCCGGCCTGGTGTTTGCTGGGGGTCTGATTTTCTTGCTGCGTCGCTACGGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCAc                                                                                                                                          <  2:619236/227‑1 (MQ=255)
                                                                               aagtattaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTTGTTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAACAAACCGCCGTTCTGGACGCGTCTTGAGCTGAACCTTGCCGCTATCGGCGGGGCGTTTTTGccccgac                                                                                                                                      >  1:445393/1‑222 (MQ=255)
                                                                               aagtattaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGGGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTCTTGGCCGCGTTTTGCGCTGATCCTTTCCGCAATCGGCGTTGCGGTTTTGCACcgcg                                                                                                                                      >  1:143293/1‑228 (MQ=255)
                                                                                    ttaagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGTCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTAATGGTGCGGCACCAAGAAACGCGCCCGTTTTCCCCTGACCCCCGCGGAGCGTGATAAGCAAGAAGGCCAGAAAAACACGCCCCTCTGGACCCCGATTGCCCTGGTTCTTTCCGCCACTCGGCGGGCCTtttttccccgcgcggac                                                                                                                                  >  1:420012/1‑215 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTTCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATACACCTGACCCCAGCGGAGGGTGAAAAGAAAGACGGCCAGAAAAAACCGCCGTTCTCGACCCGGATTGCGCTGTTCACTTCCGCCATCGGAGTGGGCGTTTTccccggccgcgaccactgtca                                                                                                                          >  1:409137/1‑208 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGGGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAAACGGCCAGAAAAAGCCCCCGTTCTGGACGCGGCTTTGGCCGAACCCTTCCGCCATCGGCCTGGGGGTTTGCGCCCGGGGGGAACGAAggcg                                                                                                                         >  1:132693/1‑226 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACCCGCATTGCGCTGATCCTTTCCCCTATCGGCGTCGCGTTTTCTCACCGGGCGGACACTGGTcc                                                                                                                        >  2:488332/1‑227 (MQ=255)
                                                                                            tttCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACCTGGGCGGCCCCCAGAAACGCGCCCGTTTCCCCCTTACCCCCCCGGGGCGTGAAAAGAAAGACGGCCAGAAAAAACCCGCGGTCTGGACACCGATTGCGCGTATCCTTTCCCCCATCGCGGTGGGGTTTTTGccgcgcgcccacaagggtca                                                                                                                         >  1:588956/1‑209 (MQ=255)
                                                                                              tcggttctcgttgtgTTTCCCTTAGTGTCGGCCTGGTTGTGGCGGGCGGTTCGGTTTGCGTGCTGCTTGCTTACTGGGGCGGCACAAATGAACGCGCCTGGTCCCACCTGCCCCCAGAGGAGCGTAAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGa                                                                                                                     <  1:509163/214‑1 (MQ=19)
                                                                                                   tctcTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCCGCGGAGCGTGGAAAGAAAGACGGCAAGAAAAAGCCGCCCGTCTGGACGCGTATTGCGCTGATCCTTTCCCCTATAGGCGTGGCGTTTTCTGACgggccggaacaggtgcaaaaaagc                                                                                                                  >  1:457517/1‑206 (MQ=255)
                                                                                                    tctggaccttttCCCCTGTGTGCGGCCGGGTGTTTGCTGGCGGTCTGATTTTCGTGCTGCGTCGCTACGGGAGCGGCACAAAGAAACGCGCCGGTATCCACCTGACCCAAGGGGAGCGGGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCa                                                                                                               <  1:488332/219‑1 (MQ=37)
                                                                                                      ctGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTAACCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAc                                                                                                                                                                                          <  2:38057/155‑1 (MQ=255)
                                                                                                        gtttgTTTCCCCTTTTGTGGGCGTGGTGTTTGCGGGCGGTTGTATTTTCGTGCTGGGCCGCTAGGGGGGCGGCCACAAAAAACGCGCCCGTATCCACCGGCCCCCAGCGGGGGTTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTgg                                                                                                           <  1:614542/224‑1 (MQ=25)
                                                                                                         tttgTTCCCCCTTTTGTCGGGGTTGTGTTTGCTGGCGGTCTGATTTTGTTGCTGCGGCGCTACTGGAGCGGCACAAAGAAACGCGCCCGTACCCCCCTGCCCCCAGCGGAGTGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGt                                                                                                          <  1:540536/225‑1 (MQ=38)
                                                                                                         tatgttttccctgtttccgggCCGGGTGTTTGCTGGCTGTCTGTTTTTTTTGCTGCGTCCCTACTGGAGCGGCACCAAAAAACGCTCCCGTACCCACCTGACCCCGGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTgg                                                                                                           <  1:131728/209‑1 (MQ=37)
                                                                                                         aTCTTTTCCCCGTTTGGCGGCCTGGTTTTGGGTGGCGGCTGGATTTCTTTGCGCGGTCGCTCCTGGGGGGGCACCAGAAAACGGGCCGATATCCACCGGACCCCAGGGGGGGTGGAAAAGAAAAACGGCAAGAAAAGGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGt                                                                                                          <  1:206869/228‑1 (MQ=19)
                                                                                                          ttgTTTCCCCTATTGGCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGGGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGt                                                                                                          <  2:558355/225‑1 (MQ=255)
                                                                                                           atgtgtataagagacgtgcccTGGTGTTTGCTGTCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTc                                                                                                         <  2:154064/209‑1 (MQ=255)
                                                                                                               tCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAACCCCCGGTCTGGACGCCGATTGCGCTGATCCTTTCCGCCATCGGCGTGGGGGTTTTCCCCGGCGCGACCGACGGGCAGAAAAGCacttgggtggt                                                                                                     >  1:558355/1‑217 (MQ=255)
                                                                                                                            ggCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGGGTTTTTTCACGGCCCGCACGATGGTCCGAAAAGGCTTTGGCTGGGTTTGTTGGGattt                                                                                        >  2:204376/1‑227 (MQ=255)
                                                                                                                                        gCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCCAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGGGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTCGACGCGGATTGCGCTGATCCTTTCCCCTATCGGCGTGGGCTTTTTGCCCGGGGCGCACGGAGGGTAAGAAAGCAAGGGTCGGGGTCTGTTTGTTATGGTTGGCgcccgg                                                                            >  1:569551/1‑223 (MQ=255)
                                                                                                                                              gggttgtttttCCTGTTGCGTCGCTATGGGAGCGGCCACAAGAAGCGCGCCCGTCTCCCCCTGACCCGAGGGGGGGGTGAAAAGAAAGACGGCAAGAAAAGGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCcagc                                                                       <  1:418567/220‑1 (MQ=37)
                                                                                                                                                     gtttgtgggcgcgcggacaggggAGGGGCACCAAGAAACGGCCCGGAAACCAGCTGACCCCAGCGGAGGGAAAAAAGAAGGTAGTGAAGAAAAGGCCGCCGGTCGGGACGAGTATTGCGCTGAGCCTTGCCGGGATCGGCGGGGCGTTTTCGCATGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCg                                                               <  1:461648/207‑1 (MQ=255)
                                                                                                                                                           ttctggggtgCTACGGGGGCGGCAACAAGAAAGGCGCCCTGAACCACCTGACCCCAGGGGGCGTTGAAAAGAAATACGGCAAGAAAAAGGCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCgtggtg                                                          <  1:34107/218‑1 (MQ=37)
                                                                                                                                                                         tGGAGCGGTCCCAAGAAACGCCGCCTGTCCCACTCGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCTTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACt                                           <  2:454136/228‑1 (MQ=37)
                                                                                                                                                                             ggcgaaacAAAAACCGCCCCCGTAACCCCCGGACCCCGGGGGGCGGGAAAAAGAAAGACGGCAAAAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCt                                       <  1:204376/221‑1 (MQ=25)
                                                                                                                                                                              gggCACAAAGAAACGCCCCCGTAACCCCCTGACCCCGGCGGAGGTGGAAAAGAAAGGCGGCAAGAAAAGGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTAGGGGAGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTa                                      <  2:588956/227‑1 (MQ=255)
                                                                                                                                                                               ggCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCCTCCTGGTGAACATGAATGCCACTGGCTAc                                     >  2:119258/1‑228 (MQ=255)
                                                                                                                                                                                  aCCAAGAAACGCGCCCGTAACCACCGGACCCCAGCGGAGCGTGAAAAAAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGaa                                  <  2:445393/228‑1 (MQ=255)
                                                                                                                                                                                     aaGAAACGCGCCCGTATCCACCTGGCCCCGGGGGGGTGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATc                               <  2:601671/228‑1 (MQ=255)
                                                                                                                                                                                     aaGAAACGCGCCCGTATCCACCTGCCCCCAGGGGAGCGTGAAAAAAAAGACGGCAAGAAAAATCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATc                               <  2:497587/228‑1 (MQ=255)
                                                                                                                                                                                      aaaaaaggaccccGGATCCAACTGGCCGCATGGGAGCGGGAAAAGAAAGACGGAAAAAAAAAGCCGCCGTTCGGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATc                               <  2:596543/218‑1 (MQ=37)
                                                                                                                                                                                             ccccgcggtttacctgtggccccCGCGGAGGGTAAAAAGAAAGACGGCAAAAAAAAGCCGCCGTTCCGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgt                         <  1:622560/209‑1 (MQ=37)
                                                                                                                                                                                              cccccGTACCCCCCGAACCCAGGGGGGGCGTGAAAAAAAAGCCGGAAAGAAAAGGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgtac                       <  2:420012/226‑1 (MQ=37)
                                                                                                                                                                                                       tCCACCTGACCCGAGCGGAGCGTGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATg               <  2:132693/228‑1 (MQ=255)
                                                                                                                                                                                                          aCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCCGGGTTCCGGGGGAAAATGAAAGCCCCATGGTACCAAATACACCGCACCCGGGATGcac            >  1:557752/1‑226 (MQ=37)
                                                                                                                                                                                                            cTGACCCCAGCGGAGCGTGAAAAGAAAGCCGGCAAGACAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCAGTTTCGCACGGCGCGCACGGTGGTCAAAACGGCATTGGTCTGGGTATGTTGGGCTTTGTTGGGGTCACGCCAGCGGGCTTCGTGGTGAACATGAACGCCACTGGGCAGGGAATCacacggcccctggagtacaca          >  1:99864/1‑212 (MQ=255)
                                                                                                                                                                                                            cTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGGAAAAGCCGCCGTTCTGGa                                                                                                                                                                                           >  1:328001/1‑55 (MQ=255)
                                                                                                                                                                                                               aCCCCAGCGGAGCGTGGAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACCGCGCGAACAGTGGTCAGAAAGGCCTTGGGCTGGTTATGTTGGTATTGATTTGCGTCCCGCCAGCAGGCTTCGTGGGGAACATGAATGCACCTGGCTTCGAAAACAACCGGACACCGGATGGCCTcaaa       >  1:249235/1‑227 (MQ=255)
                                                                                                                                                                                                                    aGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCCCGCCAGCAGGCTTCGTGGGGAACATGAATGCCACTGGCTACGCAATCACCCCGACCCGTGATGCCATCAACAACGt  >  2:606356/1‑228 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
ATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGT  >  NC_000913/3638013‑3638449

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