Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,555 Δ1 bp 5.3% coding (1126/1452 nt) uxaB ← altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,609,5550G. 5.3% 376.4 / 13.8 114coding (1126/1452 nt)uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (60/48);  new base . (4/2);  total (64/50)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 120 total aligned reads displayed.

CCAGAGTTGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGG  >  NC_000913/1609406‑1609702
                                                                                                                                                     |                                                                                                                                                    
ccAGAGTTGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcag                                                                                                                                                      >  2:31045/1‑150 (MQ=255)
ccAGAGTTGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcag                                                                                                                                                      >  2:1821867/1‑150 (MQ=255)
 cAGAGTTGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagc                                                                                                                                                     <  2:651434/150‑1 (MQ=255)
 cAGAGTTGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagc                                                                                                                                                     <  1:1348802/150‑1 (MQ=255)
  agagTTGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGCTAGCACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagcc                                                                                                                                                    >  2:1605429/1‑149 (MQ=255)
       tGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑t                                                                                                                                              <  1:1790891/150‑1 (MQ=255)
        gCTGAAAACGGACCAGCCAGTGTGCATCGTCCGGTACCGTATATGCTTCGCCGTTACGCTCACCACGATAGAACGCAATCAAAGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                                                                             <  2:350941/150‑1 (MQ=255)
         cTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑Tgg                                                                                                                                            >  2:796626/1‑150 (MQ=255)
         cTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑Tgg                                                                                                                                            <  1:1554678/150‑1 (MQ=255)
         cTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑Tgg                                                                                                                                            >  1:758295/1‑150 (MQ=255)
          ggATACCGTTCCAGCCCGTGTGCATCGTCCTGTACCGGATATGGTTCGCCGTTACGCTCACCACGATAGACCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGc                                                                                                                                           <  1:335162/149‑1 (MQ=255)
           gATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGTCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCa                                                                                                                                          <  2:925039/150‑1 (MQ=255)
           gATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCa                                                                                                                                          >  1:856852/1‑150 (MQ=255)
                cGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggca                                                                                                                                     <  1:1218112/150‑4 (MQ=255)
                 gTTCCATCCAGAGTGCATCGTCCTGTACCGGATATGTTTCTCAGTTACGCTCACCACGATAGAACGCAATCAAAGCCGCTATTCCGAAAGTAAGCCGCGCCTGAAGTATGACGTTATCCATCTGCCATGCCAGAAGCTC‑TGGCAGGAGGc                                                                                                                                    >  1:290059/1‑150 (MQ=255)
                    ccagccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                                                                 <  1:505900/150‑1 (MQ=255)
                    ccagccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                                                                 >  1:1670400/1‑150 (MQ=255)
                      agccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCATGATGCGGGTg                                                                                                                               >  2:1676951/1‑150 (MQ=255)
                       gccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTGACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGc                                                                                                                              <  1:1291768/150‑1 (MQ=255)
                              tgCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACtt                                                                                                                       >  1:1790918/1‑150 (MQ=255)
                               gCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTg                                                                                                                      >  2:15826/1‑150 (MQ=255)
                                cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                                                                     >  2:1591295/1‑150 (MQ=255)
                                 ctCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                                                                    <  2:437566/149‑1 (MQ=255)
                                  tCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTc                                                                                                                   <  1:1671756/150‑1 (MQ=255)
                                   cGTCCTGTACCGGATATGATTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCa                                                                                                                  <  2:632514/150‑1 (MQ=255)
                                    gTCCTGTACCGGATATGTTTCGCCGTTTCGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGGGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGGGCGGAACTTGGTCAt                                                                                                                 <  1:541874/150‑1 (MQ=255)
                                    gTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCAt                                                                                                                 >  1:662180/1‑150 (MQ=255)
                                       cTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                                                                              >  2:213948/1‑150 (MQ=255)
                                       cTGTACCGGATATGTTTCGCCGTTACGCCCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGGCTTCTGCCCTGCCGGCAGCTG‑TGCCAGGATGCGGGTGCGGAACTTGGTCAAAcc                                                                                                              >  2:751927/1‑150 (MQ=255)
                                       cTGTACCGGATATGTTTCGCCGGTACGCTCACCACGATAGAACGCAATCAATGCCGCTCATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                                                                              <  1:1132991/150‑1 (MQ=255)
                                       cTGTACCGGATAAGTTTCGCCGTTACGCGCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                                                                              <  1:301216/150‑1 (MQ=255)
                                       cTGGACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCACTGCCGCTAATGCGACAGTACGGCGCGCCGGAAGGGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                                                                              <  2:20410/150‑1 (MQ=255)
                                         gTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATGCCGt                                                                                                            >  2:652976/1‑150 (MQ=255)
                                          tACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGtt                                                                                                           >  2:977325/1‑150 (MQ=255)
                                           aCCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTg                                                                                                          >  2:553392/1‑150 (MQ=255)
                                             cGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAg                                                                                                        >  1:1430380/1‑150 (MQ=255)
                                              ggATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgc                                                                                                       >  1:1860501/1‑150 (MQ=255)
                                                 tatGTTTCGCCGTTACGCTCACCACGATAGAACGCACTCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTCCCTGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                                                                    >  1:82586/1‑150 (MQ=255)
                                                     tttCGCCGGTACGCGCACCACGATAGCACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTC‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                                                                <  2:316448/150‑1 (MQ=255)
                                                      ttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTTCGGAACTTGGTCATACCGTTGAGCGCGATAGa                                                                                               >  1:519762/1‑150 (MQ=255)
                                                      ttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGa                                                                                               >  1:1060599/1‑150 (MQ=255)
                                                      ttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGa                                                                                               >  1:433078/1‑150 (MQ=255)
                                                         gCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcag                                                                                            <  2:1355232/150‑1 (MQ=255)
                                                         gCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcag                                                                                            >  1:428851/1‑150 (MQ=255)
                                                           cGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCCGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagca                                                                                          >  1:1430612/1‑150 (MQ=255)
                                                            gTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagcag                                                                                         <  2:1860501/150‑1 (MQ=255)
                                                              tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCt                                                                                       >  2:111102/1‑150 (MQ=255)
                                                                cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGa                                                                                     >  2:653644/1‑150 (MQ=255)
                                                                  cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                                                                   >  2:753475/1‑150 (MQ=255)
                                                                       cacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                                                                              >  2:441772/1‑150 (MQ=255)
                                                                        acGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAAt                                                                             >  1:207754/1‑150 (MQ=255)
                                                                         cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATg                                                                            <  1:896228/150‑1 (MQ=255)
                                                                              gaaaccaaacaaaGCCGCTAATCCGAAAGTAAGGAGCCCCTGAAGTGAGCCGTTTTACTTATGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                                                                       >  2:1762172/6‑150 (MQ=255)
                                                                                    aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCTATAGACAGCAGCTGATGCNTAATGTACGGGTTACg                                                                 <  1:1517877/150‑1 (MQ=255)
                                                                                    aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                                                                 <  2:53085/150‑1 (MQ=255)
                                                                                    aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                                                                 >  2:1467693/1‑150 (MQ=255)
                                                                                      tcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGa                                                               <  2:688605/150‑1 (MQ=255)
                                                                                      tcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGa                                                               <  2:593679/150‑1 (MQ=255)
                                                                                        aatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGCCCGCAGCTGATGCTTAATGTACGGGTTACGGaaa                                                             >  2:899486/1‑150 (MQ=255)
                                                                                        aatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaaa                                                             <  2:1722553/150‑1 (MQ=255)
                                                                                           gCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGc                                                          >  2:604487/1‑150 (MQ=255)
                                                                                            ccGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc                                                         <  1:373134/150‑1 (MQ=255)
                                                                                                  aTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGa                                                   >  1:666480/1‑150 (MQ=255)
                                                                                                   tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc                                                  >  1:483462/1‑150 (MQ=255)
                                                                                                    gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGGGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTCAGGAATCGCCCGGTGTCt                                                 >  2:220772/1‑150 (MQ=255)
                                                                                                    gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACt                                                 >  2:859287/1‑150 (MQ=255)
                                                                                                     cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTg                                                >  1:1059929/1‑150 (MQ=255)
                                                                                                     cGAAAGGAACGCGCGCCGGACGTGGCCCGTTTGCCCTCTGCCCTGCCAGCCGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCTTACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTGACGGAAACGCCCGGTGACTg                                                <  2:297053/150‑1 (MQ=255)
                                                                                                         aGTAAGGCGCGCCGGAAGTGGGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACt                                            <  1:1124090/150‑1 (MQ=255)
                                                                                                            aaGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGAAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTGGAGCGCGATAGACAGCAGCTGATGCTTAATGGACGGGTTACGGAATCGCCCGTTGACTGCACTGGc                                         >  1:274619/1‑150 (MQ=255)
                                                                                                            aaGGCGCGCCGGAAGGGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGc                                         <  2:1080380/150‑1 (MQ=255)
                                                                                                              ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGa                                       >  1:146718/1‑150 (MQ=255)
                                                                                                              ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTAAGGGTTACGGAAACGCCCGGTGACTGCACTGGCGa                                       >  2:608995/1‑150 (MQ=255)
                                                                                                               gcgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGaa                                      >  1:1831290/1‑150 (MQ=255)
                                                                                                                cgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGaaa                                     <  1:486434/150‑1 (MQ=255)
                                                                                                                cgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGaaa                                     >  1:418574/1‑150 (MQ=255)
                                                                                                                 gcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAg                                    <  1:1311481/150‑1 (MQ=255)
                                                                                                                 gcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAg                                    <  2:1860410/150‑1 (MQ=255)
                                                                                                                   gcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGTACTTGGTAATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGAt                                  <  1:1676951/150‑1 (MQ=255)
                                                                                                                    ccGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGAtt                                 >  1:614437/1‑150 (MQ=255)
                                                                                                                      ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTcc                               <  1:1694435/150‑1 (MQ=255)
                                                                                                                         aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTTAGCTCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGt                            <  1:1699915/150‑1 (MQ=255)
                                                                                                                              gCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtcatc                       >  2:503355/1‑150 (MQ=255)
                                                                                                                              gCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtcatc                       >  2:72079/1‑150 (MQ=255)
                                                                                                                               ccGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtcatca                      <  2:428851/150‑1 (MQ=255)
                                                                                                                                   ttGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAg                  <  1:652976/150‑1 (MQ=255)
                                                                                                                                       cccgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaa              >  1:1272665/8‑150 (MQ=255)
                                                                                                                                       cccgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaa              >  1:1763403/8‑150 (MQ=255)
                                                                                                                                          gcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATc           >  2:31192/5‑150 (MQ=255)
                                                                                                                                           tgccctgccAGCAGCTG‑TGGCAGGATGCGGGGGCGGAACTTGGTCATACCGTTCAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGGTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc          <  2:964350/150‑1 (MQ=255)
                                                                                                                                           catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc          >  1:866871/4‑150 (MQ=255)
                                                                                                                                            gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCa         >  2:1671110/1‑150 (MQ=255)
                                                                                                                                             ccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg        >  1:204779/1‑150 (MQ=255)
                                                                                                                                             ccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg        <  2:758295/150‑1 (MQ=255)
                                                                                                                                             ccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg        <  1:796626/150‑1 (MQ=255)
                                                                                                                                             ccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg        >  1:1046787/1‑150 (MQ=255)
                                                                                                                                              cctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGt       >  1:621087/1‑150 (MQ=255)
                                                                                                                                               ctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGCTTCCAGTTCATCACGAGGCAAATCCAGTa      >  2:1104883/1‑150 (MQ=255)
                                                                                                                                                tgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAc     <  1:987039/150‑1 (MQ=255)
                                                                                                                                                   cagcagCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCgg  >  2:484469/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGAGTTGCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGG  >  NC_000913/1609406‑1609702

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