Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,187 A→T 68.5% D309V (GAC→GTC)  ygeW → putative carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,1870AT68.5% 104.7 / 77.2 89D309V (GAC→GTC) ygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (11/17);  new base T (28/33);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

AGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTG  >  NC_000913/3007089‑3007280
                                                                                                  |                                                                                             
aGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGAc                                                                                              >  2:1855393/1‑100 (MQ=255)
  tGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTg                                                                                            >  1:1467532/1‑100 (MQ=255)
   gATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTgg                                                                                           >  2:1425228/1‑100 (MQ=255)
   gATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTgg                                                                                           <  2:219673/100‑1 (MQ=255)
    aTGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGc                                                                                          <  2:929377/100‑1 (MQ=255)
        aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTg                                                                                      <  1:336924/100‑1 (MQ=255)
         agagCGTACTGAATTGCTGCGTGCGAACGATCACGAACGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGt                                                                                     <  2:1507226/100‑1 (MQ=255)
           agCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTAc                                                                                   <  1:293148/100‑1 (MQ=255)
            gCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACt                                                                                  <  1:1620766/100‑1 (MQ=255)
             cGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTg                                                                                 <  2:812642/100‑1 (MQ=255)
              gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTAc                                                                                   <  1:712493/97‑1 (MQ=255)
              gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                <  2:1502407/100‑1 (MQ=255)
              gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                >  2:1037337/1‑100 (MQ=255)
              gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGa                                                                                <  1:474281/100‑1 (MQ=255)
              gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTAGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                >  2:1243911/1‑100 (MQ=255)
              cttccgatctgctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                <  1:51973/91‑1 (MQ=255)
                aCTGAATTGCTGCGTGCGAACGATCAAGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAag                                                                              >  2:912419/1‑100 (MQ=255)
                  tGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAagag                                                                            >  1:1118089/1‑100 (MQ=255)
                   gAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAgaga                                                                           >  2:1279400/1‑100 (MQ=255)
                     aTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAgatg                                                                         >  1:1913979/1‑100 (MQ=255)
                      ttGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAgatga                                                                        <  2:887853/100‑1 (MQ=255)
                      ttGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAgatga                                                                        <  2:585629/100‑1 (MQ=255)
                        gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAagatcggaa                                                                       >  2:51973/1‑93 (MQ=255)
                        gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGatgatg                                                                      <  1:1279406/100‑1 (MQ=255)
                         ctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATgg                                                                     <  1:942600/100‑1 (MQ=255)
                         ctgcGTAAGAACGATCACGAAAGCCAAAACGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGTATTGTACTGAAGAGATGATgg                                                                     <  2:2240854/100‑1 (MQ=255)
                          tgcgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGa                                                                    >  1:272092/1‑100 (MQ=255)
                             gtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACt                                                                 <  2:1731550/100‑1 (MQ=255)
                             gtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACt                                                                 <  1:963809/100‑1 (MQ=255)
                               gcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGa                                                               >  1:1145151/1‑100 (MQ=255)
                               gcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGa                                                               >  1:106635/1‑100 (MQ=255)
                                cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                              <  2:831962/100‑1 (MQ=255)
                                cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                              >  2:1589081/1‑100 (MQ=255)
                                cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                              <  1:298690/100‑1 (MQ=255)
                                cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                              >  2:966287/1‑100 (MQ=255)
                                 gAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAcc                                                             <  1:65653/100‑1 (MQ=255)
                                  aaCGATTACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                              <  2:849654/98‑1 (MQ=255)
                                   aCGATCACGAAGTCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                           <  1:33016/100‑1 (MQ=255)
                                   aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                           <  2:272093/100‑1 (MQ=255)
                                    cGATCACGAAGGCTTAAAAGCACTGGAAAATCAATGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGGACTGACCCGt                                                          >  2:1470636/1‑100 (MQ=255)
                                    cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                          >  1:17700/1‑100 (MQ=255)
                                    cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                          <  1:2155012/100‑1 (MQ=255)
                                    cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAAATGACCAGt                                                          >  1:1212258/1‑100 (MQ=255)
                                     gATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTg                                                         <  1:1594577/100‑1 (MQ=255)
                                      aTCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCAAAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTg                                                         <  1:1904548/99‑1 (MQ=255)
                                       tCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGAt                                                       <  2:503693/100‑1 (MQ=255)
                                       tCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGAt                                                       <  2:769264/100‑1 (MQ=255)
                                         aCGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGAGAt                                                       >  1:780063/1‑98 (MQ=255)
                                         aCGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATg                                                      >  2:1987027/1‑99 (MQ=255)
                                           gAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCg                                                   >  1:2206008/1‑100 (MQ=255)
                                                 ttAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGccc                                             <  1:996765/100‑1 (MQ=255)
                                                    aaaGCACTGGAAAAACAGTGTATGGCACAGAACGCGCAACACAAAGTCTGGTATTGTACTGAAGAGATGATGGAACTGACTCGTGATGGCGAAGCCCTGt                                          <  2:132555/100‑1 (MQ=255)
                                                       gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACa                                       <  2:813577/100‑1 (MQ=255)
                                                       gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGCCCCGTGATGGCGAAGCCCTGTACa                                       >  2:1092944/1‑100 (MQ=255)
                                                       gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACa                                       <  2:801221/100‑1 (MQ=255)
                                                        cACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACAt                                      >  1:1359292/1‑100 (MQ=255)
                                                        cACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACAt                                      >  1:1154329/1‑100 (MQ=255)
                                                        cACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACAt                                      >  1:5401/1‑100 (MQ=255)
                                                           tGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGc                                    <  2:2127332/99‑1 (MQ=255)
                                                            ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                  >  2:1647968/1‑100 (MQ=255)
                                                            ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                  >  2:2109656/1‑100 (MQ=255)
                                                            ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                  >  1:505142/1‑100 (MQ=255)
                                                               aaaaCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgc                               >  2:2216709/1‑100 (MQ=255)
                                                               aaaaCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgc                               >  2:2216752/1‑100 (MQ=255)
                                                                 aaCAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcct                             <  1:530265/100‑1 (MQ=255)
                                                                  aCAGTGTCTGGCACAGAACGCGCAACACGAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGACGCCCTGTACATGCActgcctg                            >  2:844803/1‑100 (MQ=255)
                                                                  aCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctg                            <  1:113228/100‑1 (MQ=255)
                                                                      tgtCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAg                        <  1:561731/100‑1 (MQ=255)
                                                                        tCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGcn                      <  2:1686449/100‑2 (MQ=255)
                                                                        tCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCt                      <  2:1555955/100‑1 (MQ=255)
                                                                         cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                     <  1:1345448/100‑1 (MQ=255)
                                                                         cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCATGCCAGCTg                     <  2:1141238/100‑1 (MQ=255)
                                                                          tGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGa                    <  2:2130546/100‑1 (MQ=255)
                                                                           ggCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGat                   >  1:1844914/1‑100 (MQ=255)
                                                                                aGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAg              >  1:754508/1‑100 (MQ=255)
                                                                                 aaaCGCGCAACACAAAGTCTGGCAAAGTACAGAAGAGATGATGGAACTGACCCGAGATGACGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCaac             >  2:172759/2‑98 (MQ=255)
                                                                                   aCGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcgg           >  2:1280634/1‑100 (MQ=255)
                                                                                     gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg         <  2:983154/100‑1 (MQ=255)
                                                                                     gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg         <  1:393927/100‑1 (MQ=255)
                                                                                     gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg         <  2:741490/100‑1 (MQ=255)
                                                                                     gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcgg           <  1:969648/98‑1 (MQ=255)
                                                                                     gcgcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg         <  1:874232/100‑1 (MQ=255)
                                                                                      cccAACACAAAGTATGGCATTGTACTGAAGAGATGATGGAACGGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGt        >  1:1515278/3‑100 (MQ=255)
                                                                                          acacAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc    <  1:1243916/100‑1 (MQ=255)
                                                                                          acacAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc    <  2:734922/100‑1 (MQ=255)
                                                                                          acacAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc    >  2:1215667/1‑100 (MQ=255)
                                                                                          acacAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc    <  1:1486816/100‑1 (MQ=255)
                                                                                            acaAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTg  >  1:1252598/1‑100 (MQ=255)
                                                                                            acaAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTg  >  1:1815132/1‑100 (MQ=255)
                                                                                                  |                                                                                             
AGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTG  >  NC_000913/3007089‑3007280

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: ‑