Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 388,710 T→C intergenic (+64/+43) tauD → / ← hemB taurine dioxygenase, 2‑oxoglutarate‑dependent/5‑aminolevulinate dehydratase (porphobilinogen synthase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913388,7100TC100.0% 449.0 / NA 115intergenic (+64/+43)tauD/hemBtaurine dioxygenase, 2‑oxoglutarate‑dependent/5‑aminolevulinate dehydratase (porphobilinogen synthase)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (59/56);  total (59/56)

Only 100 of 115 total aligned reads displayed.

TGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTT  >  NC_000913/388562‑388853
                                                                                                                                                    |                                                                                                                                               
tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCAc                                                                                                                                               >  4:71197/1‑151 (MQ=255)
 gCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCAc                                                                                                                                               <  3:76187/150‑1 (MQ=255)
 gCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCAc                                                                                                                                               <  1:581289/150‑1 (MQ=255)
      tGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCa                                                                                                                                          >  1:535522/1‑150 (MQ=255)
       gCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAAc                                                                                                                                        >  1:975972/1‑151 (MQ=255)
          gATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAgg                                                                                                                                     >  1:1050802/1‑151 (MQ=255)
          gATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAgg                                                                                                                                     <  2:1279947/151‑1 (MQ=255)
               ccTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAAc                                                                                                                                <  2:72145/151‑1 (MQ=255)
                cTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACt                                                                                                                               >  2:1162813/1‑151 (MQ=255)
                 tGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTc                                                                                                                              >  1:596621/1‑151 (MQ=255)
                     acaGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTccccc                                                                                                                          <  1:722772/151‑1 (MQ=255)
                      caGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCt                                                                                                                         <  2:541689/151‑1 (MQ=255)
                         cGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTc                                                                                                                      >  3:365307/1‑151 (MQ=255)
                           aCGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTg                                                                                                                    <  1:69511/151‑1 (MQ=255)
                           aCGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTg                                                                                                                    <  1:612431/151‑1 (MQ=255)
                            cGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTga                                                                                                                   <  4:299037/151‑1 (MQ=255)
                              gATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTgaga                                                                                                                 <  4:330964/151‑1 (MQ=255)
                               aTAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTgag                                                                                                                  >  1:203828/1‑149 (MQ=255)
                                 aaTGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAgg                                                                                                               <  2:1146201/150‑1 (MQ=255)
                                  aTGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaa                                                                                                             <  2:1249533/151‑1 (MQ=255)
                                    gCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaaa                                                                                                            <  2:923509/150‑1 (MQ=255)
                                    gCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAAc                                                                                                           >  1:341287/1‑151 (MQ=255)
                                    gCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAAc                                                                                                           <  1:699767/151‑1 (MQ=255)
                                        cGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaaa                                                                                                        >  1:378555/1‑150 (MQ=255)
                                         gggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAAt                                                                                                      <  3:259343/151‑1 (MQ=255)
                                           gCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAAtt                                                                                                     >  1:1153184/1‑150 (MQ=255)
                                           gCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTa                                                                                                    <  3:437271/151‑1 (MQ=255)
                                            cgacgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTa                                                                                                    >  3:83620/1‑150 (MQ=255)
                                               cgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                                                                                                <  4:166603/151‑1 (MQ=255)
                                                gaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                                                                                                >  2:196271/1‑150 (MQ=255)
                                                gaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                                                                                                >  2:879362/1‑150 (MQ=255)
                                                  tCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                                                                                                >  4:30381/1‑148 (MQ=255)
                                                    cTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAt                                                                                           >  2:1219789/1‑151 (MQ=255)
                                                      tGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtct                                                                                         >  3:119757/1‑151 (MQ=255)
                                                        gggATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTGACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                                                                                       <  4:131310/151‑1 (MQ=255)
                                                        gggATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                                                                                       <  1:946487/151‑1 (MQ=255)
                                                           aTAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttc                                                                                    >  3:209407/1‑151 (MQ=255)
                                                                 cGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCTCTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGcc                                                                              <  1:932782/151‑1 (MQ=255)
                                                                 cGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGcc                                                                              >  2:43310/1‑151 (MQ=255)
                                                                     ttATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAAt                                                                          >  4:355319/1‑151 (MQ=255)
                                                                      tATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAAt                                                                          <  3:27778/150‑1 (MQ=255)
                                                                       aTCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATcc                                                                        <  2:509588/151‑1 (MQ=255)
                                                                       aTCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATcc                                                                        >  2:86847/1‑151 (MQ=255)
                                                                       aTCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATcc                                                                        >  2:393507/1‑151 (MQ=255)
                                                                         cgggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAg                                                                      <  1:1067901/151‑1 (MQ=255)
                                                                          gggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAgc                                                                     <  1:1069403/151‑1 (MQ=255)
                                                                          gggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAgc                                                                     <  1:1005308/151‑1 (MQ=255)
                                                                           ggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAgc                                                                     >  3:178869/1‑150 (MQ=255)
                                                                            gcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAgcgc                                                                   <  1:1056823/151‑1 (MQ=255)
                                                                             cgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCa                                                                  <  1:297003/151‑1 (MQ=255)
                                                                              ggggTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCaa                                                                 <  2:41460/151‑1 (MQ=255)
                                                                                ggTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCaa                                                                 >  1:224358/1‑149 (MQ=255)
                                                                                   aaTACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTa                                                             <  1:991069/150‑1 (MQ=255)
                                                                                     tACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCt                                                          >  2:898597/1‑151 (MQ=255)
                                                                                      aCGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTg                                                         >  2:341845/1‑151 (MQ=255)
                                                                                       cGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGa                                                        <  2:54273/151‑1 (MQ=255)
                                                                                        gagaGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGaa                                                       <  2:123523/151‑1 (MQ=255)
                                                                                        gagaGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGa                                                        >  1:523307/1‑150 (MQ=255)
                                                                                         agagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGaaa                                                      <  1:913807/151‑1 (MQ=255)
                                                                                          gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGaaaa                                                     <  1:533219/151‑1 (MQ=255)
                                                                                             tGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAAt                                                    >  2:695840/1‑149 (MQ=255)
                                                                                              ggACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCa                                                  <  1:154793/150‑1 (MQ=255)
                                                                                               gACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATc                                                   >  2:408055/1‑148 (MQ=255)
                                                                                                aCGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATc                                                   >  4:364739/1‑147 (MQ=255)
                                                                                                aCGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGAt                                               <  2:564790/151‑1 (MQ=255)
                                                                                                 cGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATc                                                   >  1:146644/1‑146 (MQ=255)
                                                                                                  ggTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTGTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATcc                                             >  2:826494/1‑151 (MQ=255)
                                                                                                  ggTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATcc                                             <  1:629364/151‑1 (MQ=255)
                                                                                                     ccccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgc                                           <  4:353104/150‑1 (MQ=255)
                                                                                                      cccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcac                                         <  1:94481/151‑1 (MQ=255)
                                                                                                          cGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccgc                                     >  2:606178/1‑151 (MQ=255)
                                                                                                           gCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccgc                                     >  2:94530/1‑150 (MQ=255)
                                                                                                           gCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccgc                                     >  4:451525/1‑150 (MQ=255)
                                                                                                           gCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccgc                                     >  2:723188/1‑150 (MQ=255)
                                                                                                            cccccTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccgcac                                   <  4:6287/151‑1 (MQ=255)
                                                                                                            cccccTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccgcac                                   <  1:36945/151‑1 (MQ=255)
                                                                                                               ccTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCt                                >  2:1037070/1‑151 (MQ=255)
                                                                                                                cTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCtt                               >  2:742651/1‑151 (MQ=255)
                                                                                                                  tGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTa                              <  2:690729/150‑1 (MQ=255)
                                                                                                                    gggAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATc                           >  2:107824/1‑151 (MQ=255)
                                                                                                                     ggAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATca                          <  2:851096/151‑2 (MQ=255)
                                                                                                                     ggAGAGGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCg                          >  2:334186/1‑151 (MQ=255)
                                                                                                                       agagGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGa                         >  2:101851/1‑150 (MQ=255)
                                                                                                                       agagGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGa                         >  4:86679/1‑150 (MQ=255)
                                                                                                                        gagGGTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAAc                       >  1:647113/1‑151 (MQ=255)
                                                                                                                          gggTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCt                     >  4:177975/1‑151 (MQ=255)
                                                                                                                           ggTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCt                     >  1:593537/1‑150 (MQ=255)
                                                                                                                           ggTTAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCt                     >  3:202453/1‑150 (MQ=255)
                                                                                                                              tAGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGc                 >  1:1074974/1‑151 (MQ=255)
                                                                                                                               aGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCt                <  2:535424/151‑1 (MQ=255)
                                                                                                                               aGGGTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCt                >  1:1237773/1‑151 (MQ=255)
                                                                                                                                 ggTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCttt              <  1:1213917/151‑1 (MQ=255)
                                                                                                                                  gTGAGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTc             <  3:413382/151‑1 (MQ=255)
                                                                                                                                     aGGGGGCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAg          <  2:28273/151‑1 (MQ=255)
                                                                                                                                      gggggCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGc         >  4:89438/1‑151 (MQ=255)
                                                                                                                                          gCGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACg      >  2:747209/1‑150 (MQ=255)
                                                                                                                                           cGTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGAc    <  4:269814/151‑1 (MQ=255)
                                                                                                                                            gTTCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACt   >  1:641634/1‑151 (MQ=255)
                                                                                                                                             ttCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACtt  >  1:787351/1‑151 (MQ=255)
                                                                                                                                             ttCACCGCACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACtt  >  1:901322/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTT  >  NC_000913/388562‑388853

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