Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,460,381 G→A T30M (ACG→ATG)  fadI ← beta‑ketoacyl‑CoA thiolase, anaerobic, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,460,3810GA99.3% 556.8 / ‑4.9 146T30M (ACG→ATG) fadIbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (64/81);  total (64/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 146 total aligned reads displayed.

GCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCGTCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATT  >  NC_000913/2460238‑2460528
                                                                                                                                               |                                                                                                                                                   
gcgcAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTg                                                                                                                                              <  1:1181668/151‑1 (MQ=255)
 cgcAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTg                                                                                                                                              >  1:861049/1‑150 (MQ=255)
 cgcAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGa                                                                                                                                             <  4:1619275/151‑1 (MQ=255)
     aTGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCCGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCCTGAAAAACCATCCCCTGACggg                                                                                                                                         >  1:1643869/1‑151 (MQ=255)
     aTGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACggg                                                                                                                                         <  1:1486224/151‑1 (MQ=255)
     aTGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACggg                                                                                                                                         >  4:716771/1‑151 (MQ=255)
       gTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCCTCCCTAAATCAACCGCGGGAAGGCCATTAAAAACCAGCGCCTGACGGGc                                                                                                                                        >  1:893489/1‑150 (MQ=255)
        ttGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAACAGCCATCGCCTGACGGGc                                                                                                                                        >  1:658477/1‑149 (MQ=255)
        ttGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCa                                                                                                                                       >  1:1403438/1‑150 (MQ=255)
         tGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCaa                                                                                                                                      <  1:1130658/150‑1 (MQ=255)
          ggggtCTTCAGGCATTTGTACGACCTGTCCAAAGACCAGTTGTTCAATCACTCCGGCGGGGATCCCGTTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAACCATCCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCaaaa                                                                                                                                    <  1:390848/151‑1 (MQ=255)
          gggggCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCaaaa                                                                                                                                    <  1:186199/151‑1 (MQ=255)
          gggggCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCaaa                                                                                                                                     >  1:956863/1‑150 (MQ=255)
              gCTTCAGGCATTTGTACGACCTGACCAAAGACCAGATGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGCTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGGCATGAAAAGCCATCGCCTGACGGGCaaaa                                                                                                                                    >  1:384727/1‑147 (MQ=255)
                 tCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGcgta                                                                                                                              <  4:1328126/150‑1 (MQ=255)
                     gCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTaa                                                                                                                         <  4:75444/151‑1 (MQ=255)
                        tttGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAAc                                                                                                                       <  3:536149/150‑1 (MQ=255)
                        tttGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAAc                                                                                                                       <  1:1558181/150‑1 (MQ=255)
                           gTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCt                                                                                                                   <  1:470106/151‑1 (MQ=255)
                             aCGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTTGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTa                                                                                                                  >  1:665234/1‑150 (MQ=255)
                             aCGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTa                                                                                                                  >  4:1367307/1‑150 (MQ=255)
                             aCGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTa                                                                                                                  <  1:116825/150‑1 (MQ=255)
                              cGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAAc                                                                                                                <  1:1446005/151‑1 (MQ=255)
                                 ccTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACaa                                                                                                              >  4:728168/1‑150 (MQ=255)
                                  cTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATg                                                                                                            >  4:1165784/1‑151 (MQ=255)
                                    gACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATg                                                                                                            <  4:227073/149‑1 (MQ=255)
                                    gACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATg                                                                                                            <  1:824899/149‑1 (MQ=255)
                                     aCCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCg                                                                                                         <  3:487358/151‑1 (MQ=255)
                                      ccAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGa                                                                                                        <  1:1199327/151‑1 (MQ=255)
                                      ccAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGa                                                                                                        <  4:769034/151‑1 (MQ=255)
                                        aaaGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATa                                                                                                      <  4:1506823/151‑1 (MQ=255)
                                         aaGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAAACGCTAACAATGGCGGTa                                                                                                      >  3:447052/1‑150 (MQ=255)
                                          aGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACg                                                                                                    <  1:226711/151‑1 (MQ=255)
                                              cAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCg                                                                                                <  1:1491220/151‑1 (MQ=255)
                                                gttgttCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGcc                                                                                              >  1:1335275/1‑151 (MQ=255)
                                                  tgttCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGccc                                                                                             <  3:344055/150‑1 (MQ=255)
                                                      cAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGc                                                                                         <  4:1460091/150‑1 (MQ=255)
                                                      cAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGc                                                                                         <  4:682390/150‑1 (MQ=255)
                                                        aTCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCggg                                                                                      >  4:1393820/1‑151 (MQ=255)
                                                          cACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGt                                                                                     <  4:583411/150‑1 (MQ=255)
                                                           aCTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTaa                                                                                   <  4:426518/151‑1 (MQ=255)
                                                            cTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTaa                                                                                   >  4:970301/1‑150 (MQ=255)
                                                             ttCGGCGGGGGTCTCGCTGCGGGCCAGCAGTTTGCCTACCACCATCTTTCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGGACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGGaa                                                                                   >  1:886123/1‑149 (MQ=255)
                                                               cggcggGGAACTCCCTGGGTGCCAGCAGGTCGCCTACCACCCTCTTCCCTAAATCCACCGCGGGAATGGCATGAAAAACCATCACCTGACGGGCAAAAGGCGTACGTAAAACGCTAACAATGGCGATACGATCGCCCCGGCGGGGaa                                                                                   >  1:1250993/1‑147 (MQ=255)
                                                                  cggGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCgg                                                                            <  1:1091853/151‑1 (MQ=255)
                                                                    gggATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTa                                                                          <  4:1535468/151‑1 (MQ=255)
                                                                    gggATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTa                                                                          <  4:1254483/151‑1 (MQ=255)
                                                                     ggATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTa                                                                          >  2:29083/1‑150 (MQ=255)
                                                                        tctcGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAAc                                                                      >  1:812410/1‑151 (MQ=255)
                                                                         ctcGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAAcc                                                                     >  4:95618/1‑151 (MQ=255)
                                                                           cGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTg                                                                   <  1:1046029/151‑1 (MQ=255)
                                                                           cGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTg                                                                   <  1:284939/151‑1 (MQ=255)
                                                                            gCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGa                                                                  <  4:76762/151‑1 (MQ=255)
                                                                                cGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGAccc                                                               >  1:1481708/1‑150 (MQ=255)
                                                                                cGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGAccc                                                               >  1:466724/1‑150 (MQ=255)
                                                                                      agcagTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAacac                                                        <  4:469237/151‑1 (MQ=255)
                                                                                         agTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTcc                                                     >  4:1760960/1‑151 (MQ=255)
                                                                                         agTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCCAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTcc                                                     >  1:1546073/1‑151 (MQ=255)
                                                                                             cGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCt                                                  <  2:408418/150‑1 (MQ=255)
                                                                                              gCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGt                                                <  3:25319/151‑1 (MQ=255)
                                                                                              gCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGt                                                <  1:626589/151‑1 (MQ=255)
                                                                                               ccTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTa                                               <  2:538772/151‑1 (MQ=255)
                                                                                                 tACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGt                                                >  1:53390/1‑148 (MQ=255)
                                                                                                 tACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTaa                                              >  4:1538848/1‑150 (MQ=255)
                                                                                                  accaccATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTaaa                                             >  1:1791632/1‑150 (MQ=255)
                                                                                                  accaccATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCGGACCCATAACACTCCCCTGTaaa                                             <  1:1255812/150‑1 (MQ=255)
                                                                                                   ccaccaTCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTaaaaa                                           <  4:934314/151‑1 (MQ=255)
                                                                                                   ccaccaTCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTaaaa                                            >  1:1192925/1‑150 (MQ=255)
                                                                                                    caccaTCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCACAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTaaa                                             <  4:896136/148‑1 (MQ=255)
                                                                                                      ccaTCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCGTCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAaaat                                        <  2:118779/151‑1 (MQ=255)
                                                                                                      ccaTCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAaaat                                        >  4:1179740/1‑151 (MQ=255)
                                                                                                       caTCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAaaat                                        >  1:779067/1‑150 (MQ=255)
                                                                                                             cccTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAataaataa                                  >  4:1112998/1‑150 (MQ=255)
                                                                                                              ccTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAg                                <  2:264784/151‑1 (MQ=255)
                                                                                                              ccTAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAg                                <  4:1729537/151‑1 (MQ=255)
                                                                                                                tAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTg                              >  3:418821/1‑151 (MQ=255)
                                                                                                                tAAATCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTg                              >  1:1637200/1‑151 (MQ=255)
                                                                                                                 aaaTCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTgg                             <  3:488174/151‑1 (MQ=255)
                                                                                                                   aTCAACCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTc                           >  2:38356/1‑151 (MQ=255)
                                                                                                                      aaCCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTGCGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGa                        <  1:1869803/151‑1 (MQ=255)
                                                                                                                      aaCCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTg                         >  3:13514/1‑150 (MQ=255)
                                                                                                                      aaCCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGa                        <  2:513814/151‑1 (MQ=255)
                                                                                                                       aCCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGAc                       >  1:177138/1‑151 (MQ=255)
                                                                                                                       aCCGCGGGAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGAc                       >  4:493630/1‑151 (MQ=255)
                                                                                                                            gggAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGAt                  >  4:1917883/1‑151 (MQ=255)
                                                                                                                            gggAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGAt                  <  3:265848/151‑1 (MQ=255)
                                                                                                                            gggAATGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGAt                  >  4:322477/1‑151 (MQ=255)
                                                                                                                                aTGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATa              <  2:71974/151‑1 (MQ=255)
                                                                                                                                 tGCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCAt               >  4:1637281/1‑149 (MQ=255)
                                                                                                                                  gCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAg             <  1:750828/150‑1 (MQ=255)
                                                                                                                                  gCCATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGt            <  4:1075119/151‑1 (MQ=255)
                                                                                                                                   ccATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTc           >  3:328350/1‑151 (MQ=255)
                                                                                                                                   ccATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTc           >  1:1928957/1‑151 (MQ=255)
                                                                                                                                   ccATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTc           <  4:557603/151‑1 (MQ=255)
                                                                                                                                   ccATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTc           <  4:280660/151‑1 (MQ=255)
                                                                                                                                   ccATGAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCGGACCTGATCATAGTc           <  3:201078/151‑1 (MQ=255)
                                                                                                                                       gAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTaa       <  4:1664401/151‑1 (MQ=255)
                                                                                                                                       gAAAAGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTa        <  1:1334944/150‑1 (MQ=255)
                                                                                                                                         aaaGCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAAcc     <  4:237468/151‑1 (MQ=255)
                                                                                                                                            gCCATCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCAtt  <  4:741532/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCGCAATGTTGGGGGCTTCAGGCATTTGTACGACCTGACCAAAGACCAGTTGTTCAATCACTTCGGCGGGGATCTCGCTGCGTGCCAGCAGTTCGCCTACCACCATCTTCCCTAAATCAACCGCGGGAATGCCATGAAAAGCCGTCGCCTGACGGGCAAAAGGCGTACGTAAACCGCTAACAATGGCGATACGATCGCCCTGGCGGGTAACCAGCGGTAAAACCTGACCCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATT  >  NC_000913/2460238‑2460528

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