Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 417,391 G→C 100% A84P (GCC→CCC)  phoB → response regulator in two‑component regulatory system with PhoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,3910GC56.4% 11.6 / 66.7 39A84P (GCC→CCC) phoBresponse regulator in two‑component regulatory system with PhoR
Reads supporting (aligned to +/- strand):  ref base G (12/5);  new base C (10/12);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-02
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

CGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAG  >  NC_000913/417112‑417665
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                  
cGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAgggg                                                                                                                                                                                                                                                                           >  2:293279/1‑289 (MQ=255)
                   caGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTgc                                                                                                                                                                                                                                                         >  2:40765/1‑288 (MQ=255)
                       gCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTgcg                                                                                                                                                                                                                                                        <  2:112336/285‑1 (MQ=255)
                                       cGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGcgcg                                                                                                                                                                                                                                     <  2:294207/288‑1 (MQ=255)
                                          aTTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGAt                                                                                                                                                                                                                                  <  1:38775/288‑1 (MQ=255)
                                           tgcTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGa                                                                                                                                                                                                                                <  1:126966/287‑1 (MQ=255)
                                                                          cgAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGa                                                                                                                                                                                                  >  1:156207/1‑288 (MQ=255)
                                                                                  tCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGgc                                                                                                                                                                                         <  1:40765/289‑1 (MQ=255)
                                                                                   cTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAATCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGgcg                                                                                                                                                                                        <  1:293279/289‑1 (MQ=255)
                                                                                              tCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCGGGATTTAATTCTCCGCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAg                                                                                                                                                                              <  1:128732/288‑1 (MQ=255)
                                                                                                gAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCgg                                                                                                                                                                           >  2:55928/1‑289 (MQ=255)
                                                                                                            tttCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTa                                                                                                                                                               >  1:311298/1‑289 (MQ=255)
                                                                                                                aGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTgcgc                                                                                                                                                                                                                                                       >  1:6698/1‑197 (MQ=255)
                                                                                                                aGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTgcgc                                                                                                                                                                                                                                                       <  2:6698/197‑1 (MQ=255)
                                                                                                                          aGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCg                                                                                                                                                 <  2:162200/289‑1 (MQ=255)
                                                                                                                                       gACAGTGCTGTTAATCAACTGACTGAACCCTGTCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTa                                                                                                                                                                         <  2:1692/252‑1 (MQ=255)
                                                                                                                                         cAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGAtt                                                                                                                                  >  2:41751/1‑289 (MQ=255)
                                                                                                                                           gTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGAGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGCGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGCGGTGATTga                                                                                                                                >  1:284891/1‑289 (MQ=255)
                                                                                                                                            tGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAACGCGGTAATGCGCCGTATGTCGCCAATGGCTGTGGGAGAGGTGATTgag                                                                                                                               >  1:167433/1‑289 (MQ=255)
                                                                                                                                            tGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCGGTATTTCGCCAATGGCGGTGGACGAGGTGATTgag                                                                                                                               >  1:60933/1‑289 (MQ=255)
                                                                                                                                            tGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGc                                                                                                                                                  <  1:125428/270‑1 (MQ=255)
                                                                                                                                            tGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGc                                                                                                                                                  >  2:125428/1‑270 (MQ=255)
                                                                                                                                            tGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGACGAGGTGATTgag                                                                                                                               >  1:284272/1‑289 (MQ=255)
                                                                                                                                            tGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTgag                                                                                                                               >  2:68116/1‑289 (MQ=255)
                                                                                                                                                  gAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCa                                                                                                                          <  1:68116/288‑1 (MQ=255)
                                                                                                                                                                  cccTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGAcc                                                                                                         <  1:41751/289‑1 (MQ=255)
                                                                                                                                                                                      cctcGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATg                                                                                     >  2:180107/1‑289 (MQ=255)
                                                                                                                                                                                       ctcGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGGAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCCATGGCGGTGGAAGAGGTGCTTTAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATg                                                                                     >  1:154176/1‑288 (MQ=255)
                                                                                                                                                                                                            ggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGgag                                                                <  2:185669/288‑1 (MQ=255)
                                                                                                                                                                                                                           ttCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGTCGTATTTCGCCAATGGCGGTGGTAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATATCACCGAGTGATGGGCGGCGAAGAGCCGCTGGGGCTGGGTCCGACAGAAt                                                >  1:133816/1‑289 (MQ=255)
                                                                                                                                                                                                                              atcaAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTa                                             <  2:156207/289‑1 (MQ=255)
                                                                                                                                                                                                                                caAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTaaa                                           >  2:269704/1‑289 (MQ=255)
                                                                                                                                                                                                                                     acCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAAct                                         >  1:50918/1‑286 (MQ=255)
                                                                                                                                                                                                                                        tCAAGGCCGAGCCGAGGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCa                                   <  2:265743/289‑1 (MQ=255)
                                                                                                                                                                                                                                               cgAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGTGATGGGGCCGACAGAATTTAAACTGCTGCACtttttt                            >  1:5916/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      aTGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGCGCCGCTGGAGATGGGGCCGACAGAATGTAAACTGCTGCACTTTTTTCTGACg                      >  2:262215/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                         tggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAg  >  1:223556/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         tggtgATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGc                                                         <  2:160080/234‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tggtgATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGc                                                         >  1:160080/1‑234 (MQ=255)
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                  
CGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGCCGGTCGAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAG  >  NC_000913/417112‑417665

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