Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,294 G→C 15.7% D292H (GAT→CAT)  ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,2940GC15.7% 409.0 / 63.1 159D292H (GAT→CAT) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (65/68);  new base C (9/16);  total (74/85)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-02

Only 100 of 166 total aligned reads displayed.

CGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181196‑3181393
                                                                                                  |                                                                                                   
cGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatga                                                                                                    >  1:1465372/1‑100 (MQ=255)
 gTGGAAGAGCATTATCTCCAACAAGGTACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   <  1:3091097/100‑1 (MQ=255)
 gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   <  2:875230/100‑1 (MQ=255)
   ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                                 <  1:2226882/100‑1 (MQ=255)
      agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                              <  1:3521266/100‑1 (MQ=255)
       gagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCGCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCg                                                                                             <  1:2266224/100‑1 (MQ=255)
        agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                            >  2:2413576/1‑100 (MQ=255)
        agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCAAAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                            >  1:2141955/1‑100 (MQ=255)
         gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCa                                                                                           <  1:1365050/100‑1 (MQ=255)
           aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACGTGCTGCGCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                         >  2:2753019/1‑100 (MQ=255)
           aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaa                                                                                         <  1:2703536/100‑1 (MQ=255)
            ttATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTTCGGATGATCATCTTCCGCaaaa                                                                                        <  2:2267974/100‑2 (MQ=255)
              aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                      <  1:10454/100‑4 (MQ=38)
               tctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataa                                                                                     <  1:1038749/100‑5 (MQ=38)
               tctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGa                                                                                     >  1:1109156/1‑100 (MQ=255)
                ctcCAACAAGGCACTTCTATCGCTACTGTGGGAGATGTTCCAGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataag                                                                                    <  1:2094486/100‑6 (MQ=37)
                ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaa                                                                                    >  1:2770867/1‑100 (MQ=255)
                  ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                                  <  2:850024/100‑8 (MQ=37)
                   caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcgg                                                                                 >  2:3232766/1‑92 (MQ=37)
                    aacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcggg                                                                                <  2:3765576/100‑10 (MQ=37)
                     acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCTAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                               <  2:1423731/100‑1 (MQ=255)
                     acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                               <  2:3648224/100‑1 (MQ=255)
                      caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACTCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAaaat                                                                              <  2:665609/100‑1 (MQ=255)
                       aaGGCACGTGTACCGCTACTGTGGGAGATGTGCCTGAATCACCCGGACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcgg                                                                                 <  2:693678/96‑9 (MQ=21)
                        aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                            >  1:2590057/1‑100 (MQ=255)
                           naCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                         >  2:1329203/2‑100 (MQ=255)
                           cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTgg                                                                          <  1:585509/99‑1 (MQ=255)
                           cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                         <  1:801982/100‑1 (MQ=255)
                            aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                        <  1:398688/100‑1 (MQ=255)
                             cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                        <  1:740424/99‑1 (MQ=255)
                             cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTa                                                                       >  1:888925/1‑100 (MQ=255)
                              ttCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTgg                                                                          <  1:585928/96‑1 (MQ=255)
                                 tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   >  1:3173023/1‑100 (MQ=255)
                                 tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   <  2:2706731/100‑1 (MQ=255)
                                 tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   <  1:1329217/100‑1 (MQ=255)
                                  aCCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCg                                                                  <  1:449039/100‑1 (MQ=255)
                                   ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                 <  1:1780586/100‑1 (MQ=255)
                                   ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                 >  2:3436087/1‑100 (MQ=255)
                                   ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                 >  2:3586752/1‑100 (MQ=255)
                                    cGCTACTGTGTGAGATGTTCCCGAATTACCCGAAGCTGCTGCACGAGTATTTTGCGGAAGATCTTCATAAGCAAATGGAAAAATATGTGGTTAAAcc                                                                   <  2:3226642/97‑1 (MQ=255)
                                     gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                               <  1:2442229/100‑1 (MQ=255)
                                     gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                               <  2:294252/100‑1 (MQ=255)
                                     gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                               <  2:34094/100‑1 (MQ=255)
                                      cTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                               >  1:1326487/1‑99 (MQ=255)
                                         cTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtcttct                                                           <  1:1238059/100‑1 (MQ=255)
                                          tgtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTtctc                                                          >  1:3380419/1‑100 (MQ=255)
                                           gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                         <  1:2947247/100‑1 (MQ=255)
                                           gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                         >  1:911850/1‑100 (MQ=255)
                                             gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                       <  1:329354/100‑1 (MQ=255)
                                             gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                       >  1:1520593/1‑100 (MQ=255)
                                              gtaaaTGTTCCCGAATCAACCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                       <  2:343418/95‑1 (MQ=255)
                                               gagaTGTTCCCTAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGTAAAAATATGTGGTTAAACCGATCTTCTCCCGt                                                      <  2:1914643/99‑1 (MQ=255)
                                               gagaTGTGCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCAAAAAAAAAAAAAAAAATATGTGGTTAAACCGATCTTCTcc                                                         <  2:900469/96‑1 (MQ=255)
                                                   tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                 >  2:2562762/1‑100 (MQ=255)
                                                    gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCAGCAAAAGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                                >  1:17983/1‑100 (MQ=255)
                                                     ttCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                               >  1:2223224/1‑100 (MQ=255)
                                                     ttCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                               <  1:2028640/100‑1 (MQ=255)
                                                     ttCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                                >  1:3144170/1‑99 (MQ=255)
                                                      tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGACGATGATCATCCACAAATAGAAAAATATAAGGTTAAACCGATCTTCTCCCGTGAAGgcac                                              >  1:2330633/1‑98 (MQ=255)
                                                      tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                              >  1:2618784/1‑100 (MQ=255)
                                                      tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                              >  2:1118454/1‑100 (MQ=255)
                                                      tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                              >  1:2943231/1‑100 (MQ=255)
                                                        ccGCATCACCCGAACCTGCTGCTCGCTTATTTTGCGGAAGGTGATCATCCGCAAAGGGAAAAATATGTGCGTCAACCGATCTTCTCACGTGAAGGCGCaa                                            >  2:56015/1‑100 (MQ=255)
                                                        ccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                             <  1:333244/99‑1 (MQ=255)
                                                          gAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAAc                                          <  2:2946210/100‑1 (MQ=255)
                                                             tCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtg                                       <  1:3468973/100‑1 (MQ=255)
                                                             tCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtg                                       <  1:2494669/100‑1 (MQ=255)
                                                              cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                      >  2:1678810/1‑100 (MQ=255)
                                                               aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     >  2:3856088/1‑100 (MQ=255)
                                                               aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     >  1:3706014/1‑100 (MQ=255)
                                                               aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     <  1:3277430/100‑1 (MQ=255)
                                                               aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     >  2:2863223/1‑100 (MQ=255)
                                                                cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  1:3910798/1‑100 (MQ=255)
                                                                cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  1:3612417/1‑100 (MQ=255)
                                                                cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  1:2935738/1‑100 (MQ=255)
                                                                 ccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGa                                   >  1:1686341/1‑100 (MQ=255)
                                                                  cGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGTAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGAt                                  <  1:171720/100‑1 (MQ=255)
                                                                   gAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATc                                 <  2:2609030/100‑1 (MQ=255)
                                                                      ccTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCACAAAAAAAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAATCGTGTCGATCAtt                              <  2:2408925/100‑1 (MQ=255)
                                                                       ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  1:2853526/100‑1 (MQ=255)
                                                                       ctgctgCCCGCTTATTTTGCAGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  1:459981/100‑1 (MQ=255)
                                                                        tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                            <  1:1620308/100‑1 (MQ=255)
                                                                        tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCAATCTTCTCACGTGAAGGCGCAAACGTGTCGATCATTga                            >  1:1462600/1‑100 (MQ=255)
                                                                          ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgaga                          <  2:250793/100‑1 (MQ=255)
                                                                            gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                        >  1:107662/1‑100 (MQ=255)
                                                                            gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                        >  2:1658817/1‑100 (MQ=255)
                                                                             cccGCTTATTCTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCTGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                       <  2:2572623/100‑1 (MQ=255)
                                                                              ccGCTTATTTTGCGGAAGATGATTATCCGTAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACgg                      <  1:79556/100‑1 (MQ=255)
                                                                              ccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACgg                      <  1:2894901/100‑1 (MQ=255)
                                                                                   tATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                 >  2:3062292/1‑100 (MQ=255)
                                                                                     atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTAAGAACGGCAAAAcc               >  2:3495579/2‑100 (MQ=39)
                                                                                       ttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAt             <  2:1827612/100‑1 (MQ=255)
                                                                                       ttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAt             >  2:1017531/1‑100 (MQ=255)
                                                                                       ttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAt             <  1:2975453/100‑1 (MQ=255)
                                                                                        tGCGGTAGATGATCATCCGCAAATGGAAAAATATGTGGATAAACCGATGTTCTCCCGTGAAGGCGCAAACGTGGCGATCATTGAGAACGGCAAAACCAt             >  2:488567/1‑99 (MQ=255)
                                                                                        tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            <  1:1877307/100‑1 (MQ=255)
                                                                                        tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            <  2:3910759/100‑1 (MQ=255)
                                                                                                atgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagcagc    >  1:2798889/1‑100 (MQ=255)
                                                                                                 tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg   <  2:1589063/100‑1 (MQ=255)
                                                                                                  gatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCgg  >  2:2023283/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181196‑3181393

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: