Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,198,529 A→G T272A (ACA→GCA)  hemE → uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,198,5290AG98.2% 388.1 / ‑2.7 110T272A (ACA→GCA) hemEuroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (57/51);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 113 total aligned reads displayed.

ACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCA  >  NC_000913/4198392‑4198665
                                                                                                                                         |                                                                                                                                        
aCATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCa                                                                                                                                        >  1:955433/1‑140 (MQ=255)
    gCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACgg                                                                                                                                    >  1:85372/1‑140 (MQ=255)
    gCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACgg                                                                                                                                    >  2:771544/1‑140 (MQ=255)
    gCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACgg                                                                                                                                    >  2:345614/1‑140 (MQ=255)
     cATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGCCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGa                                                                                                                                   <  1:205229/140‑1 (MQ=255)
      aTAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCa                                                                                                                                        >  1:348058/1‑134 (MQ=255)
      aTAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGAc                                                                                                                                  >  1:1037653/1‑140 (MQ=255)
      aTAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGAc                                                                                                                                  >  1:599169/1‑140 (MQ=255)
      aTAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGAc                                                                                                                                  >  2:408940/1‑140 (MQ=255)
      aTAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGAc                                                                                                                                  >  2:789519/1‑140 (MQ=255)
      aTAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGAc                                                                                                                                  >  2:891369/1‑140 (MQ=255)
        aaaaTTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACAt                                                                                                                                >  2:408449/1‑140 (MQ=255)
        aaaaTTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACAt                                                                                                                                >  2:373668/1‑140 (MQ=255)
        aaaaTTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACAt                                                                                                                                >  1:970958/1‑140 (MQ=255)
        aaaaTTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACAt                                                                                                                                >  1:846433/1‑140 (MQ=255)
        aaaaTTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACAt                                                                                                                                >  2:589444/1‑140 (MQ=255)
                           tGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCgt                                                                                                             <  2:85372/140‑1 (MQ=255)
                           tGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCgt                                                                                                             <  1:420338/140‑1 (MQ=255)
                            gCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCgtg                                                                                                            >  2:118401/1‑140 (MQ=255)
                            gCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCgt                                                                                                             >  1:1183684/1‑139 (MQ=255)
                            gCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCgt                                                                                                             <  2:1183684/139‑1 (MQ=255)
                             cGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCgtgt                                                                                                           >  2:1009800/1‑140 (MQ=255)
                              gTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGtgtg                                                                                                          >  1:758136/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAGAGGCGGTGGACAGTGTCTGGGAGCCATGGCAGAAACCGGTTGCGATGCGCGGAGCCTCGACTGGGCAACTGACATCGCCGATGGGCGCCGCGGTGTgg                                                                                                         >  2:725549/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCTGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGACTCGACTGGGCAACGGAAATCGCCGATGCGCGCCGCCGTGTgg                                                                                                         >  2:903085/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGcgcc                                                                                                                  <  2:26791/131‑1 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGcgcc                                                                                                                  >  1:26791/1‑131 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTgg                                                                                                         >  2:387681/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTgg                                                                                                         >  1:888569/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTgg                                                                                                         >  2:670888/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTgg                                                                                                         >  2:565131/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTgg                                                                                                         >  2:536571/1‑140 (MQ=255)
                               tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTgg                                                                                                         >  2:457066/1‑140 (MQ=255)
                                gAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGccgccg                                                                                                              >  1:368540/1‑134 (MQ=255)
                                gAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGccgccg                                                                                                              <  2:368540/134‑1 (MQ=255)
                                gAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTggg                                                                                                        <  1:711165/140‑1 (MQ=255)
                                gAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTggg                                                                                                        <  2:758136/140‑1 (MQ=255)
                                gAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTggg                                                                                                        <  1:408449/140‑1 (MQ=255)
                                      gacgGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATaa                                                                                                  <  1:891369/140‑1 (MQ=255)
                                           tCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTcg                                                                                             <  2:599169/140‑1 (MQ=255)
                                           tCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTcg                                                                                             <  1:789519/140‑1 (MQ=255)
                                             gccgcGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTcgcg                                                                                           >  1:1154132/1‑140 (MQ=255)
                                             gccgcGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTcgcg                                                                                           <  2:1154132/140‑1 (MQ=255)
                                                gcgTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTg                                                                                        <  2:924337/140‑1 (MQ=255)
                                                gcgTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTg                                                                                        <  1:565131/140‑1 (MQ=255)
                                                gcgTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTg                                                                                        <  1:761029/140‑1 (MQ=255)
                                                 cgTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGc                                                                                       <  1:250315/140‑1 (MQ=255)
                                                  gTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCa                                                                                      <  2:888569/140‑1 (MQ=255)
                                                  gTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCa                                                                                      <  1:457066/140‑1 (MQ=255)
                                                     ccGGTCACGCTGTTTTCCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAggg                                                                                   <  1:1042028/140‑1 (MQ=255)
                                                     ccGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAggg                                                                                   <  1:807453/140‑1 (MQ=255)
                                                       ggTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTa                                                                                 >  2:79327/1‑140 (MQ=255)
                                                        gTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTaa                                                                                >  1:1222159/1‑140 (MQ=255)
                                                         tCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAat                                                                               >  1:1134043/1‑140 (MQ=255)
                                                             gCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCTCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATgg                                                                           >  2:436918/1‑140 (MQ=255)
                                                                gTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATc                                                                        <  2:1037653/140‑1 (MQ=255)
                                                                   tACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGt                                                                     <  2:1134043/140‑1 (MQ=255)
                                                                    aCCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTc                                                                    <  1:536571/140‑1 (MQ=255)
                                                                         aGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGc                                                               >  2:398109/1‑140 (MQ=255)
                                                                           gcggcgGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTTGGCCTCGACTGGGCCACGGACATCGCCGATGCGCGCCGCCGCGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGGTg                                                             >  2:622856/1‑140 (MQ=255)
                                                                           gcggcgGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTg                                                             >  1:1223233/1‑140 (MQ=255)
                                                                           gcggcgGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTg                                                             >  2:721772/1‑140 (MQ=255)
                                                                           gcggcgGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTg                                                             >  2:468743/1‑140 (MQ=255)
                                                                           gcggcgGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTg                                                             >  2:474757/1‑140 (MQ=255)
                                                                           gcggcgGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTg                                                             >  2:81723/1‑140 (MQ=255)
                                                                           gcggcgGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTg                                                             >  1:951924/1‑140 (MQ=255)
                                                                             ggcggACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTa                                                           >  1:1135179/1‑140 (MQ=255)
                                                                               cggACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTAcg                                                         <  1:408940/140‑1 (MQ=255)
                                                                                     gTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGcgccgc                                                   <  1:1009800/140‑1 (MQ=255)
                                                                                      tGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCgccgcc                                                  <  2:951924/140‑1 (MQ=255)
                                                                                              aGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGcc                                              <  1:245243/136‑1 (MQ=255)
                                                                                              aGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGcc                                              >  2:245243/1‑136 (MQ=255)
                                                                                              aGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCa                                          <  1:771544/140‑1 (MQ=255)
                                                                                              aGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCa                                          <  1:539335/140‑1 (MQ=255)
                                                                                              aGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCa                                          <  2:1222159/140‑1 (MQ=255)
                                                                                              aGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCa                                          >  2:515188/1‑140 (MQ=255)
                                                                                              aGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGAACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGGACGCGCCGCCTGCCCGCa                                          >  2:1186050/1‑140 (MQ=255)
                                                                                                ccATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCAtt                                        >  1:158280/1‑140 (MQ=255)
                                                                                                ccATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCAtt                                        >  2:658087/1‑140 (MQ=255)
                                                                                                  aTGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGcc                                              <  2:639673/132‑1 (MQ=255)
                                                                                                  aTGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTga                                      >  2:1088569/1‑140 (MQ=255)
                                                                                                       aGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTgaagaag                                 <  2:1223233/140‑1 (MQ=255)
                                                                                                           aCCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCa                                          <  2:727172/127‑1 (MQ=255)
                                                                                                           aCCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCa                                          >  1:727172/1‑127 (MQ=255)
                                                                                                           aCCGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGt                             <  1:398109/140‑1 (MQ=255)
                                                                                                           aCCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGt                             <  1:589444/140‑1 (MQ=255)
                                                                                                            ccGGTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTa                            <  2:331803/140‑1 (MQ=255)
                                                                                                              ggTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGc                          <  1:314206/140‑1 (MQ=255)
                                                                                                              ggTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGc                          <  1:670888/140‑1 (MQ=255)
                                                                                                               gTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGt                             >  1:419613/1‑136 (MQ=255)
                                                                                                               gTTGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCg                         >  1:414721/1‑140 (MQ=255)
                                                                                                                ttGCGATGCGCTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGa                        <  1:81723/140‑1 (MQ=255)
                                                                                                                        cgcTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACtt                <  1:468743/140‑1 (MQ=255)
                                                                                                                        cgcTGGGCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACtt                >  2:904145/1‑140 (MQ=255)
                                                                                                                             ggCCTCGACTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAgg           >  2:666465/1‑140 (MQ=255)
                                                                                                                                    aCTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGt    <  1:1088569/140‑1 (MQ=255)
                                                                                                                                    aCTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGt    <  1:904145/140‑1 (MQ=255)
                                                                                                                                     cTGGGCAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTc   <  1:903085/140‑1 (MQ=255)
                                                                                                                                     cTGGGCAACGGACATCGCCGATGAGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTc   <  1:725549/140‑1 (MQ=255)
                                                                                                                                      tGGGCAACAGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCa  <  2:491936/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
ACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGGTAATATGGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCA  >  NC_000913/4198392‑4198665

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