Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,419,395 C→G 6.3% V132L (GTG→CTG)  yfcE ← phosphodiesterase activity on bis‑pNPP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,419,3950CG 6.3% 626.4 / 26.5 174V132L (GTG→CTG) yfcEphosphodiesterase activity on bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (77/86);  new base G (5/6);  total (82/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

Only 100 of 177 total aligned reads displayed.

GACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATGCAG  >  NC_000913/2419249‑2419542
                                                                                                                                                  |                                                                                                                                                   
      cgcAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCAGTGGCAGaa                                                                                                                                             <  1:1459020/149‑2 (MQ=255)
        cAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATgg                                                                                                                                           <  2:848342/149‑1 (MQ=255)
           tGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATg                                                                                                                                            >  2:564747/1‑145 (MQ=255)
            gATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTAt                                                                                                                                       >  4:1089331/1‑149 (MQ=255)
              tACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGa                                                                                                                                     >  2:1072587/1‑149 (MQ=255)
                 tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATg                                                                                                                                            >  4:1350612/1‑139 (MQ=255)
                 tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATg                                                                                                                                            <  3:1350612/139‑1 (MQ=255)
                 tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGAc                                                                                                                                    <  2:317105/147‑1 (MQ=255)
                 tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCg                                                                                                                                  >  2:1272571/1‑149 (MQ=255)
                     aTCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATgg                                                                                                                                           <  2:360511/136‑1 (MQ=255)
                     aTCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTaca                                                                                                                              >  1:1423365/1‑149 (MQ=255)
                        aTTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCa                                                                                                                           >  4:1286851/1‑149 (MQ=255)
                             gTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGcaca                                                                                                                      >  3:393974/1‑149 (MQ=255)
                                         tAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGttt                                                                                                          >  3:100221/1‑149 (MQ=255)
                                         tAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGttt                                                                                                          >  3:366406/1‑149 (MQ=255)
                                          aaGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTa                                                                                                         <  1:835765/149‑1 (MQ=255)
                                            gTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaaa                                                                                                       >  1:507477/1‑149 (MQ=255)
                                            gTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaaa                                                                                                       <  2:1423365/149‑1 (MQ=255)
                                            gTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaaa                                                                                                       <  4:393974/149‑1 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGcac                                                                                                                       >  4:654322/1‑129 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGcac                                                                                                                       <  3:654322/129‑1 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCt                                                                                                               <  3:947603/137‑1 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCt                                                                                                               >  4:947603/1‑137 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   >  2:506592/1‑149 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTCCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGa                                                                                                                                              >  1:212793/1‑106 (MQ=255)
                                                 tCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGtt                                                                                                           <  4:1410788/140‑1 (MQ=255)
                                                   aTTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   <  1:1168685/146‑1 (MQ=255)
                                                   aTTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   >  2:1168685/1‑146 (MQ=255)
                                                       tCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAtt                                                                                            >  1:254490/1‑149 (MQ=255)
                                                         cAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAtttt                                                                                          >  1:1290168/1‑149 (MQ=255)
                                                           gCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTcc                                                                                        <  1:653085/149‑1 (MQ=255)
                                                           gCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTcc                                                                                        <  2:254490/149‑1 (MQ=255)
                                                                 cATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATgg                                                                                                                                           <  4:470777/92‑1 (MQ=255)
                                                                 cATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCa                                                                                  >  1:360916/1‑149 (MQ=255)
                                                                  aTAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATgg                                                                                                                                           >  1:660253/1‑91 (MQ=255)
                                                                     aCTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGTAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaaa                                                                              <  4:100221/149‑1 (MQ=255)
                                                                     aCTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaaa                                                                              <  2:507477/149‑1 (MQ=255)
                                                                      cTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAg                                                                             <  3:1390382/149‑1 (MQ=255)
                                                                       tCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGa                                                                            <  2:1405475/149‑1 (MQ=255)
                                                                          ccGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGa                                                                         <  3:1089331/149‑1 (MQ=255)
                                                                             gATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCGCCCCGCTGTTCTGCCACTGGCAGATGGGt                                                                                                                                         >  4:1363363/1‑82 (MQ=255)
                                                                              aTTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCAt                                                                     >  4:221575/1‑149 (MQ=255)
                                                                                tACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGa                                                                            >  3:429666/1‑140 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCAGTGGCAGaa                                                                                                                                             <  3:902201/73‑2 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCAGTGGCAGaa                                                                                                                                             >  4:902201/1‑72 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCAGTGGCAGaa                                                                                                                                             >  2:1149188/1‑72 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaaa                                                                                                       >  1:1411740/1‑111 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGt                                                                 >  2:241163/1‑149 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGt                                                                 >  2:919327/1‑149 (MQ=255)
                                                                                   cGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGAt                                                                           <  2:472205/138‑1 (MQ=255)
                                                                                    gCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGAc                                                                        >  3:585030/1‑140 (MQ=255)
                                                                                    gCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGAc                                                                        >  1:929062/1‑140 (MQ=255)
                                                                                     ccTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACt                                                                                                                                                    >  1:293065/1‑63 (MQ=255)
                                                                                        ttCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCaaaaa                                                           <  1:1150109/149‑1 (MQ=255)
                                                                                         tCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaaa                                                                              >  4:1115168/1‑129 (MQ=255)
                                                                                         tCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaaa                                                                              <  3:1115168/129‑1 (MQ=255)
                                                                                            gAATACTCGCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGAc                                                                        <  3:1171431/132‑1 (MQ=255)
                                                                                            gAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGAc                                                                        >  4:1171431/1‑132 (MQ=255)
                                                                                                  tCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttgt                                                 <  1:617573/149‑1 (MQ=255)
                                                                                                   cACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttgtt                                                <  2:360916/149‑1 (MQ=255)
                                                                                                    aCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTc                                                                                         <  1:356052/107‑1 (MQ=255)
                                                                                                    aCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttt                                               <  3:1417610/149‑1 (MQ=255)
                                                                                                    aCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttt                                               <  3:221575/149‑1 (MQ=255)
                                                                                                     ccGAGCCGGCGTTGAAGCGGTACATCACCCCCCGCTGTTCTGCTACAGGCAGATGGGTATTACCGTACACCGACACATCGTTCTTGTTTAAAGCAGGCAGATTTTCCGGGCCAACATGAGCACCATCCGTCAAAAACAGACGGTGtttt                                              <  4:584943/149‑1 (MQ=255)
                                                                                                       gAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTc                                            <  1:1018261/149‑1 (MQ=255)
                                                                                                       gAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTc                                            <  2:1314292/149‑1 (MQ=255)
                                                                                                             gggTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   <  1:475380/88‑1 (MQ=255)
                                                                                                             gggTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTc                                                                                         <  1:469990/98‑1 (MQ=255)
                                                                                                             gggTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTaa                                      <  1:919327/149‑1 (MQ=255)
                                                                                                                ttGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATAc                                   <  3:1286851/149‑1 (MQ=255)
                                                                                                                     gTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGtt                              >  3:1017438/1‑149 (MQ=255)
                                                                                                                       ggAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGc                            >  1:1081002/1‑149 (MQ=255)
                                                                                                                       ggAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGc                            >  2:877653/1‑149 (MQ=255)
                                                                                                                        gAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttt                                               <  2:1366619/129‑1 (MQ=255)
                                                                                                                        gAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttt                                               >  1:1366619/1‑129 (MQ=255)
                                                                                                                        gAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCa                                          <  2:16708/134‑1 (MQ=255)
                                                                                                                        gAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATa                                    >  4:363101/1‑140 (MQ=255)
                                                                                                                          aaaTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCggg                                                                                     >  1:727887/1‑89 (MQ=255)
                                                                                                                          aaaTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCggg                                                                                     <  2:727887/89‑1 (MQ=255)
                                                                                                                            aTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGca                                                                                                                        <  3:1222160/52‑1 (MQ=255)
                                                                                                                            aTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGca                                                                                                                        >  4:1222160/1‑52 (MQ=255)
                                                                                                                              ctcCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCa                                          <  4:373329/128‑1 (MQ=255)
                                                                                                                                ccccccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaaa                                                                              >  3:576493/1‑90 (MQ=255)
                                                                                                                                    ccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTc                                            >  1:104420/1‑120 (MQ=255)
                                                                                                                                    ccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTc                                            <  2:104420/120‑1 (MQ=255)
                                                                                                                                    ccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTa               >  4:1237198/1‑149 (MQ=255)
                                                                                                                                       ttGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCg            <  2:261093/148‑1 (MQ=255)
                                                                                                                                         catCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAg                                                                             >  2:1135807/3‑82 (MQ=255)
                                                                                                                                         catCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAg                                                                             <  1:1135807/80‑1 (MQ=255)
                                                                                                                                          ttCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGa         >  1:308288/1‑149 (MQ=255)
                                                                                                                                          atCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGa                                                       <  4:1307071/102‑1 (MQ=255)
                                                                                                                                          atCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGa                                                       <  3:345977/102‑1 (MQ=255)
                                                                                                                                          atCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGa                                                       >  3:1307071/2‑103 (MQ=255)
                                                                                                                                            cTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGaaa       <  4:1017438/149‑1 (MQ=255)
                                                                                                                                             tGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATc             <  2:435959/142‑1 (MQ=255)
                                                                                                                                             tGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATc             >  1:435959/1‑142 (MQ=255)
                                                                                                                                                cACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATAc                                   <  4:599184/117‑1 (MQ=255)
                                                                                                                                                cACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGgcgc                    <  2:65945/132‑1 (MQ=255)
                                                                                                                                                cACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGt                 >  3:758762/1‑135 (MQ=255)
                                                                                                                                                cACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGt                 <  4:758762/135‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATGCAG  >  NC_000913/2419249‑2419542

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