Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,778 T→A 5.4% M136L (ATG→TTG)  yjhP ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,7780TA 5.4% 244.4 / 3.8 74M136L (ATG→TTG) yjhPputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (37/33);  new base A (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAACCGGCCGCATCGTT  >  NC_000913/4532640‑4532907
                                                                                                                                          |                                                                                                                                 
ccccGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCGGGGCTATTTCTGCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGACGATAAGCCTGATcccgcc                                                                                                                          >  8:385031/4‑148 (MQ=255)
 cacGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATcccgccc                                                                                                                         >  5:101810/1‑148 (MQ=255)
    gTCATAGCCGAGATCGTCGAACGCTCCGAGAAGGCCTGGCAGCGTCTGGAAATCGCTTGGCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGt                                                                                                                      <  6:62747/148‑1 (MQ=255)
           ccGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGc                                                                                                               >  2:200975/1‑148 (MQ=255)
            cGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCt                                                                                                              >  7:5344/1‑148 (MQ=255)
               gaTCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGtt                                                                                                                     <  8:327787/138‑1 (MQ=255)
               gaTCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGtt                                                                                                                     >  7:327787/1‑138 (MQ=255)
                 tcgtcgAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTgcg                                                                                                         <  8:5344/148‑1 (MQ=255)
                     cgAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTag                                                                                                     >  5:146303/1‑148 (MQ=255)
                        aCGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTagcag                                                                                                  >  3:301133/1‑148 (MQ=255)
                              cGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTCGCTGACGCCAGTAGGGTGCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGctcc                                                                                              <  1:82569/146‑1 (MQ=255)
                               gACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTagca                                                                                                   <  4:271951/140‑1 (MQ=255)
                               gACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTagca                                                                                                   >  3:271951/1‑140 (MQ=255)
                               gACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCtcctcc                                                                                           <  4:301133/148‑1 (MQ=255)
                                aCAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCg                                                                                          >  4:291080/1‑148 (MQ=255)
                                           cAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAAc                                                                               <  6:146303/148‑1 (MQ=255)
                                            aGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAAcc                                                                              <  5:227169/148‑1 (MQ=255)
                                                 cAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCtt                                                                                                             >  1:116845/1‑112 (MQ=255)
                                                 cAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCtt                                                                                                             <  2:116845/112‑1 (MQ=255)
                                                 cAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCg                                                                         >  2:222485/1‑148 (MQ=255)
                                                   ggAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGc                                                                       <  1:200975/148‑1 (MQ=255)
                                                       aTCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCtcctcc                                                                                           <  2:48150/124‑1 (MQ=255)
                                                       aTCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCtcctcc                                                                                           >  1:48150/1‑124 (MQ=255)
                                                        tCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATcc                                                                  <  6:101810/148‑1 (MQ=255)
                                                                  aGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGcc                                                        >  8:281032/1‑148 (MQ=255)
                                                                   gCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCg                                                       <  3:291080/148‑1 (MQ=255)
                                                                       aCGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACg                                                    <  4:234909/147‑1 (MQ=255)
                                                                       aCGCCGCAGGCCGGGGCTATTTCTTCTGTCGCTGGTCGCTGACGCCAGTAGGGGTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCGTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCCCCGGCAATCCATGTCGCGCCGACGc                                                   <  7:385031/148‑1 (MQ=255)
                                                                        cgccgcAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTTGCGCCGACg                                                    >  3:336549/1‑146 (MQ=255)
                                                                           cgcCGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTcgcg                                                          >  3:69859/1‑137 (MQ=255)
                                                                           cgcAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTcgcg                                                          >  7:87738/1‑137 (MQ=255)
                                                                           cgcAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTcgcg                                                          <  8:87738/137‑1 (MQ=255)
                                                                           cgcAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTcgcg                                                          <  4:69859/137‑1 (MQ=255)
                                                                                    gggCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCACGCCGCCACGTc                                      >  6:2013/1‑148 (MQ=255)
                                                                                     ggCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCAAGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTca                                     >  5:340809/1‑148 (MQ=255)
                                                                                     ggCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCAAGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTca                                     >  5:335246/1‑148 (MQ=255)
                                                                                          tttCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACAttt                                >  7:364116/1‑148 (MQ=255)
                                                                                           ttcttcTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTagcag                                                                                                  <  4:82136/81‑1 (MQ=255)
                                                                                           ttcttcTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTagcag                                                                                                  >  3:82136/1‑81 (MQ=255)
                                                                                           ttcttcTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCAcc                                                                          >  1:109033/1‑105 (MQ=255)
                                                                                           ttcttcTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCAcc                                                                          <  2:109033/105‑1 (MQ=255)
                                                                                           ttcttcTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACAtttt                               >  4:246334/1‑148 (MQ=255)
                                                                                           ttcttcTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACAtttt                               >  8:296428/1‑148 (MQ=255)
                                                                                               tctGTTGATGGTAGCTGACGCCAGAAGGGTTCGCCGATAAGCATGATCCAGCCCGGATTCATGCTTTTCTCTAGCAGCTCCTCCGACCCGGCAAACCCACAGGCAATCCATGTCGCGCAGACGCAGCCCGCCACGTCACATTTTTCAt                           >  6:286014/1‑148 (MQ=255)
                                                                                                ctGTTGCTCGTAGCTGACGCCAGTAGGGTTCGCCGATATGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCAtt                          >  7:396576/1‑148 (MQ=255)
                                                                                                   ttGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCa                                                                      <  7:265890/101‑1 (MQ=255)
                                                                                                   ttGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCa                                                                      >  8:265890/1‑101 (MQ=255)
                                                                                                    tGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGc                       >  1:210023/1‑147 (MQ=255)
                                                                                                     gCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCg                      <  1:222485/147‑1 (MQ=255)
                                                                                                      cTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGa                     <  7:281032/147‑1 (MQ=255)
                                                                                                      cTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCAAGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGAc                    <  6:340809/148‑1 (MQ=255)
                                                                                                      cTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCAAGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGAc                    <  6:335246/148‑1 (MQ=255)
                                                                                                         gTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCACGCCGCCACGTCACATTTTTCATTTGCGACGTa                 <  5:2013/148‑1 (MQ=255)
                                                                                                          tAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTTGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTaa                <  4:336549/148‑1 (MQ=255)
                                                                                                          tAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTaa                <  1:141168/148‑1 (MQ=255)
                                                                                                           aGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAAc               >  3:146107/1‑148 (MQ=255)
                                                                                                            gCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACg                                                    <  2:393635/110‑1 (MQ=255)
                                                                                                            gCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACg                                                    >  1:393635/1‑110 (MQ=255)
                                                                                                                       tAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAACCGGCCGCATCGt   >  5:24777/1‑148 (MQ=255)
                                                                                                                        aGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACg                                                    >  4:308119/1‑98 (MQ=255)
                                                                                                                        aGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACg                                                    <  3:308119/98‑1 (MQ=255)
                                                                                                                        aGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTaa                >  7:73827/1‑134 (MQ=255)
                                                                                                                        aGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTaa                <  8:73827/134‑1 (MQ=255)
                                                                                                                        aGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTATCCGGCCGCATCGt   >  1:338874/1‑147 (MQ=255)
                                                                                                                        aGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAACCGGCCGCATCGt   >  7:49706/1‑147 (MQ=255)
                                                                                                                        aGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAACCGGCCGCATCGt   >  8:45671/1‑147 (MQ=255)
                                                                                                                         gggTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAACCGGCCGCATCGtt  <  8:364116/147‑1 (MQ=255)
                                                                                                                         gggTTAGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCTTCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAACCGGCCGCATCGtt  <  5:286014/147‑1 (MQ=255)
                                                                                                                             tCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCAtt                          <  6:70883/119‑1 (MQ=255)
                                                                                                                             tCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCAtt                          >  5:70883/1‑119 (MQ=255)
                                                                                                                                    tAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGc                                                   <  3:18366/87‑1 (MQ=255)
                                                                                                                                    tAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGc                                                   >  4:18366/1‑87 (MQ=255)
                                                                                                                                       gCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGctcctc                                                                                            <  7:127908/43‑1 (MQ=255)
                                                                                                                                       gCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGctcctc                                                                                            >  8:127908/1‑43 (MQ=255)
                                                                                                                                          |                                                                                                                                 
CCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGGCTTTGCGCTAGCAGCTCCTCCGCCCCGGCAAACCCACCGGCAATCCATGTCGCGCCGACGCAGGCCGCCACGTCACATTTTTCATTTGCGACGTAACCGGCCGCATCGTT  >  NC_000913/4532640‑4532907

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