Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,321,743 T→A 7.3% G567G (GGT→GGA atoS → sensory histidine kinase in two‑component regulatory system with AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,321,7430TA 7.3% 187.1 / 2.0 78G567G (GGT→GGAatoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (4/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

GATTAACGCCGATCGTGAACTGCTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGAAACACATTGTGCG  >  NC_000913/2321518‑2321967
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
gATTAACGCCGATCGTGAACTGCTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGCCAACGGCTGTGGCATTGTTCCCTCGCCGCAAAAAAAAATCTTCGCACCCTTTTTTCCCCCCAAAACCCCAAGAAACGGGCTTGGTcc                                                                                                                                                                                                                                >  2:425760/1‑227 (MQ=255)
          gATCGTGAACTGCTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGACTGGCGTTAAGt                                                                                                                                                                                                                      >  1:168031/1‑228 (MQ=255)
                 aaCAGCTCAAACAGGTACTAATGAATATCCTGATCAATGCCGTGCAGGCTAACAGCGCACGAGGGAAAATTCGCATTCAAACCTTGCAATACAGCTACTCCAAACAGGCCATTTCGATAGAGGACAACGGATGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTAAACGca                                                                                                                                                                                                               <  2:94211/228‑1 (MQ=255)
                      cTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTCCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGTCGTTAAGTCACCTCACCAtt                                                                                                                                                                                                          >  2:223597/1‑228 (MQ=255)
                      cTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCACACGCTGCAAAAAAAAAATTTCGATCCCTTTTTCACCCCCAAAGCCTCAGGAACCAGACTTTGACAGTCCTTAAGTCAACCcacaatg                                                                                                                                                                                                          >  2:274286/1‑223 (MQ=255)
                         aaacgtgTACTACGGAATTTCCCGAACAACGCCGTCCAGGCTATCCGCGCACGAGGGAAAATTCGCATTCAAACCGGGCAATACAGCGACTCACAACAGGCCATTTAGATAGAGGACAACGGCTGTTGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAAt                                                                                                                                                                                                       <  2:168031/222‑1 (MQ=255)
                         aaaCAGGTCCTACTGAATATCCTGAACAATGCCGTCCAGGCTAACAGCGCACGAGGGAAAATTCTGATTCAAACCTGGCAATACAACGACTCACAACAGGCCATTTCGATAAAGGACAACTGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAAt                                                                                                                                                                                                       <  1:34257/228‑1 (MQ=255)
                         aaaCAGGGACTACCGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAAt                                                                                                                                                                                                       <  1:223597/228‑1 (MQ=255)
                          agagacagccactGAATACCCTGATCAATGCCGGCCAGGCTAACAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAAt                                                                                                                                                                                                       <  1:676821/217‑1 (MQ=255)
                          aaCAGGTCCTACTGAATATCCTGATCAATGCCGTACAGGCTATCAGCGCACGAGGGAAAGTTCGCATTCAAACCTTGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAAt                                                                                                                                                                                                       <  1:251440/227‑1 (MQ=255)
                              gtgacttatGAATAGTCTGACCAATGCCGTCCCGGCTATCAGCGCACGAGGGAAAATTCGAATTCAAACCGGGCAATACAGCGACTCACAACAGGCCAATTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCa                                                                                                                                                                                                  <  2:121506/220‑1 (MQ=255)
                               gTACTACGAAAAACCCTGATCAAGGCCGTCCAGGCTAACAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGCATCACCACAGGCCATTTCGATAGAGGACAACGGCTGTTGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCAt                                                                                                                                                                                                 <  2:656690/228‑1 (MQ=255)
                                tactacTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCAt                                                                                                                                                                                                 >  1:475714/1‑227 (MQ=255)
                                tactacTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAAAGGCTGTGGCATTGATCTCTCGCTCCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAATCCTCAGGAACCGGGCATGGTCTCGCGTTAATCCAACGCATCACTAATGCCCAt                                                                                                                                                                                                 >  1:633795/1‑227 (MQ=255)
                                tactacTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCACACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCAt                                                                                                                                                                                                 >  1:694654/1‑227 (MQ=255)
                                          tatCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGCCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGGCCGACGTTAACTAAACGCCTCATTAATGCCCACCCGGGGGat                                                                                                                                                                                        >  1:597983/1‑226 (MQ=255)
                                             ccTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCCAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGAAACATTAATGACCATCAGGGGGGTATTcg                                                                                                                                                                                   >  2:473093/1‑228 (MQ=255)
                                             ccTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTAACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCTTCATTAATGCCCATCAGGGTGATATTcg                                                                                                                                                                                   >  1:546607/1‑228 (MQ=255)
                                             ccTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGATGTGGCAATGATCTCTCACTCCAAAAAAAGATCTTCGATCCCTTTTTCACACCCAAAGCCTCACGAACCGGGATTGGGCTGGCGTTTCTTAAAACCAAAATTAACGCCCACCAgggcgggtcgct                                                                                                                                                                                   >  2:676821/1‑219 (MQ=255)
                                                     aTGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGCGGACATTCGCGTCGcc                                                                                                                                                                            >  1:224431/1‑227 (MQ=255)
                                                       gCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGt                                                                                                                                                                                                                                  <  1:204559/171‑1 (MQ=255)
                                                           tCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGCTAGAGGACACCGGCTGTAGCATTGATCTCTCTCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAATTCAACCCATCATTCATGCCCATCAGCGTGTTATTCGCGTGCCCAttttac                                                                                                                                                                     >  2:697942/1‑226 (MQ=255)
                                                           tCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGGGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTATGCCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGGAATTCGCGACGCCAGTTTg                                                                                                                                                                      >  1:683973/1‑227 (MQ=255)
                                                                gatATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGc                                                                                                                                                                <  1:605774/226‑1 (MQ=255)
                                                                    tCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAACGATCTTCGAACCCTTTTTCACCACCAAAGCCCCAGGAACCGGGATTGGTCTGCCGTTAGGTCACCGCACCATTAATCCCCACCACGGTGATTTTTCGGTGGCCGTTTTCCCGGGCtca                                                                                                                                                             >  2:583346/1‑225 (MQ=255)
                                                                     cAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGCTCTTCGCTCCCTTTTTCACCACCCAAGCCTCAGGAACCAGGCTTGGTCTGGCGTTCAGCCAACGCACCATTAATGCCCATACGGGCGCAATTCGCGCCGCCAGTTTTCCGGGCTCCg                                                                                                                                                            >  1:203136/1‑227 (MQ=255)
                                                                        cgcACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGGGATATTCTGGTCGCCGGTTTGTCGGGCTACGgc                                                                                                                                                          >  1:677221/1‑226 (MQ=255)
                                                                                    aaTTTGCATTCAAAACTGGCATAACAGCGACTCACAACAGGCCATTTCGATAGGGGACAACGGCTGTGGCATTGATCTCCCGCTGCAAAAAAAGATCTTCGAACCCTTTTTCACCCCCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCAc                                                                                                                                               <  2:26391/225‑1 (MQ=255)
                                                                                    aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAAAAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCAc                                                                                                                                               <  1:68234/225‑1 (MQ=255)
                                                                                                   cgggCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATc                                                                                                                               <  1:10739/224‑1 (MQ=255)
                                                                                                     ggCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGACTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCGGTTTGCCGGGCTACGGGCCAACCTTCACCCTTAgtttccagcccaacc                                                                                                                           >  2:655672/1‑217 (MQ=255)
                                                                                                             aGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGa                                                                                                                   >  1:114093/1‑228 (MQ=255)
                                                                                                               cGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAAc                                                                                                                                                                                                                  >  1:399781/1‑131 (MQ=255)
                                                                                                               cGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAAACGGGCTTGGTCTGCCGTTAAGTCAACGCATCATCAATGCCAACCGGGGGGATATTGCGGTCGCCAGTTTCCCGGGCTAAGCCCCCCCCTTACAGCTTTTTTTACCCACCGAccccaccgggaa                                                                                                                 >  2:272422/1‑219 (MQ=255)
                                                                                                               cGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAAATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGGCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGGGATATTCACGTACCCAGTTTGCCGGTCTACGGCGCCACCTTCACGCTTATTTTACCGTCCACCCCCCCGGGaac                                                                                                                 >  2:599386/1‑227 (MQ=255)
                                                                                                                gACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACg                                                                                                                                                                                                                 <  2:399781/131‑1 (MQ=255)
                                                                                                                caCTCACAACAGGCCAATTCAATAGAGGACAACGGCTGTGGCATTGATCTCTCTCTGCAAAAAAAGATCTGCGATCCCTTTTACACCACCAAAGCCTCAGGAACCGGGCTTGGACTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAAt                                                                                                                <  1:655672/227‑1 (MQ=255)
                                                                                                                 aCTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGCTATTCGCGTCGCCCGTTTGCCGGGCTACGGCGCAACCTTCCCCCTTATTTTACCCATCAACCCGCCGGGaaac                                                                                                                >  2:467213/1‑226 (MQ=255)
                                                                                                                     gcaacaGGGCATTTCGATAGAGGCACACGGCTGGGGCATTGATCTCACGCTGCAAAAAAAGATCGTCGAGCCCTTTTTCACCACCAATGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTGGCCGGGCTACGGCGCAACCGTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGa                                                                                                            <  1:408002/226‑1 (MQ=255)
                                                                                                                          agaccATTTCGAAGAGGGCCAACGGCTTTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGGTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTAt                                                                                                      <  2:133583/225‑1 (MQ=255)
                                                                                                                           ggCCATTTCGATAGAGGACCACGGCTGTGGCATTGATCTCCCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTAt                                                                                                      <  1:467213/227‑1 (MQ=255)
                                                                                                                            gCCATTTCTATAGTGGCCAACGGCTTGGGCATTTAATCCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCCACACCAAAGCCTCAGGACCCGGGCTTGGTCTGGCGTTAAGTCAACGACTCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGa                                                                                                    <  2:570450/228‑1 (MQ=255)
                                                                                                                            gCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGa                                                                                                    <  1:274286/228‑1 (MQ=255)
                                                                                                                            gCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATACCTTTTTCACCACCAAAACCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCCGGCTACGGCTCAACCTTCACGCTTATTTTACCGATCAACCTGCAGGGCAATCAGACTGTATg                                                                                                     >  1:697198/1‑227 (MQ=255)
                                                                                                                            gCCATTTCGATAGAGGACAACGGCCGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCCCAGGGAAATCAGACTGTa                                                                                                       >  1:133583/1‑225 (MQ=255)
                                                                                                                                 ttCGATAGCGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCTGGCTACGGCGCAACCTTCACGCTTCTTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGc                                                                                                >  1:518197/1‑227 (MQ=255)
                                                                                                                                        gagGACAACGGCTGGGGCATTGCGCTCTCCGTGCAAAAAAAGATCTTCGATCACTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAAt                                                                                         <  1:630142/227‑1 (MQ=255)
                                                                                                                                        gagGAACACGGCTGTGGCATTGATCTCCCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAAt                                                                                         <  2:114093/227‑1 (MQ=255)
                                                                                                                                         agGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGAGCCCTTTTTCACCACCAAAGCCGCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCg                                                                                       <  1:543963/228‑1 (MQ=255)
                                                                                                                                            gACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGCCTGTATGACTGCTATTAATCGCa                                                                                     >  1:600714/1‑228 (MQ=255)
                                                                                                                                              aaaCGGCTGTTGCATTTATCTCTCCCTGCAAAAAAATATTTGTGGTCCCTTTGTCCCCCCCAAAGCCTGAGGAACCGGGCTTGGTCTTGCTTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATc                                                                                   <  2:224431/227‑1 (MQ=255)
                                                                                                                                                 cGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGCAACCGGGCTTGGCCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTCTTTGACTTCTTTTAATCGCATCCtt                                                                                >  1:402619/1‑228 (MQ=255)
                                                                                                                                                   gCTGTGGCATTGATCTCTCCCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTAt                                                                              <  1:425760/228‑1 (MQ=255)
                                                                                                                                                      gtgGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATtgt                                                                           >  1:418525/1‑228 (MQ=255)
                                                                                                                                                      gtgGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGATTTACTGCTATTAATCGCATCCTTTTtgg                                                                           >  1:564950/1‑227 (MQ=255)
                                                                                                                                                      gtgGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGGGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTCTTTTACCGTTCAACCCGCAGGGTAATCAGACTGTATGACTGCTATTAATCGCATCCTTAtttt                                                                           >  1:242282/1‑226 (MQ=255)
                                                                                                                                                       tgGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGCCTGCTATTATTCGCATCCTTATtggg                                                                          >  1:655571/1‑226 (MQ=255)
                                                                                                                                                            ttGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCa                                                                                                                                                                                                  >  1:369340/1‑103 (MQ=255)
                                                                                                                                                             tGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCAt                                                                                                                                                                                                 <  2:369340/103‑1 (MQ=255)
                                                                                                                                                                   cccgttGCAAAAAAAGATCTTCGATCCCTTTTTAACCCACCAAGCCTCAGGAACCGGGCTTGGTCTGGCCGTAAGTCAACGAATCATTAATGCCCATCAGTGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGAt                                                              <  2:694654/223‑1 (MQ=255)
                                                                                                                                                                        tGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGt                                                         <  1:253144/228‑1 (MQ=255)
                                                                                                                                                                         gCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAACGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGTCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATTACTGCTATTAAACCCATCCTTATTGGGGATGATGAAGATAAt                                                           >  2:564570/1‑225 (MQ=255)
                                                                                                                                                                         gCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATCATGAAGATAATGtt                                                        >  1:470182/1‑228 (MQ=255)
                                                                                                                                                                                         gATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGc                                         >  1:387617/1‑227 (MQ=255)
                                                                                                                                                                                          gtCCCTTTTTCCCCCCCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCAc                                       <  1:497713/227‑1 (MQ=255)
                                                                                                                                                                                          aTCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCAc                                       <  1:564570/228‑1 (MQ=255)
                                                                                                                                                                                          aTCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCAc                                       <  1:697942/228‑1 (MQ=255)
                                                                                                                                                                                          aTCCCTTTTCCACCACCAAGGCCTCAGGAACCGGGCTTGGGCTGGCCTTAAGTCAACGTATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCAc                                       <  1:110791/228‑1 (MQ=255)
                                                                                                                                                                                            cccTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCGGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGTTGAAAATAATGTTCGCCGTATTCTGAGCCCCg                                     >  2:482133/1‑228 (MQ=255)
                                                                                                                                                                                                ttttCACCACCAAAGCCTCAGGAACCGGGCTTGGACTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGCTCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGAGGATGAAGCTAAAGTTCGCCGAATGATGAGCACCCCttt                                 >  1:502167/1‑228 (MQ=255)
                                                                                                                                                                                                  ttCACCACCAAAGCCTCAGGAACCTGGCTTGGTCTGGCTTTAAGTCCACGCATCATTAATGCCCATCAGGGTTATTTTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTg                               <  1:693239/228‑1 (MQ=255)
                                                                                                                                                                                                  ttCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCCTTAAGGCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTg                               <  1:677122/228‑1 (MQ=255)
                                                                                                                                                                                                        accaAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCCTCCTTTTTGTGGATGATGAAGATAATGTTACCCGAATGCGTAGCCCCCCTTTTGccccac                         >  2:601052/1‑223 (MQ=255)
                                                                                                                                                                                                          caATGCCTCCGCATCCGGGCTTGGTCTGACGTTAATGCATCGCATCATGACTGCCACACAGGTTGATATTCCCGTCGCCAGGTTGCCGGGCTACAGCGCAACCTTCACGCTTCTTTTACCCATCAACCCGCAGGGAAATCAGACTGTATGACTGGTATTAATCGCATCCTTAATGTGGATGATGAAGATAATTGACGCCGTATTCTGAGCACCTCTTTGGCACTACaa                       <  2:103016/228‑1 (MQ=255)
                                                                                                                                                                                                                  ccggAACCGGGTTTGGGCTGGGGTTAAGGCAACGCCTTATTAAGGCCCATCAGGGTGATATTCGCGGCGCCAGTTTGCCGGGCTACGGGGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGa               <  2:655571/226‑1 (MQ=255)
                                                                                                                                                                                                                   aGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGa               >  1:75723/1‑227 (MQ=255)
                                                                                                                                                                                                                   aGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGa               <  2:75723/227‑1 (MQ=255)
                                                                                                                                                                                                                             cTTGGACTTGAGTTAAGACACCGCATAATTAATGCCCACCAGGGTGATAGTCGCGTGCCCAGTTTGCCGGGGTACGGCGCAACCTTCACGCTTATTTTACCGATTAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGAAACACATtgtg    <  2:470182/228‑1 (MQ=255)
                                                                                                                                                                                                                               ttgTCTGGCGTTAAGTCATCGCTTCATTAATGCCCATCCGGGTGATTTTCGCGTCGCCAGTTTGCCGTGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGAAACACATTGTGCg  <  2:683973/226‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
GATTAACGCCGATCGTGAACTGCTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGGCATTGATCTCTCGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGAAACACATTGTGCG  >  NC_000913/2321518‑2321967

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