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.3% 907.9 / ‑4.7 296S324P (TCC→CCC) speAbiosynthetic arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (168/126);  total (169/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 305 total aligned reads displayed.

GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACA  >  NC_000913/3082871‑3083049
                                                                                          |                                                                                        
 cGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTg                                                                                       >  2:1661624/1‑93 (MQ=255)
 cGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTg                                                                                       <  1:5360336/93‑1 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGc                                                                                      >  1:5724507/1‑93 (MQ=255)
   gCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCg                                                                                     >  2:3122150/1‑93 (MQ=255)
    cAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGa                                                                                    >  1:4117796/1‑93 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                   >  1:3076567/1‑93 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                   <  1:6064269/93‑1 (MQ=255)
         cGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAg                                                                               >  2:3758574/1‑93 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                              >  1:289318/1‑93 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGACGCCGTAGTTCACCGAACAGTCGGGATGCGAACGAga                                                                              >  1:2464511/1‑92 (MQ=255)
           ttttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTa                                                                             >  2:2801245/1‑93 (MQ=255)
           gtttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTa                                                                             <  1:2331524/92‑1 (MQ=255)
            tttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTAc                                                                            >  1:493992/1‑93 (MQ=255)
             ttcttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTAcc                                                                           >  1:2617944/1‑93 (MQ=255)
               cttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCtt                                                                         >  1:6267882/1‑93 (MQ=255)
                ttcACACGCATCGCCAAGCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACACTCGGGCTCCGAACGAGTACCTTc                                                                        <  2:3686100/93‑1 (MQ=255)
                   acacGCACCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATa                                                                     <  1:7301734/93‑1 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                    >  2:1978826/1‑93 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                    <  1:1837393/93‑1 (MQ=255)
                     acGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAAt                                                                   >  2:4971171/1‑93 (MQ=255)
                     acGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAAt                                                                   >  1:1859266/1‑93 (MQ=255)
                        cATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCa                                                                <  2:6611419/93‑1 (MQ=255)
                         aTCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCAc                                                               >  1:6853041/1‑93 (MQ=255)
                          tCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACg                                                              >  2:4840097/1‑93 (MQ=255)
                            gCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTTCGAACGAGTACCTTCATAATCCACGcc                                                            >  2:7688262/1‑93 (MQ=255)
                             ccAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGccc                                                           <  2:5802679/93‑1 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                          >  1:7917770/1‑93 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                          >  1:4623914/1‑93 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                          <  2:3051023/93‑1 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                          >  1:5168296/1‑93 (MQ=255)
                               aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAg                                                         >  1:616618/1‑93 (MQ=255)
                               aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAg                                                         >  1:822521/1‑93 (MQ=255)
                                aTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGa                                                        <  1:1557567/93‑1 (MQ=255)
                                aTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGa                                                        <  1:7558123/93‑1 (MQ=255)
                                 tCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAc                                                       >  2:3463154/1‑93 (MQ=255)
                                    cccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgc                                                    >  2:1366178/1‑93 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                   <  2:7090891/93‑1 (MQ=255)
                                      cAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccg                                                  >  2:1161886/1‑93 (MQ=255)
                                      cAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccg                                                  >  2:6875032/1‑93 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                 <  2:6954775/93‑1 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                 <  1:2161150/93‑1 (MQ=255)
                                        gATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                                                >  2:2544501/1‑93 (MQ=255)
                                         aTAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgccg                                               <  2:2381045/93‑1 (MQ=255)
                                         aTAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgccg                                               >  2:5014819/1‑93 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgccgc                                              >  1:3643102/1‑92 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                              <  2:7315634/93‑1 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                              <  1:1217380/93‑1 (MQ=255)
                                            aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACg                                            >  2:2157082/1‑93 (MQ=255)
                                             tgttgtTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGATCGCCGCCGACGt                                           >  2:7912499/1‑93 (MQ=255)
                                             tgttgtTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGt                                           <  2:4769526/93‑1 (MQ=255)
                                              gttgttGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTc                                          >  2:6719263/1‑93 (MQ=255)
                                               ttgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCg                                         <  2:6551913/93‑1 (MQ=255)
                                                  ttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGCCGAAg                                      >  1:295456/1‑93 (MQ=255)
                                                   tgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGc                                     >  2:5908253/1‑93 (MQ=255)
                                                   tgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGc                                     >  1:878683/1‑93 (MQ=255)
                                                    ggCGTATTCATTGAGTCCGGAGGTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCa                                    <  1:4198777/93‑1 (MQ=255)
                                                     gCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCAc                                   <  2:2884024/93‑1 (MQ=255)
                                                      cGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACt                                  <  2:4052154/93‑1 (MQ=255)
                                                      cGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACt                                  <  2:1524451/93‑1 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                >  1:6592614/1‑93 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                >  2:1747805/1‑93 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                >  1:2737636/1‑93 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCAAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                >  1:1407098/1‑93 (MQ=255)
                                                         aTTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa                               >  1:5981262/1‑93 (MQ=255)
                                                         aTTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa                               <  2:4159666/93‑1 (MQ=255)
                                                         aTTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa                               <  1:3316004/93‑1 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGGCGCCGACGTCGAAGCACTGAATAtt                           >  2:5614831/1‑93 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                           >  1:5792698/1‑93 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                           >  2:3330743/1‑93 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                           >  1:1029177/1‑93 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                          >  2:5727736/1‑93 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                          >  1:1384607/1‑93 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                          >  2:617295/1‑93 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                         <  2:745431/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                         <  1:3532619/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                         <  1:284373/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                         <  1:7826638/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                         <  2:5095245/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATACACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                         <  1:5252960/93‑1 (MQ=255)
                                                                   gCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGcc                     >  2:954774/1‑93 (MQ=255)
                                                                   gCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGcc                     <  1:1502858/93‑1 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAg                  >  2:5203796/1‑93 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAg                  <  1:32093/93‑1 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAg                  >  2:7214390/1‑93 (MQ=255)
                                                                       tAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGc                 >  2:4422795/1‑93 (MQ=255)
                                                                          ttCACCGAACAGTCGGGCTGCGACCGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtg              <  2:4298911/93‑1 (MQ=255)
                                                                          ttCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGGATATTGACGCCCAGCttt              >  2:3387154/1‑92 (MQ=255)
                                                                             aCCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGc           >  2:6771922/1‑93 (MQ=255)
                                                                              ccGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCa          <  1:774657/93‑1 (MQ=255)
                                                                              ccGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCa          <  2:1213087/93‑1 (MQ=255)
                                                                              ccGAACAGACGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGCAGCACTGAATATTGACGCCCAGCTTGTGCa          >  1:917140/1‑93 (MQ=255)
                                                                               cGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAg         >  1:6111415/1‑93 (MQ=255)
                                                                               cGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAg         >  1:3883120/1‑93 (MQ=255)
                                                                                  aCAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTc      <  1:7119147/93‑1 (MQ=255)
                                                                                  aCAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTc      <  1:4500250/93‑1 (MQ=255)
                                                                                    aGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCca    <  1:5055491/93‑1 (MQ=255)
                                                                                    aGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCca    <  2:7515516/93‑1 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac   >  2:6313003/1‑93 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac   <  1:5559542/93‑1 (MQ=255)
                                                                                      tCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcaca  <  1:1959350/93‑1 (MQ=255)
                                                                                          |                                                                                        
GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACA  >  NC_000913/3082871‑3083049

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