Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 387,924 (G)5→4 intergenic (+54/+53) tauD → / ← hemB taurine dioxygenase/5‑aminolevulinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913387,9200G.99.1% 467.3 / ‑3.6 106intergenic (+50/+57)tauD/hemBtaurine dioxygenase/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (94/11);  total (94/12)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGC  >  NC_000913/387828‑388010
                                                                                             |                                                                                           
ggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAg                                                                                             >  2:6245038/1‑93 (MQ=255)
ggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAg                                                                                             >  2:1919951/1‑93 (MQ=255)
 gCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg                                                                                            <  1:4777887/93‑1 (MQ=255)
 gCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg                                                                                            <  2:707244/93‑1 (MQ=255)
 gCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg                                                                                            >  2:4391219/1‑93 (MQ=255)
 gCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCCTGGGGAGAGGGTTAGGGTGAgg                                                                                            >  2:1290588/1‑93 (MQ=255)
 gCGACGATCCGTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg                                                                                            <  1:4365535/93‑1 (MQ=255)
   gacgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGGTAGGGTGAgggg                                                                                          >  2:5133448/1‑93 (MQ=255)
    acgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggg                                                                                         >  2:5468727/1‑93 (MQ=255)
    acgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGCGAggggc                                                                                         >  1:1692466/1‑92 (MQ=255)
    acgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCCTGGGGAGGGGGGTAGGGTGAggggc                                                                                         >  2:3472948/1‑92 (MQ=38)
     cgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggcg                                                                                        <  1:5123372/93‑3 (MQ=255)
      gaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGGTAGGGTGAggggcgt                                                                                       >  2:3747398/1‑90 (MQ=255)
       aTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggcgtt                                                                                      >  1:4927224/1‑89 (MQ=255)
        tCCTGGGGTATAAACCGTTTTATGGGGCGGGGTAATAGGAGAGTGGACGGCCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTc                                                                                   <  2:3179973/93‑1 (MQ=25)
         ccTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCa                                                                                  >  1:3367577/1‑93 (MQ=255)
          cTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGCGAGGGGCGTTCAc                                                                                 >  1:5154668/1‑93 (MQ=255)
           ttGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCGCTCGCCCCCCTGGGGAGAGGGTTAGGGTGAGGGGCGGTCAcc                                                                                >  1:1825369/1‑93 (MQ=37)
           ttGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCAcc                                                                                >  1:5037561/1‑93 (MQ=255)
           ttGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGGTAGGGTGAGGGGCGTTCAcc                                                                                >  2:286789/1‑93 (MQ=255)
             ggggATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGt                                                                              >  1:3422318/1‑93 (MQ=255)
                gATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACt                                                                           >  2:5063322/1‑93 (MQ=255)
                 aTAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACtt                                                                          >  2:2770506/1‑93 (MQ=255)
                  tAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACttt                                                                         >  1:4700038/1‑93 (MQ=255)
                  aaaaCAGTTTTATCGGGAGGGGTAATACGAGAGTGGACGGCCCCCTCGCCCCTTGGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACttt                                                                         <  2:6309533/92‑1 (MQ=25)
                   aaaCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTc                                                                        >  2:3059067/1‑93 (MQ=255)
                    aaCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGGGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCa                                                                       >  1:5267653/1‑93 (MQ=255)
                    aaCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCa                                                                       >  2:366410/1‑93 (MQ=255)
                     aCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCaa                                                                      >  2:4376105/1‑93 (MQ=255)
                     aCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCaa                                                                      >  1:821818/1‑93 (MQ=255)
                       cGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACa                                                                    >  2:3219335/1‑93 (MQ=255)
                       cGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACa                                                                    >  1:4925564/1‑93 (MQ=255)
                       cGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACa                                                                    >  1:6226085/1‑93 (MQ=255)
                        gTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAg                                                                   >  1:575164/1‑93 (MQ=255)
                         ttttATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAgg                                                                  >  1:1408462/1‑93 (MQ=255)
                         ttttATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAgg                                                                  >  2:239580/1‑93 (MQ=255)
                          tttATGGGGCGGGTTAATACGAGAGTGGACGGCCCCCTCGCCCCTTGGGGGAGAGGGTAGGGGTGGGGGGCGTTCACCGTACTTTCAACAGGt                                                                 <  2:1032719/93‑1 (MQ=255)
                                gggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTc                                                           >  1:1944886/1‑93 (MQ=255)
                                 ggcgggGTAATACGAGAGTGGACGTTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTcc                                                          <  2:1059945/93‑1 (MQ=255)
                                     gggTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGGTCACCGTACTTTCAACAGGTTAAATCCCCCt                                                      >  2:2632284/1‑93 (MQ=37)
                                      ggTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCtt                                                     >  1:5208854/1‑93 (MQ=255)
                                       gaaaTACGAGATGGGACGGCCCCCTCGCCCCCTGGGGGAGAGGGTTGGGGTGGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCttt                                                    <  2:3050022/91‑1 (MQ=255)
                                       gTAATACGAGAGGGGACGGTCCCCTCGCCCCCTTGGGGGGAGGGTTAGGGTGATGGGCGTTCACCGTACTTTGTATAGTCTAAATCCCCCttt                                                    >  2:1230171/1‑93 (MQ=21)
                                        tAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTc                                                   >  1:1120303/1‑93 (MQ=255)
                                        tAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTc                                                   >  1:3759951/1‑93 (MQ=255)
                                           tACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTga                                                >  2:5854578/1‑93 (MQ=255)
                                           tACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTga                                                >  1:2897278/1‑93 (MQ=255)
                                             cGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTgaga                                              >  1:2116590/1‑93 (MQ=255)
                                              gagaGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGagag                                             >  1:2542157/1‑93 (MQ=255)
                                                  gTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaaa                                         >  2:605964/1‑93 (MQ=255)
                                                  gTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGaaa                                         >  2:2040394/1‑93 (MQ=255)
                                                   tGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAAc                                        >  2:5917864/1‑93 (MQ=255)
                                                    ggACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACa                                       >  2:4699550/1‑93 (MQ=255)
                                                     gACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaa                                      >  2:3011577/1‑93 (MQ=255)
                                                      gcGGCCCCCCCCCCCCGGTGGGGAGGGGGGGGGGGGGGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaaa                                     <  1:4803379/92‑1 (MQ=255)
                                                       cGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACaaaa                                    >  2:3818390/1‑93 (MQ=255)
                                                         gTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAAtt                                  >  2:5184027/1‑93 (MQ=255)
                                                         gTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAAtt                                  >  2:5428208/1‑93 (MQ=255)
                                                         gTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAAtt                                  >  1:2824472/1‑93 (MQ=255)
                                                          tCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTa                                 >  2:554590/1‑93 (MQ=255)
                                                          tCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTa                                 >  1:3722488/1‑93 (MQ=255)
                                                          cccccTCGCCCCCTGGGGGAGAGGGTTAGGGTGAGGGCCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTa                                 <  2:2160282/92‑1 (MQ=37)
                                                           ccccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTaa                                >  2:3496458/1‑93 (MQ=255)
                                                             ccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACg                              >  2:4780440/1‑93 (MQ=255)
                                                             ccTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACg                              >  2:6040561/1‑93 (MQ=255)
                                                              cTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                             >  1:2906957/1‑93 (MQ=255)
                                                              cTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                             <  1:4941769/93‑1 (MQ=255)
                                                              cTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                             >  1:935423/1‑93 (MQ=255)
                                                              cTCGCCCCCTGGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                             <  1:6276987/93‑1 (MQ=255)
                                                               tCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCa                            >  2:1257673/1‑93 (MQ=255)
                                                               cggCCCCTGTGGGGAGAGGGTTGGGGGGGGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGc                             <  2:1863462/91‑1 (MQ=255)
                                                                cccccccTTGGGGGAGAGGTTTGGGGGGGGGGCCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCa                            <  2:4081540/91‑1 (MQ=255)
                                                                cGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAg                           >  1:284333/1‑93 (MQ=255)
                                                                  cccccTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGaa                         <  1:2067257/93‑1 (MQ=255)
                                                                  cccccTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGaa                         >  2:5789574/1‑93 (MQ=255)
                                                                   ccccTTGGGGAGAGGGTTGGGGTGGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAt                        <  2:5131346/93‑1 (MQ=255)
                                                                   ccccTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAt                        >  1:6128098/1‑93 (MQ=255)
                                                                    cccTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtc                       <  2:4213992/93‑1 (MQ=255)
                                                                     ccTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtct                      >  1:4289547/1‑93 (MQ=255)
                                                                     ccTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtct                      >  2:5127623/1‑93 (MQ=255)
                                                                      cTTGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtctt                     >  1:1913377/1‑93 (MQ=255)
                                                                       ttGGGGAGAGGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTTCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                    >  1:1251169/1‑93 (MQ=255)
                                                                       ttGGGGAGAGGGTTAGGGTGGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                    <  1:1062507/93‑1 (MQ=255)
                                                                       tgggggAGAGGGTGGGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttc                    <  2:1758784/91‑1 (MQ=255)
                                                                         ggggggggggTTGGGGGGGGGGGCCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttct                   <  2:3855230/93‑1 (MQ=255)
                                                                            gggggggTTGGGGGGGGGGGGCCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttc                 <  1:3627888/89‑1 (MQ=255)
                                                                            gagaGGGTTGGGGGGAGGGCCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTtctc               <  2:808100/93‑1 (MQ=255)
                                                                            gagaGGGTTAGGGTGAGGGGAGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttcta               >  1:4093407/1‑92 (MQ=255)
                                                                              gagGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAg             >  2:5445881/1‑93 (MQ=255)
                                                                               agGGTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGc            >  2:2451382/1‑93 (MQ=255)
                                                                                gggTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGcc           >  1:2649743/1‑93 (MQ=255)
                                                                                gggTTAGGGTGAGGGGCGTTCACCGCACTCTCAACAAGTTAACTCCCCCTTTCTGAGAGGAAACAGAATTAAGGCAGAATCTTCTTCTCAGcc           >  2:2435005/1‑93 (MQ=255)
                                                                                  ggTTAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCa          >  2:3719114/1‑93 (MQ=255)
                                                                                    ttAGGGTGAGGGCCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCaaa        <  1:4247105/93‑1 (MQ=255)
                                                                                     tAGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAAt       >  2:447135/1‑93 (MQ=255)
                                                                                      aGGGTGAGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATc      >  1:1516587/1‑93 (MQ=255)
                                                                                       gggTGAGGGCCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATcc     <  1:4942866/93‑1 (MQ=255)
                                                                                         ggtggggggCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCa    <  1:2093463/89‑1 (MQ=255)
                                                                                          gtggggggCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAg   <  1:6024284/90‑1 (MQ=255)
                                                                                           gggggggCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAgc  <  1:3006622/91‑1 (MQ=255)
                                                                                             |                                                                                           
GGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGC  >  NC_000913/387828‑388010

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