Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 812,778 T→C 5.9% L170P (CTG→CCG)  bioD → dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913812,7780TC 5.9% 315.5 / 7.2 102L170P (CTG→CCG) bioDdethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base T (51/45);  new base C (2/4);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 105 total aligned reads displayed.

CTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACA  >  NC_000913/812629‑812927
                                                                                                                                                     |                                                                                                                                                     
cTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGActct                                                                                                                                                       <  1:561772/150‑1 (MQ=255)
   gTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGc                                                                                                                                                    >  1:1192387/1‑150 (MQ=255)
    tttACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCg                                                                                                                                                   <  1:1571191/150‑1 (MQ=255)
     ttACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCgg                                                                                                                                                  >  2:1090706/1‑150 (MQ=255)
      tACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCggg                                                                                                                                                 >  2:460530/1‑150 (MQ=255)
       aCGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGt                                                                                                                                                >  2:1349029/1‑150 (MQ=255)
        cgccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGtt                                                                                                                                               >  1:1101710/1‑150 (MQ=255)
         gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACTCAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCAGAGTCAGTCCGGCGTGTTg                                                                                                                                              <  2:1087255/150‑1 (MQ=255)
         gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGGAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACCCGCCAGAGTCAGTCCGGCGTGTTg                                                                                                                                              <  1:78020/150‑1 (MQ=255)
         gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTg                                                                                                                                              <  2:1192387/150‑1 (MQ=255)
         gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTg                                                                                                                                              <  1:8434/150‑1 (MQ=255)
           cgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTggg                                                                                                                                            >  1:201103/1‑150 (MQ=255)
           cgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTggg                                                                                                                                            <  1:1062120/150‑1 (MQ=255)
             cTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACCCGCCAGAGTCAGTCCGGCGTGTTgtatt                                                                                                                                          <  2:1395987/150‑5 (MQ=255)
              tttCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGtag                                                                                                                                         <  2:1494130/150‑3 (MQ=255)
                tCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCg                                                                                                                                       <  2:1207584/150‑1 (MQ=255)
                  tGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGaa                                                                                                                                     >  2:395578/1‑150 (MQ=255)
                   gACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAAc                                                                                                                                    <  1:1090706/150‑1 (MQ=255)
                    acacTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACg                                                                                                                                   <  1:1549275/150‑1 (MQ=255)
                         ttcactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGtt                                                                                                                              >  1:349630/1‑150 (MQ=255)
                           cactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTAc                                                                                                                            <  1:194927/150‑1 (MQ=255)
                           cactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTAc                                                                                                                            >  2:1629419/1‑150 (MQ=255)
                            actttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGCTACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACg                                                                                                                           <  1:1659416/150‑1 (MQ=255)
                            actttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACg                                                                                                                           >  1:1436919/1‑150 (MQ=255)
                             ctttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGc                                                                                                                          <  1:1134300/150‑1 (MQ=255)
                              ttttGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGcc                                                                                                                         <  2:1177004/150‑1 (MQ=255)
                                  gCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCg                                                                                                                     >  1:1536836/1‑150 (MQ=255)
                                     gATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGa                                                                                                                  >  2:1266916/1‑150 (MQ=255)
                                         gggTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAAc                                                                                                              >  2:1261705/1‑150 (MQ=255)
                                         gcgTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAAc                                                                                                              >  1:1041069/3‑150 (MQ=255)
                                           gTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGt                                                                                                            <  1:202306/150‑1 (MQ=255)
                                           gTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGt                                                                                                            >  2:1328604/1‑150 (MQ=255)
                                            tAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTc                                                                                                           <  2:1436919/150‑1 (MQ=255)
                                            tAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTc                                                                                                           >  2:1558408/1‑150 (MQ=255)
                                             aaCACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCa                                                                                                          >  2:558534/1‑150 (MQ=255)
                                              acacAGGAACAACTGCCGGTGATACTGGTAGTTGTTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                                         <  1:422438/150‑1 (MQ=255)
                                              acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                                         <  1:915711/150‑1 (MQ=255)
                                              acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                                         <  1:921651/150‑1 (MQ=255)
                                              acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                                         <  2:1297826/150‑1 (MQ=255)
                                              acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTCTATTAATCACGAAATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                                         >  2:1651722/1‑150 (MQ=255)
                                               cacaGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACg                                                                                                        <  2:1019692/150‑1 (MQ=255)
                                                acaGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGc                                                                                                       >  2:1536371/1‑150 (MQ=255)
                                                 caGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCt                                                                                                      <  2:235090/150‑1 (MQ=255)
                                                    gAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGaa                                                                                                   <  1:1261705/150‑1 (MQ=255)
                                                     aacaacTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAat                                                                                                  <  1:646643/150‑1 (MQ=255)
                                                         acTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGGCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatatat                                                                                              <  2:457368/150‑1 (MQ=255)
                                                            gCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGAt                                                                                                                                 >  2:352012/1‑112 (MQ=255)
                                                            gCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGAt                                                                                                                                 <  1:352012/112‑1 (MQ=255)
                                                              cGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAGCACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCa                                                                                         <  2:762916/150‑1 (MQ=255)
                                                                    tACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCa                                                                                   >  2:299131/1‑150 (MQ=255)
                                                                      cTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCAcc                                                                                 <  1:393023/150‑1 (MQ=255)
                                                                       tGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCAccc                                                                                <  1:599303/150‑1 (MQ=255)
                                                                        ggTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCg                                                                               >  1:865298/1‑150 (MQ=255)
                                                                            gTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATg                                                                           <  1:1349029/150‑1 (MQ=255)
                                                                             ttGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGCCCACGCTCACCCGCATGa                                                                          >  2:1389307/1‑150 (MQ=255)
                                                                             ttGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGa                                                                          >  2:30653/1‑150 (MQ=255)
                                                                                 tgtgAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTcc                                                                      >  1:297573/1‑150 (MQ=255)
                                                                                    gAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCcgc                                                                   >  1:622000/1‑150 (MQ=255)
                                                                                    gAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCcgc                                                                   >  1:772952/1‑150 (MQ=255)
                                                                                       aCTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGcgcc                                                                <  1:395578/150‑1 (MQ=255)
                                                                                           ggCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCgctg                                                            >  1:1374638/1‑150 (MQ=255)
                                                                                            gCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCgctgc                                                           >  2:1024742/1‑150 (MQ=255)
                                                                                             cTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatat                                                                                                <  2:590449/112‑1 (MQ=255)
                                                                                             cTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatat                                                                                                >  1:590449/1‑112 (MQ=255)
                                                                                                    aaTCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGTGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGgagag                                                   >  1:652999/1‑150 (MQ=255)
                                                                                                    aaTCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGgagag                                                   <  2:349630/150‑1 (MQ=255)
                                                                                                        aCGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATcc                                               >  1:1438000/1‑150 (MQ=255)
                                                                                                            gATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTg                                           <  2:806910/150‑1 (MQ=255)
                                                                                                             aTGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTgg                                          >  1:1407072/1‑150 (MQ=255)
                                                                                                              tGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGc                                         <  1:558534/150‑1 (MQ=255)
                                                                                                               gTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCt                                        <  1:1024742/150‑1 (MQ=255)
                                                                                                                 tGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTg                                      >  1:1366255/1‑150 (MQ=255)
                                                                                                                   aCTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCa                                    <  2:772952/150‑1 (MQ=255)
                                                                                                                     tGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTCCGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCGGCTGGGAGAGATCCCCTGGCTTGCAGa                                  <  2:409490/150‑1 (MQ=255)
                                                                                                                     tGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTAACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGa                                  <  2:1536836/150‑1 (MQ=255)
                                                                                                                     ggCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGa                                  <  2:415769/149‑1 (MQ=255)
                                                                                                                      gCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaa                                 >  1:1615472/1‑150 (MQ=255)
                                                                                                                      gCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaa                                 >  1:272781/1‑150 (MQ=255)
                                                                                                                       cacaGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                                >  2:1652533/1‑150 (MQ=255)
                                                                                                                       cacaGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                                >  2:1488372/1‑150 (MQ=255)
                                                                                                                       cacaGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                                >  2:232426/1‑150 (MQ=255)
                                                                                                                         caGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAAt                              >  1:1151636/1‑150 (MQ=255)
                                                                                                                           ggTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATcc                            >  1:192898/1‑150 (MQ=255)
                                                                                                                           ggTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATcc                            >  1:1215448/1‑150 (MQ=255)
                                                                                                                           ggTAATACAACACGCCGGACTGACTCTGCCGGGTAGGCTGGCGAACGATGTTACGACTCCCGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCAGCTGCTCGCAGACATCCCCTGGCATGCAGAAAATcc                            >  1:429715/1‑150 (MQ=255)
                                                                                                                              aaTACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGa                         >  2:171881/1‑150 (MQ=255)
                                                                                                                               acccaaccCGCCAGAGTCAGTCCGGCGTGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATc                             <  2:507502/142‑1 (MQ=255)
                                                                                                                               acccaaccCGCCAGAGTCAGTCCGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGaa                        >  2:210535/4‑150 (MQ=255)
                                                                                                                                 acaacaCGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGaaaa                      >  1:382626/1‑150 (MQ=255)
                                                                                                                                  caacaCGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAAt                     >  1:132186/1‑150 (MQ=255)
                                                                                                                                     cacGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCg                  <  1:1286381/150‑1 (MQ=255)
                                                                                                                                     cacGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCg                  <  1:117438/150‑1 (MQ=255)
                                                                                                                                      acGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCgg                 <  1:1558408/150‑1 (MQ=255)
                                                                                                                                      acGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCgg                 >  2:656295/1‑150 (MQ=255)
                                                                                                                                       cGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGc                >  2:1554578/1‑150 (MQ=255)
                                                                                                                                           ggACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAAcc            >  1:817344/1‑150 (MQ=255)
                                                                                                                                               tgactCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGaa        <  1:299131/150‑1 (MQ=255)
                                                                                                                                               tgactCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGaa        <  1:1200931/150‑1 (MQ=255)
                                                                                                                                                  ctctGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGt     <  2:1369345/150‑1 (MQ=255)
                                                                                                                                                     tGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACa  <  2:201103/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACA  >  NC_000913/812629‑812927

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: