Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,008,626 G→A 5.1% V373I (GTT→ATT)  ygeX → 2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,008,6260GA 5.1% 367.7 / 3.4 114V373I (GTT→ATT) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (55/53);  new base A (4/2);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 119 total aligned reads displayed.

TGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAGAA  >  NC_000913/3008477‑3008774
                                                                                                                                                     |                                                                                                                                                    
tGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTgg                                                                                                                                                      <  2:212801/150‑1 (MQ=255)
tGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTgg                                                                                                                                                      <  1:210723/150‑1 (MQ=255)
 gCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGt                                                                                                                                                     <  1:361982/150‑1 (MQ=255)
 gCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGt                                                                                                                                                     <  2:490597/150‑1 (MQ=255)
     ggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGCTGCCGTGGTGCTGGTTATc                                                                                                                                                 >  1:215112/1‑150 (MQ=255)
     ggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATc                                                                                                                                                 >  1:194863/1‑150 (MQ=255)
     ggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATc                                                                                                                                                 <  2:1314599/150‑1 (MQ=255)
      gTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCa                                                                                                                                                <  1:466308/150‑1 (MQ=255)
        aaTCCGTACGGCAACGCCCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTTTCACCCGCACCGTCAAAGCCTGATGGAAAAACTGGCGCTGCACAAAGATGCCGTGGTGCTGGTTTTCAGc                                                                                                                                              >  2:338446/1‑150 (MQ=255)
        aaTCCGTACGGCAACGACCCGCGCATCTTCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATGAACCCCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGCACAAAGATGCCGTGGTGCTGGTTATCAga                                                                                                                                              >  1:1063010/1‑149 (MQ=255)
        aaTCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGc                                                                                                                                              >  2:1538649/1‑150 (MQ=255)
        aaTCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGc                                                                                                                                              <  2:231952/150‑1 (MQ=255)
        aaTCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGc                                                                                                                                              >  1:428640/1‑150 (MQ=255)
           ccGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACt                                                                                                                                           >  2:812980/1‑150 (MQ=255)
             gTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGa                                                                                                                                         >  2:129918/1‑150 (MQ=255)
             gTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGa                                                                                                                                         >  1:699988/1‑150 (MQ=255)
               aCGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGGCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGATAACCAGCAccacgg                                                                                                                                       <  2:90289/150‑6 (MQ=255)
                cGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAgg                                                                                                                                      >  2:1378542/1‑150 (MQ=255)
                cGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGGCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAgg                                                                                                                                      <  1:1474490/150‑1 (MQ=255)
                 ggCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGc                                                                                                                                     >  1:531888/1‑150 (MQ=255)
                  gCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCg                                                                                                                                    <  1:1220869/150‑1 (MQ=255)
                   cAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATAAGCACTGAAGGCGa                                                                                                                                   >  2:498580/1‑150 (MQ=255)
                           cgcgCATCATCTCCGGGGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAACGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACg                                                                                                                           <  2:1097948/150‑1 (MQ=255)
                            gcgcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGCCGACACCGACGt                                                                                                                          >  2:430051/1‑150 (MQ=255)
                             cgcATCATCTCCGGTGCCTCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTg                                                                                                                         <  2:115120/150‑1 (MQ=255)
                             cgcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTGCTCGCAGCGGTTCATAATCACCCGCAACGTCAAAGCCTGCTGGAAATACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTg                                                                                                                         <  2:119711/150‑1 (MQ=255)
                               catcatCTCCGGTGAATCCGGGGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGaa                                                                                                                       <  2:1015158/150‑1 (MQ=255)
                               catcatCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCACCACTGAAGGCGACACCGACGTGaa                                                                                                                       >  2:1440231/1‑150 (MQ=255)
                                 tcatcTCCGTTGAATCCCGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCTACGTCAAAGCCTGATGTAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGc                                                                                                                     <  2:671477/150‑1 (MQ=255)
                                 tcatcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGc                                                                                                                     >  1:1243343/1‑150 (MQ=255)
                                  catcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCa                                                                                                                    >  1:1141170/1‑150 (MQ=255)
                                    tctcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCATAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACt                                                                                                                  <  1:1670573/150‑1 (MQ=255)
                                         ggTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcgc                                                                                                             <  2:1666608/150‑1 (MQ=255)
                                         ggTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcgc                                                                                                             >  2:1317620/1‑150 (MQ=255)
                                            gAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCCGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGaa                                                                                                          >  1:505674/1‑150 (MQ=255)
                                            gAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGCCGACACCGACGTGAAGCACTACCGCGaa                                                                                                          >  1:378315/1‑150 (MQ=255)
                                             aaTCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAg                                                                                                         <  1:1084134/150‑1 (MQ=255)
                                               tCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGtt                                                                                                       <  2:580511/150‑1 (MQ=255)
                                                ccGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAACCTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTg                                                                                                      >  2:101281/1‑150 (MQ=255)
                                                    cgcTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTg                                                                                                  <  1:44895/150‑1 (MQ=255)
                                                    cgcTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTg                                                                                                  <  2:215112/150‑1 (MQ=255)
                                                    cgcTGTCGGTTTGGGCATTCTCGCAGCGGTGCATAATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGGGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTg                                                                                                  <  2:606400/150‑1 (MQ=255)
                                                      cTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTggg                                                                                                <  1:524927/150‑1 (MQ=255)
                                                       tGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGa                                                                                               <  1:1224370/150‑1 (MQ=255)
                                                       ggTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGa                                                                                               <  1:1683048/149‑1 (MQ=255)
                                                               tGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAac                                                                                       >  2:200549/1‑150 (MQ=255)
                                                                 ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                                                                                     >  1:1234580/1‑150 (MQ=255)
                                                                  gCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACg                                                                                    >  2:431587/1‑150 (MQ=255)
                                                                      tctcGCAGCGGAGCATTATCCCCCGCAACGTCAAACCCAGCTGGAAAAACTCGCGCGGAACAAAGATGCCGGGGTGCGGGTTATCAGCACTTAATGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGt                                                                                <  2:440078/150‑1 (MQ=255)
                                                                       ctcGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTa                                                                               >  2:736483/1‑150 (MQ=255)
                                                                       ctcGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTa                                                                               >  1:1553635/1‑150 (MQ=255)
                                                                         cGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGc                                                                             >  2:590427/1‑150 (MQ=255)
                                                                          gCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCa                                                                            >  2:750696/1‑150 (MQ=255)
                                                                            aGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCAcc                                                                          <  1:414387/150‑1 (MQ=255)
                                                                             gCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCt                                                                         <  2:1432102/150‑1 (MQ=255)
                                                                                 ttCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAAt                                                                     <  1:810074/150‑1 (MQ=255)
                                                                                    aTTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCAc                                                                  >  1:933377/1‑150 (MQ=255)
                                                                                     ttATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCGGAACAAAGATGCCGTGGGGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCAcc                                                                 <  2:1334375/150‑1 (MQ=255)
                                                                                      tATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCAGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGACGGCGCCACCGACGAGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACCCGCTGTAGCACCTTAATCAccc                                                                >  1:283187/1‑150 (MQ=255)
                                                                                          aCCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAGGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTAtt                                                            >  2:871073/1‑150 (MQ=255)
                                                                                             cGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGa                                                         <  1:878305/150‑1 (MQ=255)
                                                                                              gCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGaa                                                        <  1:15782/150‑1 (MQ=255)
                                                                                              gCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGaa                                                        <  1:200549/150‑1 (MQ=255)
                                                                                              gCAACGTCAAAGACTGATGGAAAAACTGGCTCTTAACAAAGATGCCGTGGTGCTGGTTATAATCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGaa                                                        <  2:1531011/150‑1 (MQ=255)
                                                                                               cAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAAc                                                       <  2:1302435/150‑1 (MQ=255)
                                                                                                aaCGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACt                                                      <  2:1360004/150‑1 (MQ=255)
                                                                                                  cGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGc                                                    >  2:536665/1‑150 (MQ=255)
                                                                                                    tCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTc                                                  >  1:326860/1‑150 (MQ=255)
                                                                                                     cAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTcc                                                 >  2:568250/1‑150 (MQ=255)
                                                                                                        aGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACGGTAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTACTCACCCTATTGGAACTGCTCCCTg                                              <  2:823920/150‑1 (MQ=255)
                                                                                                         gCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGc                                             <  1:1029510/150‑1 (MQ=255)
                                                                                                          ccTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCa                                            >  1:1366064/1‑150 (MQ=255)
                                                                                                             gATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCaaaa                                         >  2:1229664/1‑150 (MQ=255)
                                                                                                              aTGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACACTATTGGAACTGCTACCTGCAAAAc                                        >  1:1511945/1‑150 (MQ=255)
                                                                                                                ggAAAAACTGGCGCTGAACAAAGATGCAGTGGTGCTGGTTATCAGCAATGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACgg                                      >  1:970499/1‑150 (MQ=255)
                                                                                                                 gAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACggg                                     <  2:1312228/150‑1 (MQ=255)
                                                                                                                 gAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACggg                                     >  1:1354583/1‑150 (MQ=255)
                                                                                                                 gAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATGGAACTGCTCCCTGCAAAACgggg                                    <  2:699988/150‑1 (MQ=255)
                                                                                                                  aaaaaCTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACgggg                                    <  1:1021902/150‑1 (MQ=255)
                                                                                                                  aaaaaCTGGCGCTGAACAAAGATGCCGTGGTGCGGGTTATCAGTACTGAAGGCGACACCGACGTGAGGCACTACCGCGAAGTTGTCTGGGAGGGCAAACACGCTGTAGCACCTTTATCACCCTATTGGAACTGCTCCCTGCAAAACgggg                                    >  2:1642106/1‑150 (MQ=255)
                                                                                                                    aaaCTGGCGCTTAACAAAGATGCCGTGGGGCTGGTTATCAGCAGTGAAGGCGACACCGACGTGAATCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGTGGAg                                  <  2:283187/150‑1 (MQ=255)
                                                                                                                    aaaCTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAg                                  <  1:1296636/150‑1 (MQ=255)
                                                                                                                    aaaCTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAg                                  <  1:614887/150‑1 (MQ=255)
                                                                                                                          gcgcTGAACAAAGATGCCGTGGTGCGGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTaaaaa                            <  1:1573823/150‑1 (MQ=255)
                                                                                                                            gcTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAAt                          >  2:857181/1‑150 (MQ=255)
                                                                                                                               gAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCCCCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTg                       >  2:1300889/1‑150 (MQ=255)
                                                                                                                                aaCAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTgg                      <  1:812980/150‑1 (MQ=255)
                                                                                                                                  cAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGgag                    <  2:1223685/150‑1 (MQ=255)
                                                                                                                                     aGATGCCGTGGTGCTGGTTATCAGCACTGAAGCCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCAGGAGaaa                 <  1:430051/150‑1 (MQ=255)
                                                                                                                                       aTGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGaaaaa               >  1:1705192/1‑150 (MQ=255)
                                                                                                                                         ccttcagtgCTGATAAACAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCACTCACGCAGTAGGACCTTAATCACCCTTTTGGACCTGCTCCCTGCAAAACCGGGAGTAAAAAAACTGGAGAAAaata             >  1:981403/7‑150 (MQ=255)
                                                                                                                                           ttcagtgCTGATAACCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAAACTGGAGAAAaataat           >  2:708421/5‑150 (MQ=255)
                                                                                                                                              ggtgCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGc        <  1:61059/150‑1 (MQ=255)
                                                                                                                                              ggtgCTGGTTATCACCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGc        <  1:1440231/150‑1 (MQ=255)
                                                                                                                                              agtgCTGATAACCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGc        >  2:1517147/2‑150 (MQ=255)
                                                                                                                                                tgCTGATAACCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTa      <  1:1517147/150‑1 (MQ=255)
                                                                                                                                                 gCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGCAAAATAATGGCTaa     >  1:55589/1‑150 (MQ=255)
                                                                                                                                                  cTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAg    >  2:1012907/1‑150 (MQ=255)
                                                                                                                                                    ggTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAGaa  >  1:1674328/1‑150 (MQ=255)
                                                                                                                                                    ggTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAGaa  >  1:29985/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAGAA  >  NC_000913/3008477‑3008774

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