Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 388,676 T→G 5.3% intergenic (+30/+77) 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,6760TG 5.3% 249.2 / 2.3 82intergenic (+30/+77)tauD/hemBtaurine dioxygenase, 2‑oxoglutarate‑dependent/5‑aminolevulinate dehydratase (porphobilinogen synthase)
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base G (3/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGC  >  NC_000913/388530‑388822
                                                                                                                                                  |                                                                                                                                                  
gCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTg                                                                                                                                                   <  8:166872/148‑1 (MQ=255)
   aTTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTggg                                                                                                                                                 >  2:161727/1‑147 (MQ=255)
   aTTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAAACGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTggg                                                                                                                                                 >  2:77866/1‑147 (MQ=255)
    tttGGGATAACCGCGGGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGgag                                                                                                                                              <  4:225048/149‑1 (MQ=255)
       gggATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                                                                                                                                           >  4:8205/1‑149 (MQ=255)
        ggATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCCTCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCTCCCCAAGGGGGCGAggg                                                                                                                                          <  5:10572/149‑1 (MQ=255)
        ggATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAggg                                                                                                                                          >  7:76952/1‑149 (MQ=255)
        ggATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                                                                                                                                           <  7:294735/148‑1 (MQ=255)
           tAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGGGAGGGTTa                                                                                                                                       <  5:127281/149‑1 (MQ=255)
             aCCGCGTGACCCAGCACTATGCCAATGCCGAATACCTGCAACAGCGACTGATAATGCATCGGGAGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGAAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg                                                                                                                                     <  2:72279/149‑1 (MQ=255)
              ccGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg                                                                                                                                    <  3:23715/149‑1 (MQ=255)
              ccGCGTGACCCAGCACTATGCCAATGCCGATTACCTACCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg                                                                                                                                    <  4:22703/149‑1 (MQ=255)
                   tGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg                                                                                                                               <  2:62057/149‑1 (MQ=255)
                       ccAGCACTATGCCAATGCCGATTACCTGCCCCAGCGACGGATAATGCATCGGGCGCCGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTGGGGGAGAGGGTTAGGGTGCGGGGGc                                                                                                                           >  5:88941/1‑149 (MQ=255)
                              tataCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCAcc                                                                                                                    >  3:189750/1‑149 (MQ=255)
                                tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCAcc                                                                                                                    <  6:192582/147‑1 (MQ=255)
                                  ccAATGCCGATTACCTGCCACAGCGACGGATCATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTa                                                                                                                 <  3:222875/148‑1 (MQ=255)
                                  ccAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAggg                                                                                                                                          <  5:247864/123‑1 (MQ=255)
                                  ccAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAggg                                                                                                                                          >  6:247864/1‑123 (MQ=255)
                                   cAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGCCGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACt                                                                                                               >  4:91923/1‑149 (MQ=255)
                                   cAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCTCCCCAAGGGGGCGAggg                                                                                                                                          >  5:147231/1‑122 (MQ=255)
                                   cAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCTCCCCAAGGGGGCGAggg                                                                                                                                          <  6:147231/122‑1 (MQ=255)
                                          gATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCaa                                                                                                          >  4:223434/1‑147 (MQ=255)
                                              aCCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGt                                                                                                     <  8:76952/148‑1 (MQ=255)
                                                cTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCTTTTTATCGGGCGGGGTAATACGTGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCa                                                                                                           <  7:79998/140‑1 (MQ=255)
                                                cTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCa                                                                                                           >  8:79998/1‑140 (MQ=255)
                                                  gCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACt                                                                                                >  7:337404/1‑149 (MQ=255)
                                                  gCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACt                                                                                                >  5:273930/1‑149 (MQ=255)
                                                        gCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGTAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCt                                                                                          >  6:208949/1‑149 (MQ=255)
                                                            cGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCt                                                                                      <  8:337404/149‑1 (MQ=255)
                                                              ggtAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGAGGTGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGATGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTga                                                                                    <  5:322502/147‑1 (MQ=255)
                                                              gATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTga                                                                                    <  4:189750/149‑1 (MQ=255)
                                                                  aTGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGagag                                                                                 <  5:208949/148‑1 (MQ=255)
                                                                   tGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGTTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAgg                                                                                >  2:367953/1‑148 (MQ=255)
                                                                      aTCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaaa                                                                             >  3:217951/1‑148 (MQ=255)
                                                                      aTCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaaa                                                                             >  5:271579/1‑148 (MQ=255)
                                                                         gggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaa                                                                          <  6:196769/148‑1 (MQ=255)
                                                                          ggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaaaa                                                                        >  4:337353/1‑149 (MQ=255)
                                                                          ggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaaaa                                                                        <  1:189877/149‑1 (MQ=255)
                                                                          ggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaaaa                                                                        <  7:237305/149‑1 (MQ=255)
                                                                           gCGACGATCCTTGGGGAGAACACGTTTTCTCGGGCGGGGTAATACGAGCGTGGACGGTCCCCTCGCCCCCTTGTGGAGAGGGTTATGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaaaa                                                                        <  4:89604/148‑1 (MQ=255)
                                                                             gacgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaaa                                                                             >  8:16147/1‑141 (MQ=255)
                                                                             gacgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaaa                                                                             <  7:16147/141‑1 (MQ=255)
                                                                                  tCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACg                                                                  >  7:210560/1‑147 (MQ=255)
                                                                                     ttGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGCGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGaa                                                             >  6:198406/1‑149 (MQ=255)
                                                                                     ttGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGaa                                                             >  5:338875/1‑149 (MQ=255)
                                                                                     ttGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGaa                                                             >  2:386921/1‑149 (MQ=255)
                                                                                       ggggATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGaa                                                             <  1:161727/147‑1 (MQ=255)
                                                                                          gATAAACCGTTTTATCGGGCGTGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                                                        <  5:280898/149‑1 (MQ=255)
                                                                                          gATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                                                        <  6:338875/149‑1 (MQ=255)
                                                                                          gATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                                                        <  3:337353/149‑1 (MQ=255)
                                                                                              aaCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttct                                                    >  7:281204/1‑149 (MQ=255)
                                                                                              aaCCGTTTTATCGGGCAGGGTAATACGAGCGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGTGTTAGGGTGAGGGGGCGTTCCCCGTACATTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttc                                                     >  8:155872/1‑148 (MQ=255)
                                                                                               aCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttct                                                    >  5:286929/1‑148 (MQ=255)
                                                                                               aCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttctt                                                      <  1:385970/146‑1 (MQ=255)
                                                                                               aCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttctt                                                      >  2:385970/1‑146 (MQ=255)
                                                                                                   ttttATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTgaga                                                                                  <  5:359698/114‑1 (MQ=255)
                                                                                                   ttatATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTgaga                                                                                  >  6:359698/4‑114 (MQ=255)
                                                                                                       aTCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGCGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTACCGCAGAATCTTCTTCTCAGCCAAAt                                           >  5:258911/1‑149 (MQ=255)
                                                                                                           ggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttct                                                    >  3:153119/1‑136 (MQ=255)
                                                                                                           ggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttct                                                    <  4:153119/136‑1 (MQ=255)
                                                                                                           ggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAg                                       <  3:8205/149‑1 (MQ=255)
                                                                                                                  tactACGAGAGTGGACGGGCCCCTCGCCCCCTTGCGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                <  1:155837/146‑1 (MQ=255)
                                                                                                                  tAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                >  7:294303/1‑149 (MQ=255)
                                                                                                                    aTACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                >  6:13916/1‑147 (MQ=255)
                                                                                                                        gagaGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCt                           >  2:315628/1‑148 (MQ=255)
                                                                                                                        gagaGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCt                           >  3:320890/1‑148 (MQ=255)
                                                                                                                            gTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGaaaa                      <  3:223434/149‑1 (MQ=255)
                                                                                                                               gACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCTGCGCAAAGTAGCTGAAAATc                    <  2:17248/148‑1 (MQ=255)
                                                                                                                               gACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAAt                     >  5:151760/1‑147 (MQ=255)
                                                                                                                               gACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATc                    <  3:68343/148‑1 (MQ=255)
                                                                                                                                ccctaacccTCTCCCCAAGGGGGCGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAAt                     >  7:107026/7‑146 (MQ=255)
                                                                                                                                  ggTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAAt                     >  8:312674/1‑144 (MQ=255)
                                                                                                                                      cccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCTCAAAGTAGCTGAAAATCAGATccgc            <  6:286929/149‑1 (MQ=255)
                                                                                                                                      cccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAGCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccg             >  1:10513/1‑148 (MQ=255)
                                                                                                                                       ccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgca           <  6:271579/149‑1 (MQ=255)
                                                                                                                                         tCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacc         >  3:345369/1‑149 (MQ=255)
                                                                                                                                         tCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacc         >  1:303869/1‑149 (MQ=255)
                                                                                                                                           gCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcccccg       >  3:37340/1‑149 (MQ=255)
                                                                                                                                           gCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccg       >  2:230621/1‑149 (MQ=255)
                                                                                                                                                cTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGACAAGCAGATCCGCACCCGCACGc  >  4:35335/1‑149 (MQ=255)
                                                                                                                                                cTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGc  <  5:198406/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGC  >  NC_000913/388530‑388822

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: ‑