Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,551,917 A→C 18.3% V186G (GTG→GGG)  malP ← maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,551,9170AC18.3% 106.9 / 2.1 59V186G (GTG→GGG) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (35/10);  new base C (4/8);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 5.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCGCCGGG  >  NC_000913/3551694‑3552150
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
gTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGCCAACGGGGGGATCCCACGCTTGACCGGGAATGGTAAACTCCGGCTCCCAGGGTCCGCCTTTCGTACCTTTACCGCCAATc                                                                                                                                                                                                                             >  1:382985/1‑238 (MQ=255)
gTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGCCAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCCCCCGGACTCCGTCTTTCTTCACTTTACCCCCAATc                                                                                                                                                                                                                             >  2:774321/1‑238 (MQ=255)
gTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACAGTAGCCGACAACGGGGGGATACCACGCTTGGCCGGTAATGGTAAACTCCGCTTCCCAACGGCCGCCTTTTGCACCTTTACCGCCatac                                                                                                                                                                                                                             >  1:1250998/1‑235 (MQ=255)
 tCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGCCCGTCTTTCGTCACTTTACCGCCAATc                                                                                                                                                                                                                             >  1:842872/1‑237 (MQ=255)
      tGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCACCGTCCGTCTTTCGTCACTTTACCGCCATTCCCTAcc                                                                                                                                                                                                                       >  1:1190659/1‑238 (MQ=255)
        gATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCAAAGACGCAGCGGCTGCGCAACGCCATTACGATAGCCGACAACGGGTAGATCCCACGCTTGACCGGCAATGGTAAACTCCGGCTCCCAGCGTCCGCCTTTCGTCACTTTACTCCCTATCCCTACCGGc                                                                                                                                                                                                                    >  1:1141934/1‑239 (MQ=255)
             agaACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTCCCGCCAATCCCTACCTGCACAt                                                                                                                                                                                                                >  1:316620/1‑238 (MQ=255)
             agaACTTTGGTCAGATTATCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGAAAGCCGACAACGGGCAACTCCCACGCTTGACCGCTAATGGGAAACTCCGGCCCCCAGCGACCGGCTTTCGACACTTTCCACACAACCCCTACCCGcccacc                                                                                                                                                                                                               >  2:1303238/1‑237 (MQ=255)
                  tttGGTGAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCCCGCCATTACAATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTATTGGAGAACGCCGGCTCCCAGCGGCCGTCTTTCGGCACTTTACCGCCAACCCCTACCTGcacacacacgg                                                                                                                                                                                                          >  2:420057/1‑232 (MQ=255)
                  tttGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTCCCGCCAATCCCTACCTGCACCTCCAgcg                                                                                                                                                                                                          >  1:1031491/1‑237 (MQ=255)
                  tttGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTCCACCTCCAGt                                                                                                                                                                                                           >  1:243781/1‑238 (MQ=255)
                  tttGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAAACCCGCCTCCCAGCGCCCGTCTTTCGGCCCTTTACCCCCAAAacctacctcccaaccaaag                                                                                                                                                                                                           >  2:132968/1‑227 (MQ=255)
                  tttGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGGAATGGGAAACTCCGGCTCCCACCGTCCGCCTTTCTCCCCTTTACCCCCAAACCCTACCTGCACACCCAgg                                                                                                                                                                                                           >  2:914374/1‑237 (MQ=255)
                            ttttCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTtgtg                                                                                                                                                                                                >  1:504380/1‑239 (MQ=255)
                                       tgtGCCGTGTTGTCCGGCACGAAAAAACCCCCCGTGTTAAAATTGGGTCATACCAAGCGGAGGCGGGTGCTCCCCGGCCCACACCCGGAGGGGCTGCGCCACGCCATTACGATAGCCGACAGCGGGAAGTTCCCACGTTTGACCGTTAATGGTAAACCCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACg                                                                                                                                                                                     <  2:946063/237‑1 (MQ=255)
                                                ctgTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGCCACTTTACCGCCAATCCCTACCTGCCCATCCAGGGCTTCGTTGTGGCGGAAACACGGGGAGtt                                                                                                                                                                              >  1:254809/1‑236 (MQ=255)
                                                   ttCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACCACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTCCACATCCAGGGCTTCGTTGTGTCGGAACCACgggagcgtaacg                                                                                                                                                                          >  2:84545/1‑229 (MQ=255)
                                                        cACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCCCATCCAGTGCTTCGTTTTGGCGGAACCACCGGTAGTTACTGCGATGc                                                                                                                                                                    >  1:87422/1‑239 (MQ=255)
                                                                       cgtcgtTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCgcgcct                                                                                                                                                     >  1:339001/1‑234 (MQ=255)
                                                                          cgtTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTa                                                                                                                                                                             >  1:878642/1‑212 (MQ=255)
                                                                           gtTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGGGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCa                                                                                                                                                  >  1:366131/1‑238 (MQ=255)
                                                                                                cGGATGCGCGTGCGTCGCCCGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACaaaaa                                                                                                                            >  1:599125/1‑238 (MQ=255)
                                                                                                              gccccTGCCACAGACGAAGCGGGTGCGCATCCCCATTACGAAAGCCGACAAGGGGAAACTCCCACGCTTGCCCGTTAATGGTAACCTCCGGCTCCCGGCTTCCGTCTTTCGTCGCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACaa                                                                                                               <  1:84545/235‑1 (MQ=255)
                                                                                                                        accgccgcgggggCTGGCCCACGCCATTCCGATAGCCGACAACGGGGAGACCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGa                                                                                                     <  1:914374/228‑1 (MQ=255)
                                                                                                                         aaccgcgaGCGGCCGCGCCCCCGCCATTACAATAGCCGAAACGGGGGAGACCCCCCTGTTGACCGGGAAGGGTAAACCCCGGCTCACAGCTCCCGTTTTTCGTCCCTTCACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGa                                                                                                     <  2:68258/232‑1 (MQ=255)
                                                                                                                              gcgccgcgcgcgcCCCGCCTTTCCGATACCCGACACGGGGGATCTCCCCCGCTTTACCGGTAATGGTAAACCCCGCCCCCCAGCGTCTGTCTTCCGTCATTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTc                                                                                               <  2:339001/230‑1 (MQ=255)
                                                                                                                               gaGCGGCTGCCCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTGCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCa                                                                                              <  1:804773/238‑1 (MQ=255)
                                                                                                                                aGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAg                                                                                             >  1:68258/1‑239 (MQ=255)
                                                                                                                                aGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCACGCGCTTCAACCCGTTTGCACTACAAAAAAAACCGGGGGAACAAAACATATTGaaaatt                                                                                                >  2:1046998/1‑231 (MQ=255)
                                                                                                                                        cgccacgccaTTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGACCCCCGGGTAATTACTGCGATGCCAGTCATCCGGGGCTTCAACCTTTTTGCCATCGACAAAAGACTCGCGGAACAAACCCTAATTAAAGTTCAAAcccgaa                                                                                      >  2:993866/1‑234 (MQ=255)
                                                                                                                                          ccccccccaTTCCGAGAGCCGAAAACGGGGACTTCCCCGGCTTGCCCGTTAATGGAAAACTCCGGCCCCCGGGGCCCTTCTTCCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTGCTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAAc                                                                                     <  1:1106233/232‑1 (MQ=255)
                                                                                                                                               gccaTTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCAt                                                                                                                                     <  1:1057032/185‑1 (MQ=255)
                                                                                                                                               gccaTTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCAt                                                                                                                                     >  2:1057032/1‑185 (MQ=255)
                                                                                                                                                   gtACGATAGCCGACACGGGGGAAACCCCCCGCTTGCCCGGTAATGGTAAACTCCGCCTCCCAGCCTCCTTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGa                                                                           <  1:774321/238‑1 (MQ=255)
                                                                                                                                                      cagttgCCAACACCGGGGGGTCCCCACGGTTGTCCAGTAATGGTAACCCCCGTCCCCCGGCGTCCGTCTTCCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTTATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGATTg                                                                        <  2:87422/234‑1 (MQ=255)
                                                                                                                                                         tAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATc                                                                                                                                                                                                               <  1:181825/101‑1 (MQ=255)
                                                                                                                                                         tAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATc                                                                                                                                                                                                               >  2:181825/1‑101 (MQ=255)
                                                                                                                                                           gCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCCCAGACTGACCGa                                                                   >  1:537482/1‑239 (MQ=255)
                                                                                                                                                           gCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCAAATTGTTAATTCAAACCGTAACCCGTCGCAGaccgaccg                                                                    >  1:841655/1‑238 (MQ=255)
                                                                                                                                                               gaagcGGGGAATTCCCTCCCTTGAGCGGTAATGGTAACCCCCGGCCCCCAGCGTCCGTTTTTCGTCACTTTACCGCCAATCCCTACCTGCTCATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGt                                                               <  1:532278/235‑1 (MQ=255)
                                                                                                                                                                      ggAGATCCCCCGTTTCCCCGGTAAGGTTAAACCCCGCCCCCCCGTCTCCGTTTTCCGTCCTTTCACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCAt                                                         <  2:949179/238‑1 (MQ=255)
                                                                                                                                                                      ggAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTTATAGTTCAGACCGTAACCCCTCGCAAACTGACCGCCAGTTTCCAtt                                                        >  1:911916/1‑239 (MQ=255)
                                                                                                                                                                               aCGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCAAATTGATAGTTCAGACCGTAACCCGCCGCAGACGACCCGCCAATTGCCCtttggacgcgg                                               >  1:361331/1‑230 (MQ=255)
                                                                                                                                                                                cGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCAAATTGATAGTTCAGAACGGAACCCGTCGCAGACTGAACGACAGGTGCCATTGAGTCGAGGa                                              >  1:554220/1‑239 (MQ=255)
                                                                                                                                                                                  cTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGTCGGAACAAACCATATTGATAGTTCATACCGTCACCCGCCGCCAACTGACCGCCAATTGCCCTTTGGTCCAGGAAg                                            >  1:488807/1‑239 (MQ=255)
                                                                                                                                                                                     gACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGa                                                                                                                                                                       <  1:1280008/113‑1 (MQ=255)
                                                                                                                                                                                     gACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGa                                                                                                                                                                       >  2:1280008/1‑113 (MQ=255)
                                                                                                                                                                                     gACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCCTTGAGTCGAGGAAGCa                                          >  1:66393/1‑238 (MQ=255)
                                                                                                                                                                                              tGGTAAACCCCGGCTCCCGGCTCACGTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGa                                <  1:757526/239‑1 (MQ=255)
                                                                                                                                                                                                gtaaccccgggCTCCCGGTGTCCGTCTTACGTCCCTGTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACg                              <  1:1303238/230‑1 (MQ=255)
                                                                                                                                                                                                        cgccTCCCTGCTCCCGTTTTTCGTGTCTTTCCCCCCAACCCCTACCCGCTCATCCAGGGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGac                      <  2:504380/236‑1 (MQ=255)
                                                                                                                                                                                                        cGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTAcc                                                                                                                                                                                                                       >  1:249857/1‑46 (MQ=255)
                                                                                                                                                                                                        cGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAAACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACGGACCGACCAATGCCATTTAGTCCAGGAAGCACGGCGCCAGAACACCCAGac                      >  1:339051/1‑239 (MQ=255)
                                                                                                                                                                                                        cGGCCCCCAGGCCCCGTCTTTCCTCCCTTGTCCGCCAATCCCTTCCCGCCCATCCCGTGCTTCGTTGTTGCCGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGac                      <  2:366131/239‑1 (MQ=255)
                                                                                                                                                                                                        cGGATCCCAGGGTCCGTCTTCCGCCCTTTTACCGCCAATCCCTGCCTGCACAGCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGac                      <  1:508475/239‑1 (MQ=255)
                                                                                                                                                                                                         ggCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTTACCGACAGTTGCCATTGAGTCACGGAAGCAAGCCGCCAGACGTCCCAGac                      >  1:1078440/1‑238 (MQ=255)
                                                                                                                                                                                                                        tcTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCCTATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCCAACGTCCCAAACCACCGTTACCCAgcg       >  1:637231/1‑238 (MQ=255)
                                                                                                                                                                                                                         tttccGTCCCTTTACCGCCAATCCCTGCCGGTACAACCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAgcgc      <  1:1275855/234‑1 (MQ=255)
                                                                                                                                                                                                                            tCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGCCCGACAGTTGCCATTTAGTCGAGGAAGAACGCCGCCAGACGTCCCACACCACCGTTACCCCGCGCCggg  >  1:1360183/1‑239 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
GTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCGCCGGG  >  NC_000913/3551694‑3552150

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