Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,555 Δ1 bp 5.4% 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.4% 395.7 / 15.6 112coding (1126/1452 nt)uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (47/59);  new base . (2/4);  total (49/63)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

Only 100 of 112 total aligned reads displayed.

GCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAC  >  NC_000913/1609414‑1609699
                                                                                                                                             |                                                                                                                                                 
gCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtgt                                                                                                                                           <  4:413164/149‑1 (MQ=255)
  tGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtgtg                                                                                                                                          <  2:406451/148‑1 (MQ=255)
      aaCGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAgg                                                                                                                                     >  4:385299/1‑149 (MQ=255)
      aaCGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAg                                                                                                                                      <  2:736756/148‑1 (MQ=255)
      aaCGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgg                                                                                                                                     <  3:336530/149‑1 (MQ=255)
           tCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGc                                                                                                                                 >  2:642379/1‑148 (MQ=255)
            ccagccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCg                                                                                                                                <  2:483516/148‑1 (MQ=255)
            ccagccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgggcaga                                                                                                                                <  1:218467/148‑6 (MQ=255)
             cagccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGTCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCggg                                                                                                                              >  4:441565/1‑149 (MQ=255)
              agccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGt                                                                                                                             >  3:808276/1‑149 (MQ=255)
              agccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGt                                                                                                                             >  4:82611/1‑149 (MQ=255)
              agccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGt                                                                                                                             >  4:155681/1‑149 (MQ=255)
               gccagTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTg                                                                                                                            <  4:827908/149‑1 (MQ=255)
                   gtgtGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCgg                                                                                                                         <  1:416435/148‑1 (MQ=255)
                     gtgCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAAc                                                                                                                      >  3:860132/1‑149 (MQ=255)
                        cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACtt                                                                                                                    >  2:55108/1‑148 (MQ=255)
                        cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACtt                                                                                                                    <  1:624010/148‑1 (MQ=255)
                        cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTg                                                                                                                   <  3:385299/149‑1 (MQ=255)
                        cATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTg                                                                                                                   >  3:449064/1‑149 (MQ=255)
                         aTCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTgg                                                                                                                  <  4:256322/149‑1 (MQ=255)
                         aTCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTgg                                                                                                                  >  4:900094/1‑149 (MQ=255)
                         aTCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTgg                                                                                                                  >  3:256322/1‑149 (MQ=255)
                           cGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGt                                                                                                                 >  2:741008/1‑148 (MQ=255)
                           cGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCa                                                                                                                                       <  2:606763/126‑1 (MQ=255)
                           cGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCa                                                                                                                                       >  1:606763/1‑126 (MQ=255)
                              ccTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCAt                                                                                                              >  2:169701/1‑148 (MQ=255)
                                tGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                                                                           <  3:900094/149‑1 (MQ=255)
                                tGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATAc                                                                                                            >  2:726631/1‑148 (MQ=255)
                                   aCCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGtt                                                                                                        <  3:94751/149‑1 (MQ=255)
                                   aCCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGt                                                                                                         <  2:251713/148‑1 (MQ=255)
                                     cGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCGCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTg                                                                                                       >  2:549719/1‑148 (MQ=255)
                                     cGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGa                                                                                                      <  3:670237/149‑1 (MQ=255)
                                      ggATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAg                                                                                                     <  3:84030/149‑1 (MQ=255)
                                      ggATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAg                                                                                                     >  4:84030/1‑149 (MQ=255)
                                      ggATATGTTTCGCCGTTACGCTCACCACCATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAg                                                                                                     >  3:60060/1‑149 (MQ=255)
                                         tatGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcg                                                                                                   >  1:830157/1‑148 (MQ=255)
                                         tatGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcg                                                                                                   <  1:140092/148‑1 (MQ=255)
                                          atGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcgc                                                                                                  <  1:169701/148‑1 (MQ=255)
                                          atGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcgc                                                                                                  <  1:473937/148‑1 (MQ=255)
                                           tGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                                                                 <  2:754755/148‑1 (MQ=255)
                                             tttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                                                               <  2:124515/148‑1 (MQ=255)
                                              ttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATa                                                                                              >  2:723011/1‑148 (MQ=255)
                                               tCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                                                             <  2:321669/148‑1 (MQ=255)
                                                 gCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                                                                >  4:201137/1‑143 (MQ=255)
                                                 gCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                                                                <  3:201137/143‑1 (MQ=255)
                                                  ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcag                                                                                         >  3:378875/1‑149 (MQ=255)
                                                  ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                                                                          <  1:409154/148‑1 (MQ=255)
                                                     ttACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagca                                                                                       <  2:892430/148‑1 (MQ=255)
                                                      tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGGTAGACagcagc                                                                                     <  3:514998/149‑1 (MQ=255)
                                                      tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACagcagc                                                                                     <  3:82611/149‑1 (MQ=255)
                                                           tCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                                                                >  4:871667/1‑149 (MQ=255)
                                                           tCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                                                                >  3:71790/1‑149 (MQ=255)
                                                            caccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                                                                <  1:624559/148‑1 (MQ=255)
                                                               cacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTa                                                                            >  3:703727/1‑149 (MQ=255)
                                                                acGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTa                                                                            <  2:867945/148‑1 (MQ=255)
                                                                 cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACagcagc                                                                                     >  1:392850/1‑138 (MQ=255)
                                                                 cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACagcagc                                                                                     <  2:392850/138‑1 (MQ=255)
                                                                 cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                                                                           <  1:772812/148‑1 (MQ=255)
                                                                  gATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTg                                                                                                                   >  4:49995/1‑107 (MQ=255)
                                                                  gATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTg                                                                                                                   <  3:49995/107‑1 (MQ=255)
                                                                   aTAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATg                                                                         >  2:72296/1‑148 (MQ=255)
                                                                      gatctcaatcaatGCCGCTAATGCGAAAGTAAGGCCCGCCGGAAGTGTGCCGTTTGCCTTCTTCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCAGTCCGTTGAGCGCGACAGACAGCAGCTGATGCATAATGTAc                                                                      <  1:201229/143‑1 (MQ=255)
                                                                       aaCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                                                                    <  4:449064/149‑1 (MQ=255)
                                                                        aCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                                                                    <  4:881292/148‑1 (MQ=255)
                                                                        aCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                                                                    >  3:881292/1‑148 (MQ=255)
                                                                            aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTc                                                                                                                >  3:549639/1‑100 (MQ=255)
                                                                            aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTc                                                                                                                <  4:549639/100‑1 (MQ=255)
                                                                            aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTAAATGTACGGGTTAc                                                               >  3:359076/1‑149 (MQ=255)
                                                                                  tGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc                                                         <  4:808276/149‑1 (MQ=255)
                                                                                   gCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg                                                        >  3:600471/1‑149 (MQ=255)
                                                                                     cGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCATGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGc                                                       <  2:245733/148‑1 (MQ=255)
                                                                                      gCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCTATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGccc                                                     >  3:353371/1‑149 (MQ=255)
                                                                                           tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTg                                                 <  1:72296/148‑1 (MQ=255)
                                                                                            gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc                                               <  3:416667/149‑1 (MQ=255)
                                                                                             cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc                                               <  1:494171/148‑1 (MQ=255)
                                                                                              gAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTg                                             >  4:557506/1‑149 (MQ=255)
                                                                                               aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTg                                             >  2:560946/1‑148 (MQ=255)
                                                                                                 aGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCa                                           <  1:247583/148‑1 (MQ=255)
                                                                                                    aaGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTg                                        <  2:316475/148‑1 (MQ=255)
                                                                                                         gcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGaaa                                  <  3:539964/149‑1 (MQ=255)
                                                                                                         gcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGaa                                   <  2:241143/148‑1 (MQ=255)
                                                                                                             cGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGAtt                              >  4:517811/1‑149 (MQ=255)
                                                                                                              ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGAtt                              <  2:326667/148‑1 (MQ=255)
                                                                                                               gAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTc                             >  1:791583/1‑148 (MQ=255)
                                                                                                                aaGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCa                           <  4:860132/149‑1 (MQ=255)
                                                                                                                aaGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTAAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCa                           <  4:359076/149‑1 (MQ=255)
                                                                                                                   tgtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGtt                        <  4:600471/149‑1 (MQ=255)
                                                                                                                      gCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtca                      >  2:832777/1‑148 (MQ=255)
                                                                                                                        cGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTtcatc                    >  2:790143/1‑148 (MQ=255)
                                                                                                                           ttGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACg                 <  2:163705/148‑1 (MQ=255)
                                                                                                                           ttGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGa                <  3:745125/149‑1 (MQ=255)
                                                                                                                              ccTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGc             <  3:339310/149‑1 (MQ=255)
                                                                                                                                ttCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCa            >  1:804997/1‑148 (MQ=255)
                                                                                                                                ttCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACAGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaa           >  4:577821/1‑149 (MQ=255)
                                                                                                                                 tCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaa           >  1:845787/1‑148 (MQ=255)
                                                                                                                                   tgccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAAt         <  1:560946/148‑1 (MQ=255)
                                                                                                                                    atcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc       >  3:507828/3‑149 (MQ=255)
                                                                                                                                     ccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc       <  2:347090/148‑1 (MQ=255)
                                                                                                                                        tgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGt    >  1:278032/1‑148 (MQ=255)
                                                                                                                                          ccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAc  <  1:723011/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GCTGATAACGTTCCAGCCAGTGTGCATCGTCCTGTACCGGATATGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAC  >  NC_000913/1609414‑1609699

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: