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

GCAGTATGGTTGTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCGCCGGG  >  NC_000913/3551681‑3552148
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
ggttttttttttttgTTTGGATAGAAACTTTGGCGCAGTTTTCCCCATTTGACCCCCTCTTGTCCGGCCCGCAAAAAATCCCCGTCGTTAAATTTAGTCAGATCAACGGGCGGCGCGTGGGGCGCCTGCCACAGCCCCGGCGCCTGCCCCCGCCCATTACGCTAGCCGCCAACGGGGAGACCCCCCGTTTGACCGGTAATGGTAACCTCCGGCTCCCAGCGTCCGTCTTTCGTCACttt                                                                                                                                                                                                                                       <  2:35991/225‑1 (MQ=255)
           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)
                                                         ttgttgtccggaccccAAAAAATCCCCGTCTTTAAATTTGTCCAGAGAAAACGGATGGCGGGGGTGCGCTGGCCCAAGGCGCAGCGGCGCCGCCACGCCATTACGAACGCCGACAACGGGGAGACCCCACGCTTGACCGGTAATGGTAACCTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGt                                                                                                                                                                               <  2:1080478/224‑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)
                                                                                                                            tccccgcccccacccggtgCGGCTCCGCCACCCGATTCCGAACGCGCACAAGGGGAAGTCCCCCAGCTGACCCGGTATGGGAAAACCCCGTCTCCCAGCGTCCTTCTTTCGTCCCTTTCCCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACaa                                                                                                               <  2:964904/220‑1 (MQ=255)
                                                                                                                                cggcacaagacccgcgggcggcCCCCCGCCATTCCGATACCCAAAAAGGGGGAGATCCACGGTTTGACCGGTAAGGGTAAACTCCGGCTCCCAGCGTCCGTCTTGCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCat                                                                                                          <  1:901669/223‑1 (MQ=255)
                                                                                                                                   ccACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCCATGCCAGTCATCCGCCGCTTCCACCCGTTTCCCCTCACCaaaaaaacgcgggaaaaaaaccaatt                                                                                                       >  2:733942/1‑217 (MQ=255)
                                                                                                                                    ccaaagcgccagggcgtcgcccacgcaaTTACGATAGCCGAAAAGGGGGAAATCCCACGCTTGACCGGTAATGGAAAACCCCGGCCCCCAGCGTCCGTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCat                                                                                                          <  1:1046998/216‑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)
                                                                                                                                             cacggggtccccccccccatcgcGAGAGCCGAAACGGGGGGGTCCCCACGTTTTACCGGTATGGGTAAACCCCGCCTCCCAGCGTCCTTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCa                                                                                              <  2:1190659/217‑1 (MQ=255)
                                                                                                                                             aGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAg                                                                                             >  1:68258/1‑239 (MQ=255)
                                                                                                                                             aGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCACGCGCTTCAACCCGTTTGCACTACAAAAAAAACCGGGGGAACAAAACATATTGaaaatt                                                                                                >  2:1046998/1‑231 (MQ=255)
                                                                                                                                                     cgccacgccaTTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGACCCCCGGGTAATTACTGCGATGCCAGTCATCCGGGGCTTCAACCTTTTTGCCATCGACAAAAGACTCGCGGAACAAACCCTAATTAAAGTTCAAAcccgaa                                                                                      >  2:993866/1‑234 (MQ=255)
                                                                                                                                                      ggcccccgtggtgcGAGAGCCGACAAGGGGGAGCACCCCCTCGTCGCCGTGAAGGACAATCCCCGGCTCCCCGCGTCCGTCTTTCGTCACGTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAAc                                                                                     <  2:101801/226‑1 (MQ=255)
                                                                                                                                                      ccccccccttttccGATAGCCCACAGCGGGGAGCCCCCCCGCTTGCCCGGTAATGGTAAACCCCGGCCCCCAGCGTCCGTCTTCCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAAc                                                                                     <  2:599125/226‑1 (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                                                                                                                                                                                                               >  2:181825/1‑101 (MQ=255)
                                                                                                                                                                      tAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATc                                                                                                                                                                                                               <  1:181825/101‑1 (MQ=255)
                                                                                                                                                                       tacccggacggggggaaccccccCGGTTGACCGTTAGCTGTAACCCCCGGCCCCCTGCCTCCGTTTTCGGTCACTTTCCCGCCAATCCCTACCTGCACAGCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACGGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGAc                                                                      <  2:254809/219‑1 (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)
                                                                                                                                                                                            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)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
GCAGTATGGTTGTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAGATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGCTTCAACCTGTTTGCCATCGACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCGCCGGG  >  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: