Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,557,869 C→G 7.6% V271L (GTG→CTG)  eutB ← ethanolamine ammonia‑lyase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,557,8690CG 7.6% 302.6 / 10.7 90V271L (GTG→CTG) eutBethanolamine ammonia‑lyase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (43/40);  new base G (2/5);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ATCGTAATGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCC  >  NC_000913/2557723‑2558016
                                                                                                                                                  |                                                                                                                                                   
aTCGTAATGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTgc                                                                                                                                                  >  2:1158189/1‑150 (MQ=255)
  cGTAATGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTCTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTgcgc                                                                                                                                                <  1:1046677/150‑1 (MQ=255)
  cGTAATGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTCTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTgcgc                                                                                                                                                >  2:1046677/1‑150 (MQ=255)
       tGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCGCAGTgggcgcgga                                                                                                                                           <  2:1252103/150‑3 (MQ=255)
       tGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCGCAGTgggcgcgga                                                                                                                                           <  2:1258043/150‑3 (MQ=255)
       tGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCt                                                                                                                                           <  2:987861/150‑1 (MQ=255)
       tGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCt                                                                                                                                           <  2:1132532/150‑1 (MQ=255)
        gACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCtt                                                                                                                                          <  2:291547/150‑1 (MQ=255)
         ccgcgCCAGCCCGTAGTTACGTTCTTCCTTCGTTACCTGGTCTGCGCCGATGTTAGCGCCAGCGTATAGCGCAGAGCCTTGTCCGTTTTCTAAGTAGAGGCTTTTTTCCCCGGCTATACGGTTCAACTCCGCGCGCTGTgggcgcggagt                                                                                                                                         <  2:1330474/149‑5 (MQ=255)
         aCGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACATGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCGCAGTgggcgcggagt                                                                                                                                         <  1:1326727/150‑5 (MQ=255)
          gatggaaGCACGTAACTACGGGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTtcg                                                                                                                                        >  1:221720/7‑150 (MQ=255)
          cgcgCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTtcg                                                                                                                                        >  2:453487/1‑150 (MQ=255)
           atggaaGCACGTAACTACGGGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCACGCGCTtcgt                                                                                                                                       >  2:1402964/6‑150 (MQ=255)
            tggaaGCACGTAACTACGGGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTtcgtc                                                                                                                                      >  1:825277/5‑150 (MQ=255)
             gcCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTtcgtcg                                                                                                                                     <  2:111803/150‑1 (MQ=255)
             gcCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTtcgtcg                                                                                                                                     <  2:626252/150‑1 (MQ=255)
               cAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAg                                                                                                                                   <  1:274776/150‑1 (MQ=255)
                    cGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGGAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAGTGCGCGCGCTTCGTCGAGCATCg                                                                                                                              <  2:387748/150‑1 (MQ=255)
                    cGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCg                                                                                                                              >  2:186985/1‑150 (MQ=255)
                      tAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGcc                                                                                                                            <  2:703011/150‑1 (MQ=255)
                      tAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGcc                                                                                                                            <  1:520141/150‑1 (MQ=255)
                      tAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGcc                                                                                                                            <  1:514471/150‑1 (MQ=255)
                            cGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTcc                                                                                                                      >  2:610449/1‑150 (MQ=255)
                              tGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCTCTTTTTATTCGTTTGAACTCCGCGAGCACTTCTCGATCTTCGTCGAGAATCGTCATCtactc                                                                                                                    <  2:701672/150‑5 (MQ=255)
                                   ccATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCTAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCaa                                                                                                               >  2:178479/1‑150 (MQ=255)
                                    cATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCaaa                                                                                                              >  1:567703/1‑150 (MQ=255)
                                        gTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAActct                                                                                                          >  2:485299/1‑150 (MQ=255)
                                        gTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAActct                                                                                                          >  1:189388/1‑150 (MQ=255)
                                         ttACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCtt                                                                                                         <  1:1158189/150‑1 (MQ=255)
                                           aCCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTc                                                                                                       <  2:55281/150‑1 (MQ=255)
                                             cTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAg                                                                                                     <  1:1666232/150‑1 (MQ=255)
                                                 tCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCt                                                                                                 <  1:200996/150‑1 (MQ=255)
                                                  cTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCtt                                                                                                >  2:1481839/1‑150 (MQ=255)
                                                   tGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCttt                                                                                               >  1:1505569/1‑150 (MQ=255)
                                                     cgcCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCttttt                                                                                             >  2:1020098/1‑150 (MQ=255)
                                                     cgcCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCCCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCttttt                                                                                             >  2:1036370/1‑150 (MQ=255)
                                                      gcCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTc                                                                                            <  1:133769/150‑1 (MQ=255)
                                                          aaGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTg                                                                                        <  1:1430216/150‑1 (MQ=255)
                                                            gTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGcc                                                                                      <  2:982874/150‑1 (MQ=255)
                                                            gTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGcc                                                                                      >  1:837429/1‑150 (MQ=255)
                                                            gTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGcc                                                                                      <  2:1003615/150‑1 (MQ=255)
                                                              tAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGcc                                                                                      >  2:274172/1‑148 (MQ=255)
                                                              tAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGcc                                                                                      <  1:274172/148‑1 (MQ=255)
                                                              tAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCcac                                                                                    <  1:1514596/150‑1 (MQ=255)
                                                                      gCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATActc                                                                            >  1:581861/1‑150 (MQ=255)
                                                                       gggATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATActct                                                                           <  2:1580465/149‑1 (MQ=255)
                                                                        ggATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTg                                                                          <  2:221720/150‑1 (MQ=255)
                                                                           tAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGaa                                                                       >  2:1252661/1‑150 (MQ=255)
                                                                           tAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGaa                                                                       <  1:190021/150‑1 (MQ=255)
                                                                             gcgcAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGaaaa                                                                     <  2:837429/150‑1 (MQ=255)
                                                                              cgcAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAAt                                                                    >  2:562925/1‑150 (MQ=255)
                                                                                     ccTTGTCCGGATTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAgcccg                                                             >  2:1491220/1‑150 (MQ=255)
                                                                                       ttGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAgcccgcc                                                           <  2:19126/150‑1 (MQ=255)
                                                                                           ccGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGaaaa                                                                     <  2:1293224/136‑1 (MQ=255)
                                                                                           ccGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGaaaa                                                                     >  1:1293224/1‑136 (MQ=255)
                                                                                           ccGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGt                                                       >  2:350047/1‑150 (MQ=255)
                                                                                             ggTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTgc                                                     >  1:1136832/1‑150 (MQ=255)
                                                                                                  cGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCac                                                >  2:615550/1‑150 (MQ=255)
                                                                                                      gTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCacgacg                                            <  2:1505569/150‑1 (MQ=255)
                                                                                                      gTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCATTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCacgacg                                            >  2:381930/1‑150 (MQ=255)
                                                                                                       tAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACgg                                           >  1:804018/1‑150 (MQ=255)
                                                                                                          agGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATc                                        <  1:598716/150‑1 (MQ=255)
                                                                                                           ggCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCg                                       >  1:627295/1‑150 (MQ=255)
                                                                                                           ggCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCg                                       >  1:1180606/1‑150 (MQ=255)
                                                                                                           ggCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCg                                       >  1:577074/1‑150 (MQ=255)
                                                                                                             cAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCt                                     >  1:1407984/1‑150 (MQ=255)
                                                                                                             cAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCt                                     >  2:1326590/1‑150 (MQ=255)
                                                                                                             cAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCt                                     >  2:850623/1‑150 (MQ=255)
                                                                                                               gTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTc                                   >  1:517395/1‑150 (MQ=255)
                                                                                                                     cccGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTg                             >  1:146106/1‑150 (MQ=255)
                                                                                                                         gCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGa                                 >  1:1006641/1‑142 (MQ=255)
                                                                                                                         gCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGa                                 <  2:1006641/142‑1 (MQ=255)
                                                                                                                          cGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGggtgg                        <  2:189388/150‑1 (MQ=255)
                                                                                                                          cGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGggtgg                        <  1:850623/150‑1 (MQ=255)
                                                                                                                           gATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGggtggt                       <  2:567703/150‑1 (MQ=255)
                                                                                                                             tACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGa                     <  2:410889/150‑1 (MQ=255)
                                                                                                                             tACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGa                     <  1:1491220/150‑1 (MQ=255)
                                                                                                                               cGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACg                   >  1:86291/1‑150 (MQ=255)
                                                                                                                                    gAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTgcgc              <  1:1020098/150‑1 (MQ=255)
                                                                                                                                    acgaagcgcgcgCAGTGGGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTgcgc              >  1:616450/7‑150 (MQ=255)
                                                                                                                                      gaagcgcgcgCAGTGGGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACGCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAACTCAGCCCGCCCGGTGCGCCACGCCGGCTCGCTTCGATCGGGGTGGTGACGTGCGCCa            >  2:894087/5‑150 (MQ=255)
                                                                                                                                      aCTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCa            >  2:512300/1‑150 (MQ=255)
                                                                                                                                      aCTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTACCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCa            <  1:1481839/150‑1 (MQ=255)
                                                                                                                                       cTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAg           >  2:636062/1‑150 (MQ=255)
                                                                                                                                          cgcgCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCg                                  <  2:830210/124‑1 (MQ=255)
                                                                                                                                          cgcgCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCg                                  >  1:830210/1‑124 (MQ=255)
                                                                                                                                          cgcgCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTac        >  1:268741/1‑150 (MQ=255)
                                                                                                                                          cgcgCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTac        <  1:1237453/150‑1 (MQ=255)
                                                                                                                                           gcgcCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTaca       <  1:1252661/150‑1 (MQ=255)
                                                                                                                                                cACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAACCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGcc  >  1:202276/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCGTAATGACGCGCCAGCCCGTAGTTACGTGCTTCCATCGTTACCTGGTCTGCGCCGAAGTTAGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCTTCGTCGAGCATCGCCAGCTCCACGCCAAACTCTTTCAGCCCTTTTTCGCTGCCACAGATACTCTGGAAAATCAGCCCGCCCGGTGCGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCC  >  NC_000913/2557723‑2558016

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