Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,551,915 A→C 15.7% 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,9150AC15.7% 136.1 / 2.9 73V186G (GTG→GGG) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (39/19);  new base C (4/9);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCAGTATGGTTGTCGTTTGGATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCGCCGGG  >  NC_000913/3551681‑3552148
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
ggttttttttttttgTTTGGATAGA‑AACTTTGGCGCAGTTTTCCCCATTTGACCCCCTCTTGTCCGGCCCGCAAAAAATCCCCGTCGTTAAATTTAGTCAGATCAACGGGCGGCGCGTGGGG‑CGCCTG‑CCACAGCCCCGGCGCCTGCC‑CCCGCCCATTACGCTAGCCGCCAACGGGGAGAC‑CCCCCGTTTGACCGGTAATGGTAACCTCCGGCTCCCAGCGTCCGTCTTTCGTCACttt                                                                                                                                                                                                                                       <  2:35991/225‑1 (MQ=255)
           gTCGTTTGGATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGCCAACGGGGGGAT‑CCCACGCTTGACCGGGAATGGTAAACTCCGGCTCCCAGGGTCCGCCTTTCGTACCTTTACCGCCAATc                                                                                                                                                                                                                             >  1:382985/1‑238 (MQ=255)
           gTCGTTTGGATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGCCAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCCCCCGGACTCCGTCTTTCTTCACTTTACCCCCAATc                                                                                                                                                                                                                             >  2:774321/1‑238 (MQ=255)
           gTCGTTTGGATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACAGTAGCCGACAACGGGGGGAT‑ACCACGCTTGGCCGGTAATGGTAAACTCCGCTTCCCAACGGCCGCCTTTTGCACCTTTACCGCCatac                                                                                                                                                                                                                             >  1:1250998/1‑235 (MQ=255)
            tCGTTTGGATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGCCCGTCTTTCGTCACTTTACCGCCAATc                                                                                                                                                                                                                             >  1:842872/1‑237 (MQ=255)
                 tGGATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCACCGTCCGTCTTTCGTCACTTTACCGCCATTCCCTAcc                                                                                                                                                                                                                       >  1:1190659/1‑238 (MQ=255)
                   gATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCAAAGACGCAGCGGCTGCG‑CAACGCCATTACGATAGCCGACAACGGGTAGAT‑CCCACGCTTGACCGGCAATGGTAAACTCCGGCTCCCAGCGTCCGCCTTTCGTCACTTTACTCCCTATCCCTACCGGc                                                                                                                                                                                                                    >  1:1141934/1‑239 (MQ=255)
                        agaACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTCCCGCCAATCCCTACCTGCACAt                                                                                                                                                                                                                >  1:316620/1‑238 (MQ=255)
                        agaACTTTGG‑TCAGATTATCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGAAAGCCGACAACGGGCAACT‑CCCACGCTTGACCGCTAATGGGAAACTCCGGCCCCCAGCGACCGGCTTTCGACACTTTCCACACAACCCCTACCCGcccacc                                                                                                                                                                                                               >  2:1303238/1‑237 (MQ=255)
                             tttGG‑TGAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCCCGCCATTACAATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTATTGGAGAACGCCGGCTCCCAGCGGCCGTCTTTCGGCACTTTACCGCCAACCCCTACCTGcacacacacgg                                                                                                                                                                                                          >  2:420057/1‑232 (MQ=255)
                             tttGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTCCCGCCAATCCCTACCTGCACCTCCAgcg                                                                                                                                                                                                          >  1:1031491/1‑237 (MQ=255)
                             tttGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTCCACCTCCAGt                                                                                                                                                                                                           >  1:243781/1‑238 (MQ=255)
                             tttGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAAACCCGCCTCCCAGCGCCCGTCTTTCGGCCCTTTACCCCCAAAacctacctcccaaccaaag                                                                                                                                                                                                           >  2:132968/1‑227 (MQ=255)
                             tttGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGGAATGGGAAACTCCGGCTCCCACCGTCCGCCTTTCTCCCCTTTACCCCCAAACCCTACCTGCACACCCAgg                                                                                                                                                                                                           >  2:914374/1‑237 (MQ=255)
                                        ttttCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTtgtg                                                                                                                                                                                                >  1:504380/1‑239 (MQ=255)
                                                   tgtGCCGTGTTGTCCGGCACGAAAAAACCCCCCGTGTTAAAATTGGGTCATACCAAGCGGAGGCGGGTGC‑T‑CCCCGGCCCACACCCGGAGGGGCTGCG‑CCACGCCATTACGATAGCCGACAGCGGGAAGTT‑CCCACGTTTGACCGTTAATGGTAAACCCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACg                                                                                                                                                                                     <  2:946063/237‑1 (MQ=255)
                                                         ttgttgtccggaccccAAAAAATCCCCGTCTTTAAATTTGTCCAGAGAAAACGGATG‑GCGGGGGTGCGCTGG‑CCCAAGGCGCAGCGGCGCCG‑CCACGCCATTACGAACGCCGACAACGGGGAGAC‑CCCACGCTTGACCGGTAATGGTAACCTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGt                                                                                                                                                                               <  2:1080478/224‑1 (MQ=255)
                                                            ctgTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCT‑‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGCCACTTTACCGCCAATCCCTACCTGCCCATCCAGGGCTTCGTTGTGGCGGAAACACGGGGAGtt                                                                                                                                                                              >  1:254809/1‑236 (MQ=255)
                                                               ttCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACCACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTCCACATCCAGGGCTTCGTTGTGTCGGAACCACgggagcgtaacg                                                                                                                                                                          >  2:84545/1‑229 (MQ=255)
                                                                    cACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCCCATCCAGTGCTTCGTTTTGGCGGAACCACCGGTAGTTACTGCGATGc                                                                                                                                                                    >  1:87422/1‑239 (MQ=255)
                                                                                   cgtcgtTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCgcgcct                                                                                                                                                     >  1:339001/1‑234 (MQ=255)
                                                                                      cgtTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTa                                                                                                                                                                             >  1:878642/1‑212 (MQ=255)
                                                                                       gtTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGGGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCa                                                                                                                                                  >  1:366131/1‑238 (MQ=255)
                                                                                                            cGGATGCGCGTGCGT‑CGCCCG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACaaaaa                                                                                                                            >  1:599125/1‑238 (MQ=255)
                                                                                                                           gccccTG‑CCACAGACGAAGCGGGTGCG‑CATCCCCATTACGAAAGCCGACAAGGGGAAACT‑CCCACGCTTGCCCGTTAATGGTAACCTCCGGCTCCCGGCTTCCGTCTTTCGTCGCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACaa                                                                                                               <  1:84545/235‑1 (MQ=255)
                                                                                                                            tccccgcccccacccggtgCGGCTCCG‑CCACCCGATTCCGAACGCGCACAAGGGGAAGTC‑CCCCAGCTGACCCGGTATGGGAAAACCCCGTCTCCCAGCGTCCTTCTTTCGTCCCTTTCCCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACaa                                                                                                               <  2:964904/220‑1 (MQ=255)
                                                                                                                                cggcacaagacccgcgggcggcC‑CCCCGCCATTCCGATACCCAAAAAGGGGGAGAT‑CCACGGTTTGACCGGTAAGGGTAAACTCCGGCTCCCAGCGTCCGTCTTGCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCat                                                                                                          <  1:901669/223‑1 (MQ=255)
                                                                                                                                   ccACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCCATGCCAGTCATCCGCCGCTTCCACCCGTTTCCCCTCACCaaaaaaacgcgggaaaaaaaccaatt                                                                                                       >  2:733942/1‑217 (MQ=255)
                                                                                                                                    ccaaagcgccagggcgtcgcccacgcaaTTACGATAGCCGAAAAGGGGGAAAT‑CCCACGCTTGACCGGTAATGGAAAACCCCGGCCCCCAGCGTCCGTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCat                                                                                                          <  1:1046998/216‑1 (MQ=255)
                                                                                                                                     accgccgcgggggCTGGC‑CCACGCCATTCCGATAGCCGACAACGGGGAGAC‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGa                                                                                                     <  1:914374/228‑1 (MQ=255)
                                                                                                                                      aaccgcgaGCGGCCGCGCCCCCGCCATTACAATAGCCGAAACGGGGGAGAC‑CCCCCTGTTGACCGGGAAGGGTAAACCCCGGCTCACAGCTCCCGTTTTTCGTCCCTTCACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGa                                                                                                     <  2:68258/232‑1 (MQ=255)
                                                                                                                                           gcgccgcgcgcg‑cCCCGCCTTTCCGATACCCGACACGGGGGATCT‑CCCCCGCTTTACCGGTAATGGTAAACCCCGCCCCCCAGCGTCTGTCTTCCGTCATTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTc                                                                                               <  2:339001/230‑1 (MQ=255)
                                                                                                                                            gaGCGGCTGCC‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTGCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCa                                                                                              <  1:804773/238‑1 (MQ=255)
                                                                                                                                             cacggggtccccccccccatcgcGAGAGCCGA‑AACGGGGGGGTCCCCACGTTTTACCGGTATGGGTAAACCCCGCCTCCCAGCGTCCTTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCa                                                                                              <  2:1190659/217‑1 (MQ=255)
                                                                                                                                             aGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAg                                                                                             >  1:68258/1‑239 (MQ=255)
                                                                                                                                             aGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCACGCGCTTCAACCCGTTTGCACTACAAAAAAAACCGGGGGAACAAAACATATTGaaaatt                                                                                                >  2:1046998/1‑231 (MQ=255)
                                                                                                                                                     cg‑ccacgccaTTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGACCCCCGGGTAATTACTGCGATGCCAGTCATCCGGGGCTTCAACCTTTTTGCCATCGACAAAAGACTCGCGGAACAAACCCTAATTAAAGTTCAAAcccgaa                                                                                      >  2:993866/1‑234 (MQ=255)
                                                                                                                                                      ggcccccgtggtgcGAGAGCCGACAAGGGGGAGCA‑CCCCCTCGTCGCCGTGAAGGACAATCCCCGGCTCCCCGCGTCCGTCTTTCGTCACGTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAAc                                                                                     <  2:101801/226‑1 (MQ=255)
                                                                                                                                                      ccccccccttttccGATAGCCCACAGCGGGGAGCC‑CCCCCGCTTGCCCGGTAATGGTAAACCCCGGCCCCCAGCGTCCGTCTTCCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAAc                                                                                     <  2:599125/226‑1 (MQ=255)
                                                                                                                                                       ccccccccaTTCCGAGAGCCGAAAACGGGGACTT‑CCCCGGCTTGCCCGTTAATGGAAAACTCCGGCCCCCGGGGCCCTTCTTCCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTGCTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAAc                                                                                     <  1:1106233/232‑1 (MQ=255)
                                                                                                                                                            gccaTTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCAt                                                                                                                                     <  1:1057032/185‑1 (MQ=255)
                                                                                                                                                            gccaTTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCAt                                                                                                                                     >  2:1057032/1‑185 (MQ=255)
                                                                                                                                                                gtACGATAGCCGACACGGGGGAAAC‑CCCCCGCTTGCCCGGTAATGGTAAACTCCGCCTCCCAGCCTCCTTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGa                                                                           <  1:774321/238‑1 (MQ=255)
                                                                                                                                                                   cagttgCCAACACCGGGGGGTC‑CCCACGGTTGTCCAGTAATGGTAACCCCCGTCCCCCGGCGTCCGTCTTCCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTTATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGATTg                                                                        <  2:87422/234‑1 (MQ=255)
                                                                                                                                                                      tAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATc                                                                                                                                                                                                               >  2:181825/1‑101 (MQ=255)
                                                                                                                                                                      tAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATc                                                                                                                                                                                                               <  1:181825/101‑1 (MQ=255)
                                                                                                                                                                       tacccggacggggggaaccccccCGGTTGACCGTTAGCTGTAACCCCCGGCCCCCTGCCTCCGTTTTCGGTCACTTTCCCGCCAATCCCTACCTGCACAGCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACGGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGAc                                                                      <  2:254809/219‑1 (MQ=255)
                                                                                                                                                                        gCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCCCAGACTGACCGa                                                                   >  1:537482/1‑239 (MQ=255)
                                                                                                                                                                        gCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCAAATTGTTAATTCAAACCGTAACCCGTCGCAGaccgaccg                                                                    >  1:841655/1‑238 (MQ=255)
                                                                                                                                                                            gaagcGGGGAATT‑CCCTCCCTTGAGCGGTAATGGTAACCCCCGGCCCCCAGCGTCCGTTTTTCGTCACTTTACCGCCAATCCCTACCTGCTCATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGt                                                               <  1:532278/235‑1 (MQ=255)
                                                                                                                                                                                   ggAGAT‑CCCCCGTTTCCCCGGTAAGGTTAAACCCCGCCCCCCCGTCTCCGTTTTCCGTCCTTTCACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCAt                                                         <  2:949179/238‑1 (MQ=255)
                                                                                                                                                                                   ggAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTTATAGTTCAGACCGTAACCCCTCGCAAACTGACCGCCAGTTTCCAtt                                                        >  1:911916/1‑239 (MQ=255)
                                                                                                                                                                                            aaggtgaaaaagagacgggTAAACCCCGGCCCCCGAGGCCCGTTTTCCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGa                                                 <  1:271828/221‑1 (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)
                                                                                                                                                                                                    aCCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACCAAAGACTGGCGGAACAAACCCTATTGATAGTTCATACCGTAACCCGTCTCACGCTGACCGCCCGTTTCCAtttgattcgtggagcca                                         >  1:591752/1‑223 (MQ=255)
                                                                                                                                                                                                        gTAATGGAAAACTCCGGCTCCCAGCGGCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCCCGGGTAGTTACTGCGATTCCAGTCATCCGGCGCTTCAACCTGTTTGCCCACGTCAAAAGCATGGCGGAACAAACCCAATTAGTAATTACGACCGTAACCCGTGCCAGACTGAACGACAGTTGACATTTAGTCggaggagaacggccc                                    >  2:905260/1‑225 (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)
                                                                                                                                                                                                                        gCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACCAAAGACTGGCGGAACAAAACATATTGATAGTTTCGACCGTAACCCGTCGCAGCGTGACCCCAAGTTGCCCTTTGGTCGAGGAAGCACCCCGCCAGagcacccccacca                    >  2:793399/1‑227 (MQ=255)
                                                                                                                                                                                                                                      tcTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCCTATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCCAACGTCCCAAACCACCGTTACCCAgcg       >  1:637231/1‑238 (MQ=255)
                                                                                                                                                                                                                                       tttccGTCCCTTTACCGCCAATCCCTGCCGGTACAACCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAgcgc      <  1:1275855/234‑1 (MQ=255)
                                                                                                                                                                                                                                          tCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGCCCGACAGTTGCCATTTAGTCGAGGAAGAACGCCGCCAGACGTCCCACACCACCGTTACCCCGCGCCggg  >  1:1360183/1‑239 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
GCAGTATGGTTGTCGTTTGGATAGAGAACTTTGG‑TCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGT‑CGCCTG‑CCACAGACGCAGCGGCTGCG‑CCACGCCATTACGATAGCCGACAACGGGGAGAT‑CCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCGCCGGG  >  NC_000913/3551681‑3552148

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