Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 7.3% intergenic (‑82/+36) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,8550CA 7.3% 413.9 / 12.5 142intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (57/73);  new base A (6/5);  total (63/79)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

TATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAC  >  NC_000913/4366708‑4367004
                                                                                                                                                   |                                                                                                                                                      
tATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTg                                                                                                                                                      >  2:595499/1‑150 (MQ=255)
   acccaccACTCTGGCGCCCAAGAACATCGCCAGCAAAACTATGATGAGTTCTACAACTATCATATTCGCCTTCCTTGTTTTTTCACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACCAGAAAAAAGGCACGTCATCTGACg                                                                                                                                                   <  2:1426355/149‑1 (MQ=255)
     cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTTCTGCTCACAAGAAAAAAGGCACGTCTTCTGACGTg                                                                                                                                                 >  2:333659/1‑150 (MQ=255)
     cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                                 <  1:1505344/150‑1 (MQ=255)
     cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                                 >  2:1806583/1‑150 (MQ=255)
     cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                                 >  2:1751077/1‑150 (MQ=255)
     cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                                 <  1:1202889/150‑1 (MQ=255)
       cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGcc                                                                                                                                               <  1:851904/150‑1 (MQ=255)
       cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGcc                                                                                                                                               <  2:696471/150‑1 (MQ=255)
        ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑t                                                                                                                                             >  1:1292857/1‑150 (MQ=255)
         cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCC‑tt                                                                                                                                            <  1:1657180/150‑1 (MQ=255)
          aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCC‑ttt                                                                                                                                           <  2:1304807/150‑1 (MQ=255)
             cTGTCGCCCAAGAAGATCGACAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑tttttt                                                                                                                                        <  1:197607/150‑1 (MQ=255)
               ggcggcAAAGAAGTTCGCCATCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTAt                                                                                                                                      <  2:1358015/150‑1 (MQ=255)
                gcgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTAtt                                                                                                                                     <  2:1321666/150‑1 (MQ=255)
                 cgcCCCAGAAGCACGCCCGCAAAACTATGATGAGGTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGATGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAACGCACGTCAGATGACGTGCC‑TTTTTTCttg                                                                                                                                    <  2:268019/150‑2 (MQ=255)
                 cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTAttt                                                                                                                                    <  1:839380/150‑1 (MQ=255)
                 cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTAttt                                                                                                                                    <  2:615903/150‑1 (MQ=255)
                 cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCC‑TTTTTTCttg                                                                                                                                    <  1:551979/150‑2 (MQ=255)
                   cccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGt                                                                                                                                  >  1:1458853/1‑150 (MQ=255)
                     cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGtac                                                                                                                                >  1:823805/1‑150 (MQ=255)
                     cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGtaa                                                                                                                                >  2:1779010/1‑149 (MQ=255)
                     cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAACATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGtac                                                                                                                                >  2:829362/1‑150 (MQ=255)
                     cAAGAAGATAGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTTCA‑TTTTTTATTTGttc                                                                                                                                >  1:816825/1‑148 (MQ=255)
                          agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTAccc                                                                                                                           <  1:710119/150‑1 (MQ=255)
                           gaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCt                                                                                                                          >  2:1232462/1‑150 (MQ=255)
                           gaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCt                                                                                                                          >  2:1181502/1‑150 (MQ=255)
                            aTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTg                                                                                                                         >  1:450892/1‑150 (MQ=255)
                             tCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGt                                                                                                                        >  1:1475047/1‑150 (MQ=255)
                                ccAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTTTACTACCCTGTACg                                                                                                                     >  1:238550/1‑150 (MQ=255)
                                 cAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGa                                                                                                                    >  2:453529/1‑150 (MQ=255)
                                  aGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGAt                                                                                                                   <  1:1547303/150‑1 (MQ=255)
                                   gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGAtt                                                                                                                  >  2:1663386/1‑150 (MQ=255)
                                     aaaaCTATGATGAGTTCTACAACTAGCATATTATCCTTCCTTGTTTTTTATCAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTAc                                                                                                                <  1:912392/150‑1 (MQ=255)
                                     aaaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGATGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCA‑TTTTTTATTTGTACTACCCTGTACGATTAc                                                                                                                >  2:734381/1‑150 (MQ=255)
                                       aaCTATGATGAGTTCTACAACTTGCATATTAGCCTTCCTTGTTTTTTATCAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTg                                                                                                              <  1:1779678/150‑1 (MQ=255)
                                       aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTg                                                                                                              <  1:1786544/150‑1 (MQ=255)
                                       aaCTATGATGAGTTCTACAACTAGCATAATAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTg                                                                                                              <  2:771970/150‑1 (MQ=255)
                                           atgatgAGTTCTACAACTAGCTTATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCGGACGTGCC‑TTTTTTATTTCTACTACCCTGTACGATTACTGTTCg                                                                                                          >  2:1098875/1‑150 (MQ=255)
                                           atgatgAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCg                                                                                                          >  2:1737323/1‑150 (MQ=255)
                                              atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTATTTAACAACTTGATAGTAGATTGATATTTTTAAGTTACTTCTCACAATAAAAAAGGCACGTCATCTGACTTGCC‑TTTTTTATTTGTACTACCCTTAACGATTACTATTCTCtt                                                                                                       >  1:79869/1‑150 (MQ=255)
                                               tgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                      <  1:458509/150‑1 (MQ=255)
                                                 aGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                    >  1:1482646/1‑150 (MQ=255)
                                                  gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                   <  1:725313/150‑1 (MQ=255)
                                                   ttcttcAAATAGAATATTAGCCTTCCTTTTTTTTTAACTAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                  <  1:339092/150‑1 (MQ=255)
                                                     cTACAACTAGCATATTAGCATTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATATGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTAGCTATCATCAg                                                                                                >  2:1045621/1‑150 (MQ=255)
                                                     atACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                                <  1:441324/149‑1 (MQ=255)
                                                        cAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtat                                                                                             <  1:1124646/150‑1 (MQ=255)
                                                             aGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGt                                                                                        >  1:695639/1‑150 (MQ=255)
                                                                 tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTg                                                                                    >  1:575303/1‑150 (MQ=255)
                                                                    tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                 <  2:257928/150‑1 (MQ=255)
                                                                    tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                 <  2:1190062/150‑1 (MQ=255)
                                                                    tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                 <  2:1573185/150‑1 (MQ=255)
                                                                      gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTg                                                                               <  2:846765/150‑1 (MQ=255)
                                                                      gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTg                                                                               <  2:509265/150‑1 (MQ=255)
                                                                       ccttccttGTTTTTTCACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACACTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGt                                                                              <  1:85633/150‑1 (MQ=255)
                                                                       ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGt                                                                              >  1:1452683/1‑150 (MQ=255)
                                                                        cttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTa                                                                             <  1:1232462/150‑1 (MQ=255)
                                                                        cttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTa                                                                             <  1:1181502/150‑1 (MQ=255)
                                                                           ccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGc                                                                          >  2:812062/1‑150 (MQ=255)
                                                                           acttGTTTTTTAACAATTTGATATTAGTTTGTTATTTTTAAGTTACTGCTCAGAAGAAAAAAGGCACTTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGc                                                                          <  1:59570/149‑1 (MQ=255)
                                                                            cttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGAC‑‑TTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCg                                                                        >  1:1675200/1‑150 (MQ=255)
                                                                              tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                       >  2:1012408/1‑150 (MQ=255)
                                                                                 tttttAACAAGTTGATTTTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTg                                                                    <  2:1541084/150‑1 (MQ=255)
                                                                                 tttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTg                                                                    <  2:695639/150‑1 (MQ=255)
                                                                                    ttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                 >  1:1073070/1‑150 (MQ=255)
                                                                                    ttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                 >  2:1034399/1‑150 (MQ=255)
                                                                                     tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGTAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATc                                                                <  2:1021024/150‑1 (MQ=255)
                                                                                     tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATc                                                                <  1:898206/150‑1 (MQ=255)
                                                                                       aCAAGTAGATATTAGATTGTTATTTTTAAGTTAATGCTCACAAGAAAAAAGGCACGTCATCTCACGTGCC‑TTTTTTATATCTACTACCCTGTACGATTACTGTTCGCTTTCATCAGAATAGAGTTTTTCTTTGAAAGCCGGCTGAATCac                                                              >  2:451001/1‑149 (MQ=255)
                                                                                             tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                        <  2:668487/150‑1 (MQ=255)
                                                                                              gATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                       >  1:411946/1‑150 (MQ=255)
                                                                                               atatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                      >  1:811961/1‑150 (MQ=255)
                                                                                                tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTAc                                                     <  2:1292857/150‑1 (MQ=255)
                                                                                                tatTAGATTGTTATTTTTAACTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGGGCATCAGATTCTGTAc                                                     <  2:727408/150‑1 (MQ=255)
                                                                                                    aGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgag                                                 <  1:450447/150‑1 (MQ=255)
                                                                                                    aGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTATTGGTTTGTAAGACTGGTGCATCAGATTCTGTACGgag                                                 >  2:1110483/1‑150 (MQ=255)
                                                                                                         gTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCATTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                            >  2:1145658/1‑150 (MQ=255)
                                                                                                         gTTATTTTTAAGTTACAGCTCACAAGAAAAAAGGCACGACATCTGACGTGCC‑ATTTTAATTTGTACTACCCAGTACGATTACTGTTCGCTTTCATAAGTATAGCGATATACATTCTAAGCCGGGTGCAACAGATTCTGTACCTAGAAAat                                            >  2:676784/1‑150 (MQ=255)
                                                                                                             tttttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCg                                        >  2:1337529/1‑150 (MQ=255)
                                                                                                              ttttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGt                                       <  1:933414/150‑1 (MQ=255)
                                                                                                                ttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCa                                     <  2:1422841/150‑1 (MQ=255)
                                                                                                                ttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCa                                     <  1:734381/150‑1 (MQ=255)
                                                                                                                          gCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCtt                           >  1:276652/1‑150 (MQ=255)
                                                                                                                           cTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                          >  2:1290006/1‑150 (MQ=255)
                                                                                                                           cTCACAAGAAAAAAGGCACGTCATCGGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                          <  2:490848/150‑1 (MQ=255)
                                                                                                                             cacaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTTTTTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTAATGTTCCGCTTCa                         >  2:1024143/1‑150 (MQ=255)
                                                                                                                             cacaAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                        <  1:1582479/150‑1 (MQ=255)
                                                                                                                             cacaAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                        <  2:1000980/150‑1 (MQ=255)
                                                                                                                             cacaAGAAAAAAGGCACGTCATATGACGTGCC‑TTTTTTATTTGTACTACCATGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTTCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                        >  2:1176955/1‑150 (MQ=255)
                                                                                                                               caAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                      >  2:58915/1‑150 (MQ=255)
                                                                                                                               aaaGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTTCTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                      <  1:697173/149‑1 (MQ=255)
                                                                                                                                       aaGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc              <  1:1047672/150‑1 (MQ=255)
                                                                                                                                       aaGGCACGTCAGATGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc              >  1:652467/1‑150 (MQ=255)
                                                                                                                                              gTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTcc       <  1:812062/150‑1 (MQ=255)
                                                                                                                                                cATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg     <  1:792734/150‑1 (MQ=255)
                                                                                                                                                cATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg     <  2:450892/150‑1 (MQ=255)
                                                                                                                                                 aTCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGa    <  2:419329/150‑1 (MQ=255)
                                                                                                                                                   cTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAc  <  1:1751077/150‑1 (MQ=255)
                                                                                                                                                   cTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAc  <  2:1073070/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAC  >  NC_000913/4366708‑4367004

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