Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,684,604 G→A 100% Q339* (CAG→TAG)  yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,6040GA92.3% 262.5 / 10.0 93Q339* (CAG→TAG) yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base A (39/46);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACA  >  NC_000913/3684508‑3684699
                                                                                                |                                                                                               
gaTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAcc                                                                                              <  2:995593/100‑1 (MQ=255)
  tCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTa                                                                                            >  1:233634/1‑100 (MQ=255)
   cGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTaa                                                                                           >  1:899832/1‑100 (MQ=255)
   cGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTaa                                                                                           >  1:1344889/1‑100 (MQ=255)
     ccGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAgg                                                                                         <  1:236869/100‑1 (MQ=255)
      cGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGt                                                                                        <  2:1123146/100‑1 (MQ=255)
       gTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTg                                                                                       >  2:954640/1‑100 (MQ=255)
        tAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGa                                                                                      >  2:306206/1‑100 (MQ=255)
          gAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGAtt                                                                                    >  2:2231689/1‑100 (MQ=255)
           aaCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTg                                                                                   <  2:1974848/100‑1 (MQ=255)
             cATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCg                                                                                 <  1:1215555/100‑1 (MQ=255)
             cATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCg                                                                                 >  2:861195/1‑100 (MQ=255)
              aTCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGt                                                                                <  2:705863/100‑1 (MQ=255)
                 gccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCTCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGcc                                                                             <  2:2507386/100‑1 (MQ=255)
                 gccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGcc                                                                             <  2:2613873/100‑1 (MQ=255)
                   cacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAc                                                                           <  2:558163/100‑1 (MQ=255)
                     cgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACgg                                                                         <  1:40492/100‑1 (MQ=255)
                       cgaaTGCTACAGTGCTGACGGAGTTGAATACGTTCGATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACg                                                                          <  2:2348244/95‑1 (MQ=255)
                        caATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTc                                                                      >  1:572418/1‑100 (MQ=255)
                         aaTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTcc                                                                     >  2:2205198/1‑100 (MQ=255)
                           tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   <  2:242853/100‑1 (MQ=255)
                           tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   <  1:2279208/100‑1 (MQ=255)
                            gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  >  2:493350/1‑100 (MQ=255)
                            gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  <  2:970061/100‑1 (MQ=255)
                               aCAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTaca                                                               >  1:1717083/1‑100 (MQ=255)
                               aCAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTaca                                                               >  1:2539748/1‑100 (MQ=255)
                                  gTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCg                                                            >  2:2659343/1‑100 (MQ=255)
                                   tGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGt                                                           >  2:1063419/1‑100 (MQ=255)
                                    gCGTACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATATTGATCACTTGCTATTTTAAGGTGATTTCATGCCACGGTTCCTGTACACCGtt                                                          >  1:1929559/1‑100 (MQ=255)
                                    gCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGtt                                                          >  1:1127867/1‑100 (MQ=255)
                                     cggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTg                                                         <  1:2108552/100‑1 (MQ=255)
                                      gcactgtgTTGAATACGTTAAATCGGCAGGCGCTCGCGCATGATAGTGAGCACTCGCTAACCTAATGTGATTGCGTGCCACGGTTCCTGTACTCCGttn                                                         <  2:788323/92‑2 (MQ=255)
                                       gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                       <  1:841791/100‑1 (MQ=255)
                                       gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                       <  1:841660/100‑1 (MQ=255)
                                        acggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCa                                                      <  1:1100018/100‑1 (MQ=255)
                                         cggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCaa                                                     >  2:582146/1‑100 (MQ=255)
                                         cggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCaa                                                     >  1:580468/1‑100 (MQ=255)
                                           gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   >  1:1303247/1‑100 (MQ=255)
                                           gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   <  1:2422838/100‑1 (MQ=255)
                                           gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   <  2:2539730/100‑1 (MQ=255)
                                           gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   >  1:58723/1‑100 (MQ=255)
                                           gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   <  1:2182739/100‑1 (MQ=255)
                                             gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGAc                                                 <  2:1254048/100‑1 (MQ=255)
                                               tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                               <  1:2178728/100‑1 (MQ=255)
                                               tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                               >  2:2178713/1‑100 (MQ=255)
                                               tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                               >  1:442361/1‑100 (MQ=255)
                                               tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCAATTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                               >  2:1362352/1‑100 (MQ=255)
                                                gAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGc                                              <  1:1972420/100‑1 (MQ=255)
                                                  aTACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCaa                                            >  2:1906866/1‑100 (MQ=255)
                                                  aTACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCaa                                            <  2:745284/100‑1 (MQ=255)
                                                    aCGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAg                                          <  1:306758/100‑1 (MQ=255)
                                                    aCGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAg                                          <  1:1312413/100‑1 (MQ=255)
                                                     cGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGt                                         >  1:1790260/1‑100 (MQ=255)
                                                     cGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGt                                         <  2:99632/100‑1 (MQ=255)
                                                      gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                        >  2:1690837/1‑100 (MQ=255)
                                                      gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                        >  2:2401866/1‑100 (MQ=255)
                                                       ttCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTGAGGTGATTGCGTGCCACGGTGCCTGTACACCGTTGGCAATGACAGCAAAGTCa                                       <  2:968044/100‑1 (MQ=255)
                                                       ttCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCa                                       >  2:102297/1‑100 (MQ=255)
                                                       ttCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCa                                       <  2:1978967/100‑1 (MQ=255)
                                                        tCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCAt                                      >  1:907780/1‑100 (MQ=255)
                                                         cAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATa                                     >  2:553517/1‑100 (MQ=255)
                                                           aTCGGCAGGCGGTCGCTCATGATAGTGAGCACTTGCTAAGCTAAGGTGAGTGCGTGCCACGGTTGCTGTACACCGTTGGCAATGACAGCAAGGtaataa                                    >  2:1715470/1‑99 (MQ=255)
                                                           aTCGGCAGGCGCTCGCTCATGAGAGTGAGCATTTGCTTACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                                   <  2:1858806/100‑1 (MQ=255)
                                                              ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                >  1:1818860/1‑100 (MQ=255)
                                                              ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                >  2:446979/1‑100 (MQ=255)
                                                                  ggCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAt                            <  2:664654/100‑1 (MQ=255)
                                                                  ggCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAt                            >  1:1879124/1‑100 (MQ=255)
                                                                           tcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAg                   <  1:44850/100‑1 (MQ=255)
                                                                           tcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAg                   <  2:580467/100‑1 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  <  1:36303/100‑1 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  <  2:1818853/100‑1 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  >  1:506312/1‑100 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  >  1:2100572/1‑100 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  <  2:275888/100‑1 (MQ=255)
                                                                             aTGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGa                 >  2:1170104/1‑100 (MQ=255)
                                                                              tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt                <  1:36474/100‑1 (MQ=255)
                                                                              tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt                <  2:2690848/100‑1 (MQ=255)
                                                                              tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt                <  1:698037/100‑1 (MQ=255)
                                                                                  aGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATa            <  1:348690/100‑1 (MQ=255)
                                                                                  aGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATa            <  1:2428697/100‑1 (MQ=255)
                                                                                   gTGAGCACTTGCTGACCTAATGTGATAGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCAtgc           >  1:1597011/1‑98 (MQ=255)
                                                                                   gTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATAc           >  1:1625836/1‑100 (MQ=255)
                                                                                       gCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTgg       <  1:1204395/100‑1 (MQ=255)
                                                                                       gCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTgg       >  2:88270/1‑100 (MQ=255)
                                                                                        cACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGc      <  1:2430649/100‑1 (MQ=255)
                                                                                         aCTTGCTAACCTAAGGTGATTGCGTGCCACTGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCg     <  1:1183372/100‑1 (MQ=255)
                                                                                          cTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGa    >  1:441086/1‑100 (MQ=255)
                                                                                          cTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGa    >  1:569281/1‑100 (MQ=255)
                                                                                           ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  2:1783908/100‑1 (MQ=255)
                                                                                           ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  1:391832/100‑1 (MQ=255)
                                                                                           ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  2:139038/100‑1 (MQ=255)
                                                                                            tGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACa  >  1:528971/1‑100 (MQ=255)
                                                                                            tGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACa  >  1:694542/1‑100 (MQ=255)
                                                                                                |                                                                                               
GATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACA  >  NC_000913/3684508‑3684699

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