Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,155 T→G 100% D124E (GAT→GAG potF → putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,1550TG98.2% 418.8 / ‑2.6 122D124E (GAT→GAGpotFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (64/55);  total (66/56)
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

Only 100 of 122 total aligned reads displayed.

AGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGC  >  NC_000913/894058‑894249
                                                                                                 |                                                                                              
aGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGcc                                                                                              <  1:2092347/100‑1 (MQ=255)
 gTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGccc                                                                                             >  1:4115032/1‑100 (MQ=255)
    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                          <  2:3047718/100‑1 (MQ=255)
      cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACaa                                                                                        >  2:3189712/1‑100 (MQ=255)
      cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACaa                                                                                        >  2:623501/1‑100 (MQ=255)
      cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACaa                                                                                        >  2:3977873/1‑100 (MQ=255)
      cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACa                                                                                         <  2:347986/99‑1 (MQ=255)
       tGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAAt                                                                                       <  1:535931/100‑1 (MQ=255)
        gCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATa                                                                                      >  2:3991357/1‑100 (MQ=255)
          gggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATaaa                                                                                    >  1:467199/1‑100 (MQ=255)
           ggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAAt                                                                                   <  2:1440947/100‑1 (MQ=255)
              gTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTg                                                                                >  1:2323090/1‑100 (MQ=255)
               ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                               <  2:2059982/100‑1 (MQ=255)
               ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                               >  2:2536916/1‑100 (MQ=255)
                tttCCAGCCGCTGTACAAAAGCAAATTGCCTGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCt                                                                              <  2:3833539/100‑1 (MQ=255)
                    cAGCCGCTGGACAAAAGCAAGTTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGc                                                                          >  2:2533538/1‑100 (MQ=255)
                    cAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGc                                                                          <  1:960773/100‑1 (MQ=255)
                    cAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGc                                                                          >  1:1146423/1‑100 (MQ=255)
                     cttccgatctcAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGcc                                                                         <  1:621382/90‑1 (MQ=255)
                     agcagcTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGcc                                                                         <  2:3394578/100‑1 (MQ=255)
                     aGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGcc                                                                         >  1:1364632/1‑100 (MQ=255)
                       ccGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCt                                                                       <  1:3414707/100‑1 (MQ=255)
                        cGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCta                                                                      <  2:467195/100‑1 (MQ=255)
                        cGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCta                                                                      >  2:1379675/1‑100 (MQ=255)
                          cTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGTAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCtata                                                                    <  1:762778/100‑1 (MQ=255)
                          cTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCtata                                                                    <  1:2733024/100‑1 (MQ=255)
                            ggACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtg                                                                  <  1:2699897/100‑1 (MQ=255)
                            ggACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtg                                                                  >  2:195052/1‑100 (MQ=255)
                             gACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgt                                                                 <  1:1716184/100‑1 (MQ=255)
                             gACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgt                                                                 >  1:895580/1‑100 (MQ=255)
                              aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                <  2:1495777/100‑1 (MQ=255)
                              aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                <  2:1255414/100‑1 (MQ=255)
                               cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCACACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                               >  2:1874978/1‑100 (MQ=255)
                               cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtg                                                                  >  2:3615493/1‑97 (MQ=255)
                               cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                               >  1:1663831/1‑100 (MQ=255)
                                aaaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                              <  1:1909409/100‑1 (MQ=255)
                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                             >  2:1438721/1‑100 (MQ=255)
                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                             >  1:2515861/1‑100 (MQ=255)
                                   aGCAAATTGCCGGAGTGGAAGAATCTCGACCAGGAACTGCTGAAGCTGGTCGCCAAACACGAGCACGACAATAGATTTGCTATGCCGTATATGTGGGcg                                                            >  2:4034325/1‑99 (MQ=255)
                                    gCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgac                                                          >  1:3160752/1‑100 (MQ=255)
                                    gCAAATTGCCGGAGTGGAAGAATCTCGACCAGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgac                                                          >  1:3067195/1‑100 (MQ=255)
                                      aaaTTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgacga                                                        >  2:1566401/1‑100 (MQ=255)
                                      aaaTTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgacg                                                         >  1:928773/1‑99 (MQ=255)
                                       aaTTGCCGGATTGGAAGAAGCGGGGCCGGGAACTGCTTAATCTGGTCGTCAAACACGAGCCCGACAATAAATTTGCGATGCCCTATATGTGGGcgac                                                          <  2:3810183/97‑1 (MQ=255)
                                          tGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCgg                                                    <  2:2617073/100‑1 (MQ=255)
                                          tGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCgg                                                    <  2:2621319/100‑1 (MQ=255)
                                           gCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGAGACGACCggg                                                   >  2:2944869/1‑100 (MQ=255)
                                            ccGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGa                                                  >  1:3538960/1‑100 (MQ=255)
                                            ccGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGTCGACGACCGGGAt                                                 >  2:3922791/1‑100 (MQ=255)
                                              ggAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGAtt                                                >  1:3983175/1‑100 (MQ=255)
                                               gAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTg                                               <  1:2536930/100‑1 (MQ=255)
                                               gAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTg                                               <  1:491405/100‑1 (MQ=255)
                                                aGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTgg                                              <  2:804706/100‑1 (MQ=255)
                                                 gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                             <  1:436766/100‑1 (MQ=255)
                                                 gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                             <  1:1747508/100‑1 (MQ=255)
                                                  tGGAAGAATCTCGACCCGGAACTGCTGACGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                             >  2:3770903/1‑99 (MQ=255)
                                                  tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                             <  2:121100/99‑1 (MQ=255)
                                                  tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                            >  1:1182697/1‑100 (MQ=255)
                                                  tGGAAGAATCTCGACCCGGAACTGCGGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                            <  1:504659/100‑1 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGACGCTGGTCGCCAACCACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGCTTg                                               >  2:3881660/1‑96 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATCAATTTGCTATGCCCTATATGTGGGGGACGACCGGGATTgctt                                            >  2:2337371/1‑96 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                           >  2:3480358/1‑100 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                           >  2:251089/1‑100 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                           >  2:549051/1‑100 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                           >  2:1539656/1‑100 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTACATGTGGGCGACGACCGGGATTGGc                                             >  2:3907820/1‑98 (MQ=255)
                                                    gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                          <  2:665343/100‑1 (MQ=255)
                                                       gaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAAc                                       <  2:1983957/100‑1 (MQ=255)
                                                        aaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACg                                      <  1:23051/100‑1 (MQ=255)
                                                          tctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGtt                                    <  2:711548/100‑1 (MQ=255)
                                                          tctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGtt                                    <  2:3529671/100‑1 (MQ=255)
                                                           ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACAGGGATTGGCTATAAAGTTg                                   >  1:342173/1‑100 (MQ=255)
                                                            tcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGa                                  <  1:978578/100‑1 (MQ=255)
                                                             cGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGAt                                 >  1:919150/1‑100 (MQ=255)
                                                             cGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGAt                                 >  1:580129/1‑100 (MQ=255)
                                                              gACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATa                                >  2:2985527/1‑100 (MQ=255)
                                                               aCCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGTCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATa                                <  2:3621555/99‑1 (MQ=255)
                                                                 ccGGAACTGCTTAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAg                             <  2:3688783/100‑1 (MQ=255)
                                                                  cGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAAGAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGt                            <  2:2821410/100‑1 (MQ=255)
                                                                      aCTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTaaa                        >  2:2068032/1‑100 (MQ=255)
                                                                        tgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGc                      <  1:427670/100‑1 (MQ=255)
                                                                        tgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGc                      >  1:3460888/1‑100 (MQ=255)
                                                                         gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCg                     >  2:536285/1‑100 (MQ=255)
                                                                             aaGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGc                 <  2:2258056/100‑1 (MQ=255)
                                                                              aGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCt                >  1:3266682/1‑100 (MQ=255)
                                                                              aGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCt                >  1:83781/1‑100 (MQ=255)
                                                                                 tGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTggg             >  1:365736/1‑100 (MQ=255)
                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAGAGCGGTGCTGGGc            <  1:2775886/100‑1 (MQ=255)
                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc            <  1:1903995/100‑1 (MQ=255)
                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc            <  1:847965/100‑1 (MQ=255)
                                                                                   gTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCg           <  2:3990586/100‑1 (MQ=255)
                                                                                   gTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCg           >  2:4055031/1‑100 (MQ=255)
                                                                                    tCGCCCAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGCCGACCGGGATTGGCCAAAACGTCGATAAAGGGAAAGCGGTGcagag             >  2:490646/1‑93 (MQ=255)
                                                                                    tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc            >  2:1404299/1‑98 (MQ=255)
                                                                                     cGCCAAACACGAGCCCGAAAATAAATTTGCTATGCACTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa          >  2:3135387/1‑99 (MQ=255)
                                                                                      gCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaa        <  1:1271409/100‑1 (MQ=255)
                                                                                            cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc  <  1:3635104/100‑1 (MQ=255)
                                                                                            cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc  >  1:2274493/1‑100 (MQ=255)
                                                                                            cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc  <  2:2287091/100‑1 (MQ=255)
                                                                                            cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc  <  1:1121188/100‑1 (MQ=255)
                                                                                                 |                                                                                              
AGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGC  >  NC_000913/894058‑894249

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