Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 423,990 C→G 5.8% N491K (AAC→AAG) ‡ malZ → maltodextrin glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913423,9900CG 5.8% 531.5 / 16.3 156N491K (AAC→AAG) ‡malZmaltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base C (69/77);  new base G (4/5);  total (73/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 161 total aligned reads displayed.

TAATCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGT  >  NC_000913/423840‑424139
                                                                                                                                                      |                                                                                                                                                     
tAATCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATACGGTGATGAAGTAGGACTGGATGGCAAAAAc                                                                                                                                                       <  2:2439709/150‑1 (MQ=255)
 aaTCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAAc                                                                                                                                                       <  1:2333854/150‑1 (MQ=255)
 aaTCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAAc                                                                                                                                                       <  2:2403856/150‑1 (MQ=255)
  aTCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACg                                                                                                                                                      >  2:1889699/1‑150 (MQ=255)
    cAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGAt                                                                                                                                                    >  2:1232463/1‑150 (MQ=255)
     aGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATc                                                                                                                                                   <  2:664777/150‑1 (MQ=255)
     aGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATc                                                                                                                                                   >  1:1429940/1‑150 (MQ=255)
            cAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTg                                                                                                                                            >  2:624987/1‑150 (MQ=255)
             aGCCACGATACTGCGGGATTTAAAACGCTGGTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGc                                                                                                                                           <  2:480802/150‑1 (MQ=255)
             aGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGc                                                                                                                                           <  2:322495/150‑1 (MQ=255)
             aGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGc                                                                                                                                           <  1:1705317/150‑1 (MQ=255)
               ccACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCg                                                                                                                                         >  1:1002510/1‑150 (MQ=255)
               ccACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCg                                                                                                                                         >  2:202423/1‑150 (MQ=255)
                cACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCACTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGccat                                                                                                                                        <  1:1079698/150‑3 (MQ=255)
                 aCGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGccatc                                                                                                                                       <  1:1461392/150‑4 (MQ=255)
                 aCGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGccatc                                                                                                                                       <  1:389875/150‑4 (MQ=255)
                    aTACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAAc                                                                                                                                    <  2:2428470/150‑1 (MQ=255)
                     tACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGccatccagt                                                                                                                                   <  1:92815/150‑8 (MQ=255)
                       cTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGt                                                                                                                                 <  1:556/150‑1 (MQ=255)
                         gcgcGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGAAGGCAAAAACGATCCGTTTTGCCGTAAACCGTTc                                                                                                                               >  1:2305979/1‑150 (MQ=255)
                          cgcgATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcc                                                                                                                              >  2:2072570/1‑150 (MQ=255)
                          cgcgATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcc                                                                                                                              >  2:505927/1‑150 (MQ=255)
                           gcgATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTccc                                                                                                                             >  1:2439311/1‑150 (MQ=255)
                              aTTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTg                                                                                                                          <  1:2430660/150‑1 (MQ=255)
                                ttAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGc                                                                                                                        <  1:671186/150‑1 (MQ=255)
                                     aCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTg                                                                                                                   >  2:1325111/1‑150 (MQ=255)
                                      cGCTGCTCGGGCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTgg                                                                                                                  <  2:733940/150‑1 (MQ=255)
                                          gctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTgg                                                                                                                  <  1:946257/146‑1 (MQ=255)
                                          gctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTgg                                                                                                                  >  2:946257/1‑146 (MQ=255)
                                          gctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGaaaa                                                                                                              >  2:1621919/1‑150 (MQ=255)
                                          gctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGaaaa                                                                                                              >  2:1559624/1‑150 (MQ=255)
                                           ctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAg                                                                                                             <  2:1379217/150‑1 (MQ=255)
                                            tcggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGc                                                                                                            <  2:2310056/150‑1 (MQ=255)
                                              ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                          <  1:1089628/150‑1 (MQ=255)
                                              ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                          >  2:1149433/1‑150 (MQ=255)
                                              ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                          <  1:1536078/150‑1 (MQ=255)
                                              ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                          >  1:1830218/1‑150 (MQ=255)
                                               gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAgg                                                                                                         <  1:1567001/150‑1 (MQ=255)
                                                    gATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCATGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACg                                                                                                    >  2:1921768/1‑150 (MQ=255)
                                                        ttGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGt                                                                                                <  1:1871295/150‑1 (MQ=255)
                                                        ttGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGt                                                                                                <  1:1292644/150‑1 (MQ=255)
                                                        ttGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGt                                                                                                <  1:15502/150‑1 (MQ=255)
                                                        ttGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGt                                                                                                <  1:163491/150‑1 (MQ=255)
                                                            gcgcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTAtt                                                                                            >  1:80279/1‑150 (MQ=255)
                                                              gcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTcg                                                                                          <  1:1983108/150‑1 (MQ=255)
                                                                 tGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCgcgc                                                                                       <  1:1657076/150‑1 (MQ=255)
                                                                 tGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCgcgc                                                                                       <  2:1418539/150‑1 (MQ=255)
                                                                 tGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCgcgc                                                                                       <  2:171297/150‑1 (MQ=255)
                                                                       tGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTAcc                                                                                 <  1:838132/150‑1 (MQ=255)
                                                                        ggCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCa                                                                                <  1:928016/150‑1 (MQ=255)
                                                                         gCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAg                                                                               <  2:1925884/150‑1 (MQ=255)
                                                                           ggtggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCg                                                                             <  1:1564800/150‑1 (MQ=255)
                                                                             tggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                           <  2:112790/150‑1 (MQ=255)
                                                                              ggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAAt                                                                          <  2:537221/150‑1 (MQ=255)
                                                                              ggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAAt                                                                          >  2:1516905/1‑150 (MQ=255)
                                                                                   ggCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTg                                                                     >  2:738899/1‑150 (MQ=255)
                                                                                   ggCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTg                                                                     <  1:1547443/150‑1 (MQ=255)
                                                                                   ggCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTg                                                                     <  2:2305979/150‑1 (MQ=255)
                                                                                   ggCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTg                                                                     >  2:828236/1‑150 (MQ=255)
                                                                                    gCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTgc                                                                    <  2:2011392/150‑1 (MQ=255)
                                                                                        ttCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGAACCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTg                                                                <  2:1769844/150‑1 (MQ=255)
                                                                                             ctggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTaa                                                           <  2:80279/150‑1 (MQ=255)
                                                                                                   tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAg                                                     <  1:1468168/150‑1 (MQ=255)
                                                                                                   tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAg                                                     <  2:1455068/150‑1 (MQ=255)
                                                                                                    ggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGt                                                    <  1:1472656/150‑1 (MQ=255)
                                                                                                     gtgtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTc                                                   <  2:1830218/150‑1 (MQ=255)
                                                                                                          ccGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGgcg                                              >  1:451232/1‑150 (MQ=255)
                                                                                                           cGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAg                                                                                                             >  2:1952010/1‑86 (MQ=255)
                                                                                                           cGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGgc                                               <  1:2137218/148‑1 (MQ=255)
                                                                                                              gCATTTATTACGGTGATTAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTAc                                          <  2:858133/150‑1 (MQ=255)
                                                                                                              gCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTAc                                          >  1:127415/1‑150 (MQ=255)
                                                                                                                 tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTc                                       >  1:1622303/1‑150 (MQ=255)
                                                                                                                   tattaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTcg                                                                                          <  1:2052913/97‑1 (MQ=255)
                                                                                                                   tattaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTcg                                                                                          >  2:2052913/1‑97 (MQ=255)
                                                                                                                   tattaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCAt                                     >  1:1724094/1‑150 (MQ=255)
                                                                                                                    attaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGt                                                             >  1:696604/1‑125 (MQ=255)
                                                                                                                    attaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGt                                                             <  2:696604/125‑1 (MQ=255)
                                                                                                                       aCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                 >  1:1661938/1‑150 (MQ=255)
                                                                                                                       aCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                 >  2:1595109/1‑150 (MQ=255)
                                                                                                                       aCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                 >  1:342636/1‑150 (MQ=255)
                                                                                                                       aCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                 >  2:1402327/1‑150 (MQ=255)
                                                                                                                             atgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                           <  1:606042/150‑1 (MQ=255)
                                                                                                                             atgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                           >  1:1131670/1‑150 (MQ=255)
                                                                                                                              tgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCa                          >  1:280712/1‑150 (MQ=255)
                                                                                                                              tgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCa                          >  1:2250128/1‑150 (MQ=255)
                                                                                                                              tgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCa                          >  2:1687073/1‑150 (MQ=255)
                                                                                                                                   gaGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGc                     <  1:2206516/149‑1 (MQ=255)
                                                                                                                                    aGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCt                    >  2:2207523/1‑150 (MQ=255)
                                                                                                                                       aCTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTa                 <  1:1232463/150‑1 (MQ=255)
                                                                                                                                          ggATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGc              >  1:2171753/1‑150 (MQ=255)
                                                                                                                                          ggATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGc              <  2:239688/150‑1 (MQ=255)
                                                                                                                                          ggATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGc              <  2:1082690/150‑1 (MQ=255)
                                                                                                                                           tacggCAAAACGGATCGTTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGGTACGTCATGGCGGCTGTCAGGTGCTGTATGCg             >  2:1655249/4‑150 (MQ=255)
                                                                                                                                            aTGGCAAAAACGATCCGTTTTTCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCgg            <  1:1149433/150‑1 (MQ=255)
                                                                                                                                             tGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGa           <  1:828236/150‑1 (MQ=255)
                                                                                                                                              gtaaaaaaCGATCCGTTTTGCCGTAAACCGTTCCCATGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGaa          >  2:1635995/4‑150 (MQ=255)
                                                                                                                                              ggCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGaa          <  1:1484116/150‑1 (MQ=255)
                                                                                                                                               gCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAg         >  2:456814/1‑150 (MQ=255)
                                                                                                                                                cAAAAACGATCCGTTTTGCCGTAAATCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGa        >  2:2315057/1‑150 (MQ=255)
                                                                                                                                                      cGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACgt  >  1:1234738/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TAATCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGT  >  NC_000913/423840‑424139

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