Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,082,961 A→G S324P (TCC→CCC)  speA ← biosynthetic arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,9610AG99.4% 1047.2 / ‑3.7 333S324P (TCC→CCC) speAbiosynthetic arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (187/144);  total (187/146)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 351 total aligned reads displayed.

GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAT  >  NC_000913/3082871‑3083050
                                                                                          |                                                                                         
gCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCt                                                                                         <  2:1702892/93‑1 (MQ=255)
 cGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTg                                                                                        >  1:5217859/1‑93 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTgg                                                                                       >  1:8196382/1‑92 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGc                                                                                       <  1:119545/93‑1 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGc                                                                                       <  1:9240934/93‑1 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    <  2:2446800/93‑1 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    <  2:6740544/93‑1 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAAc                                                                                   >  2:1560858/1‑93 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAAc                                                                                   >  1:1761753/1‑93 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAAc                                                                                   >  2:4805726/1‑93 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGACGTAGTTCACCGAACAGTCAGGCTGCGAAc                                                                                   >  2:4150756/1‑93 (MQ=255)
       aCCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACg                                                                                  >  2:1309583/1‑93 (MQ=255)
       aCCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACg                                                                                  <  1:1444024/93‑1 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAga                                                                               >  2:7999669/1‑92 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                               >  1:2052731/1‑93 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                               <  2:2495786/93‑1 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                               <  1:4862842/93‑1 (MQ=255)
            tttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTAc                                                                             <  1:102863/93‑1 (MQ=255)
            tttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTAc                                                                             <  2:4541167/93‑1 (MQ=255)
                 tcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCa                                                                        >  2:3697824/1‑93 (MQ=255)
                 tcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCa                                                                        <  2:9233538/93‑1 (MQ=255)
                   acacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATa                                                                      <  2:6382058/93‑1 (MQ=255)
                   acacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATa                                                                      >  2:7131597/1‑93 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                     >  2:3111078/1‑93 (MQ=255)
                     acGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAAt                                                                    <  2:6278581/93‑1 (MQ=255)
                     acGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAAt                                                                    <  1:5885316/93‑1 (MQ=255)
                     acGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAAt                                                                    >  2:1237541/1‑93 (MQ=255)
                       gCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGTACGAGGACCTTCATAATcc                                                                  >  2:7847143/1‑93 (MQ=255)
                       gCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATcc                                                                  <  1:5125530/93‑1 (MQ=255)
                       gCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATcc                                                                  >  1:2437392/1‑93 (MQ=255)
                          tCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACg                                                               >  2:7495040/1‑93 (MQ=255)
                           cGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGc                                                              <  2:5119503/93‑1 (MQ=255)
                            gCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGcc                                                             >  1:8914089/1‑93 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                           >  2:8136901/1‑93 (MQ=255)
                               aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAg                                                          >  2:3181904/1‑93 (MQ=255)
                                 tCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAc                                                        >  2:1400470/1‑93 (MQ=255)
                                   gCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccg                                                      >  2:6603455/1‑93 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                    <  2:9716857/93‑1 (MQ=255)
                                      cAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccg                                                   <  2:5010963/93‑1 (MQ=255)
                                      cAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccg                                                   >  2:769461/1‑93 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                  <  2:3493549/93‑1 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                  <  1:167684/93‑1 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                  <  1:3806294/93‑1 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAGCAGTCGGGCTGCGAACGAGTACCTTCATAATACACGCCCAGACCGCCGCCGa                                               >  2:2617775/1‑93 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGCCCACCGCCGa                                               <  2:6106189/93‑1 (MQ=255)
                                           aaTGTTGTGGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGAc                                              <  1:4459432/93‑1 (MQ=255)
                                            aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACg                                             >  1:8459758/1‑93 (MQ=255)
                                             tgttgtTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGatgc                                            >  1:541333/1‑90 (MQ=255)
                                               ttgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCg                                          <  1:2615022/93‑1 (MQ=255)
                                               ttgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCg                                          <  2:3081369/93‑1 (MQ=255)
                                                 gttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGaa                                        <  2:8429292/93‑1 (MQ=255)
                                                  ttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAg                                       >  1:5816386/1‑93 (MQ=255)
                                                   tgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGc                                      >  1:135768/1‑93 (MQ=255)
                                                   tgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGc                                      <  2:2074087/93‑1 (MQ=255)
                                                    ggCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCa                                     >  1:7511206/1‑93 (MQ=255)
                                                     gCGTATTCATTGAGGCTGTAGTTCACCGAACAGTCGGGCGGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCTCCGACGTCGAAGCAc                                    <  2:545621/93‑1 (MQ=255)
                                                      cGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACt                                   <  2:5303670/93‑1 (MQ=255)
                                                      cGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACt                                   <  2:8540636/93‑1 (MQ=255)
                                                      cGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACt                                   >  2:3910142/1‑93 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                 <  2:7580157/93‑1 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                 >  1:9506393/1‑93 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                 >  1:4766353/1‑93 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCCACGTCGAAGCACTGa                                 >  1:2311452/1‑93 (MQ=255)
                                                          ttCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAat                               >  1:9761769/1‑93 (MQ=255)
                                                           tccttGTGGCCGTAATTCACCCGACAGTCGGGCTGCGAACGAAGACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAata                              <  1:712546/90‑1 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                            >  1:9052563/1‑93 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                            >  1:6338015/1‑93 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                           >  2:8877886/1‑93 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                          <  2:9468322/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                          <  1:7317185/93‑1 (MQ=255)
                                                                gAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGAc                         >  2:5882400/1‑93 (MQ=255)
                                                                gAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATGGAc                         >  2:7333971/1‑93 (MQ=255)
                                                                  ggCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGc                       >  1:4236112/1‑93 (MQ=255)
                                                                  ggCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGc                       >  1:7768647/1‑93 (MQ=255)
                                                                  ggCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGc                       <  1:9397387/93‑1 (MQ=255)
                                                                   gCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCTCAGACCGCCGTCGACGTCGAAGCACTGAATATTGACGcc                      >  1:3133160/1‑93 (MQ=255)
                                                                   gCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGcc                      <  2:6484103/93‑1 (MQ=255)
                                                                    ccGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGccc                     <  1:5061872/93‑1 (MQ=255)
                                                                     cGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCa                    >  1:4332181/1‑93 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCTGCCGACGTCGAAGCACTGAATATTGACGCCCAg                   >  1:9670718/1‑93 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAg                   >  1:4607125/1‑93 (MQ=255)
                                                                         gTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCtt                <  1:6276681/93‑1 (MQ=255)
                                                                         gTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCtt                >  1:326498/1‑93 (MQ=255)
                                                                         gTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGAtt                >  1:3522777/1‑93 (MQ=255)
                                                                           tCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtgt              <  2:1933175/93‑1 (MQ=255)
                                                                            cACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtgtg             >  1:2165379/1‑93 (MQ=255)
                                                                            cACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtgtg             <  1:815809/93‑1 (MQ=255)
                                                                            cACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtgtg             >  1:289566/1‑93 (MQ=255)
                                                                             aCCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGc            <  1:6210725/93‑1 (MQ=255)
                                                                              ccGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCa           <  2:8346663/93‑1 (MQ=255)
                                                                              ccGAACAGTAGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGACCAGCTCGTGCa           >  1:9893555/1‑93 (MQ=255)
                                                                               cGAACGGTCAGACTGCGAACGAGTACCTACATAATCCACGCCCAGACCTCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAg          <  2:356338/93‑1 (MQ=255)
                                                                               cGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCa                                                                 <  1:5656268/38‑1 (MQ=255)
                                                                                gAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGt         >  1:3911488/1‑93 (MQ=255)
                                                                                 aaCAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGtt        >  1:8745113/1‑93 (MQ=255)
                                                                                  aCAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTc       >  2:224450/1‑93 (MQ=255)
                                                                                    aGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCca     >  2:5765118/1‑93 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac    <  1:3571155/93‑1 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac    >  1:6788450/1‑93 (MQ=255)
                                                                                      tCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcaca   >  2:2782108/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAT  >  NC_000913/3082871‑3083050

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