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,1550TG96.6% 371.9 / 1.5 116D124E (GAT→GAGpotFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (53/58);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 116 total aligned reads displayed.

AGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCC  >  NC_000913/894058‑894250
                                                                                                 |                                                                                               
aGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGcc                                                                                               <  1:2811854/100‑1 (MQ=255)
   tGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAATGGAAGAATCTCGACCAGGAAATAATGAAGCTGGTCGCCAAACACGAGccc                                                                                              >  1:203239/1‑98 (MQ=255)
    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                           <  1:750554/100‑1 (MQ=255)
    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                           <  1:1235226/100‑1 (MQ=255)
    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                           <  2:2793598/100‑1 (MQ=255)
    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                           <  1:2549960/100‑1 (MQ=255)
    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                           <  1:2146872/100‑1 (MQ=255)
      cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACaa                                                                                         >  2:1141947/1‑100 (MQ=255)
        gCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACAAGATcggaagagca                                                                                       >  1:1566247/1‑91 (MQ=255)
         cGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATa                                                                                       <  2:1887188/99‑1 (MQ=255)
          gggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATaaa                                                                                     >  2:2787328/1‑100 (MQ=255)
           ggAGTTTTCCAGCCGCTGGACAAAAGAAAATTGCCTGAGTGGAAGAATCTCGACCCGGAACTGCTGAATCTGGTCGCCAAACACGAGCCCGACAATaaag                                                                                    <  2:1243935/100‑2 (MQ=255)
              gTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTg                                                                                 >  1:583186/1‑100 (MQ=255)
               ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                <  1:26024/100‑1 (MQ=255)
                tttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                >  1:542237/1‑99 (MQ=255)
                 ttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTa                                                                              <  2:137368/100‑1 (MQ=255)
                  tCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTAt                                                                             <  1:67714/100‑1 (MQ=255)
                   ccAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATg                                                                            >  2:1790770/1‑100 (MQ=255)
                   ccAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATg                                                                            >  1:2857847/1‑100 (MQ=255)
                    cAGCCGCTGGATAAAAGCAAATTGCCGGAGTGGAAGATTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGc                                                                           <  1:1160221/100‑1 (MQ=255)
                       ccGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGccc                                                                         >  2:539326/1‑99 (MQ=255)
                       ccGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCt                                                                        >  1:1231238/1‑100 (MQ=255)
                        cGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCta                                                                       >  1:2416344/1‑100 (MQ=255)
                        cGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCta                                                                       <  2:3025517/100‑1 (MQ=255)
                           tGGACGAAAGCAAATTGCCGGAGTGTAAGCATCTCGACCCGGAACTGCTGAAGCTTTTCCCCAAAGACGAGCCCGACAATAAATTTGCTATGccc                                                                         <  2:2728873/95‑1 (MQ=255)
                             gACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCATATAtgt                                                                  >  2:2632438/1‑100 (MQ=255)
                              aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                 <  2:2426465/100‑1 (MQ=255)
                               cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                                >  1:2671436/1‑100 (MQ=255)
                                aaaaGCAAATTGTCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                               <  1:1399960/100‑1 (MQ=255)
                                aaaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTTGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                               <  2:1681301/100‑1 (MQ=255)
                                aaaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                               <  1:1941453/100‑1 (MQ=255)
                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                              <  2:427362/100‑1 (MQ=255)
                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGACCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                              >  1:2695468/1‑100 (MQ=255)
                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAACAAGACAATAAATTTGCTAtaaaatata                                                                     >  1:2328821/1‑85 (MQ=255)
                                  aaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcg                                                             <  1:1296626/100‑1 (MQ=255)
                                   aGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcga                                                            >  1:997498/1‑100 (MQ=255)
                                   aGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcga                                                            >  2:1243913/1‑100 (MQ=255)
                                      aaaTTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgacga                                                         <  1:3032478/100‑1 (MQ=255)
                                          tGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCgg                                                     <  1:2267596/100‑1 (MQ=255)
                                            ccGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGa                                                   >  1:167013/1‑100 (MQ=255)
                                              ggAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGAtt                                                 >  2:1533048/1‑100 (MQ=255)
                                               gAGTGGAAGAATCTCGACCTGGAACTGCTGAAGCCGGTTGCCAAACACGTGCCTGACAATAAATTTGCTATGCTCTATATGTGGGCGACGACCggg                                                    <  2:2690371/96‑1 (MQ=255)
                                                 gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAAGACGAGCCCGACAATATATTTGCTATGCCCTATATGTGGGCGACGACCGGGATGGGc                                              <  2:593239/100‑1 (MQ=255)
                                                 gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                              <  2:319039/100‑1 (MQ=255)
                                                  tGGAAGAATCTCGACCGGTAACTGCTGAATATGGTCTACAAACCCGAGCCCGACGATAAATTTGTTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                             <  2:312097/100‑1 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                            >  2:2839567/1‑100 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                            >  1:2073052/1‑100 (MQ=255)
                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCAACGACCGGGATTGGCta                                            >  1:1315945/1‑100 (MQ=255)
                                                    gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                           >  1:41898/1‑100 (MQ=255)
                                                    gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                           <  1:994215/100‑1 (MQ=255)
                                                      agaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATaa                                         >  1:1933843/1‑100 (MQ=255)
                                                         aTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGt                                      >  2:654896/1‑100 (MQ=255)
                                                          tctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGtt                                     <  2:2870082/100‑1 (MQ=255)
                                                          tctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGtt                                     <  2:1519416/100‑1 (MQ=255)
                                                           ctcGACCCGGAACTGCTGATGCTGGTCGCCAAACACGTGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                    <  1:1189057/100‑1 (MQ=255)
                                                           ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                    <  2:2860221/100‑1 (MQ=255)
                                                           ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                    >  2:2665544/1‑100 (MQ=255)
                                                           ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                    >  2:241872/1‑100 (MQ=255)
                                                            tcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGa                                   <  2:2166747/100‑1 (MQ=255)
                                                            tcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGa                                   <  1:1965262/100‑1 (MQ=255)
                                                            tcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGa                                   <  1:1068567/100‑1 (MQ=255)
                                                             cGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGAt                                  >  2:2517297/1‑100 (MQ=255)
                                                             cGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGAt                                  >  2:169399/1‑100 (MQ=255)
                                                                 ccGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAg                              <  2:3001618/100‑1 (MQ=255)
                                                                 ccGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAg                              <  2:2382964/100‑1 (MQ=255)
                                                                 ccGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAg                              <  2:802627/100‑1 (MQ=255)
                                                                  cGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGt                             <  2:41898/100‑1 (MQ=255)
                                                                   ggAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGtt                            <  1:2787360/100‑1 (MQ=255)
                                                                    gAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTa                           >  1:1695627/1‑100 (MQ=255)
                                                                      aCTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTaaa                         >  1:863058/1‑100 (MQ=255)
                                                                       ctgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAg                        <  2:1128977/100‑1 (MQ=255)
                                                                        tgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGc                       <  2:1753311/100‑1 (MQ=255)
                                                                         gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCg                      >  1:1639785/1‑100 (MQ=255)
                                                                         gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGAGCGACGACCGGAATTGACTATAACGTTGATAAAGTTAAAGCg                      >  1:2675661/1‑100 (MQ=255)
                                                                          ctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCgg                     <  2:2117517/100‑1 (MQ=255)
                                                                          ctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCgg                     >  1:273896/1‑100 (MQ=255)
                                                                          ctgAAGCTGGTCGCCAAAAACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCgg                     >  2:1176788/1‑100 (MQ=255)
                                                                            gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                   >  2:1373780/1‑100 (MQ=255)
                                                                            gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                   >  1:1874326/1‑100 (MQ=255)
                                                                            gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                   >  2:1071736/1‑100 (MQ=255)
                                                                            gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                   >  1:1564651/1‑100 (MQ=255)
                                                                            gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGACTATAACGTTGATAAAGTTAAAGCgg                     >  1:1382916/1‑98 (MQ=255)
                                                                             aaGCTGGTCGCCAAACACGAGCCTGATAATATATTTGTTATGCCCTATATGTGGGGGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGc                  <  1:1505146/100‑1 (MQ=255)
                                                                               gCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTg                <  1:304364/100‑1 (MQ=255)
                                                                               gCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTg                <  2:2517901/100‑1 (MQ=255)
                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc             >  2:801066/1‑100 (MQ=255)
                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc             <  1:2306015/100‑1 (MQ=255)
                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACAGGGATTGGTTATAACGTTGATAAAGTTAAAGCGGTGCTggg              >  2:617538/1‑99 (MQ=255)
                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCACTATATGTGGGCGACGACCGGGATTGGCTATAACATTGATAAAGTTAAAGCGGTGCTggg              >  2:1870176/1‑99 (MQ=255)
                                                                                   gTCGCCAAACACGAGCCCGACAATAAATTTGCTATGTCCTATATGTGAGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc             <  1:1051618/99‑1 (MQ=255)
                                                                                    tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa           >  1:2435004/1‑100 (MQ=255)
                                                                                    tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa           >  1:514325/1‑100 (MQ=255)
                                                                                      gCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaa         >  1:2314641/1‑100 (MQ=255)
                                                                                      gCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaa         >  1:1056533/1‑100 (MQ=255)
                                                                                       ccAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaaa        >  2:2877779/1‑100 (MQ=255)
                                                                                         aaaCACGAGCCCGACAATAAATTTGCTATGCTCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAc       <  1:1418954/99‑1 (MQ=255)
                                                                                          aaCACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCgcgttttc     <  1:2699336/100‑8 (MQ=255)
                                                                                            cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc   >  1:1606014/1‑100 (MQ=255)
                                                                                            cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc   <  2:2434975/100‑1 (MQ=255)
                                                                                             acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATGACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc  <  2:2675629/100‑1 (MQ=255)
                                                                                                 |                                                                                               
AGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCC  >  NC_000913/894058‑894250

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