Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,389 A→G 6.1% F55F (TTT→TTC)  ybcI ← DUF457 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,3890AG 6.1% 203.9 / 3.1 66F55F (TTT→TTC) ybcIDUF457 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (34/32)
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.17e-01

AGCTGTGCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCG  >  NC_000913/556250‑556537
                                                                                                                                           |                                                                                                                                                    
aGCTGTGCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTg                                                                                                                                             <  6:347164/149‑1 (MQ=255)
 gCTGTGCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGc                                                                                                                                            <  1:145635/149‑1 (MQ=255)
 gCTGTGCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGc                                                                                                                                            <  5:324327/149‑1 (MQ=255)
  cTGTGCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGcc                                                                                                                                           <  2:76913/149‑1 (MQ=255)
     tgCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGCCCAAAAACATTGCCGt                                                                                                                                         >  3:66273/1‑148 (MQ=255)
          tAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCaa                                                                                                                                   >  5:120998/1‑149 (MQ=255)
            gCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAAc                                                                                                                                  <  3:129311/148‑1 (MQ=255)
                 gagaCGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCaaa                                                                                                                            <  7:315327/149‑1 (MQ=255)
                           aaaaCAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGa                                                                                                                   <  1:123999/148‑1 (MQ=255)
                                agccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACACCGGGACAACAATCGCAAACACCAGCGAATGGGTAATCCCGCGATGACCACAAACATTGCCGTAAGCAACGCCAAATTTAATCGACaaatc                                                                                                             >  1:36896/1‑146 (MQ=255)
                                   cagcagCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATATCAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATAcgtc                                                                                                          <  4:242518/149‑1 (MQ=255)
                                    agcagcGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCCCGACAACAAACGCAAACACCACCGAATGGGTAAACCCCCGATGACCCAAAACATTGCCGTAAGCCACGCCAAATTTAAACGACAAGAcgtcg                                                                                                         >  5:32083/1‑149 (MQ=255)
                                             tCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGgc                                                                                                <  1:2382/149‑1 (MQ=255)
                                                       ccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGc                                                                                        >  3:180148/1‑147 (MQ=255)
                                                        cTGAACCATCGGCGGCCAATTAACACACATACCAGCGGGACAACAACCGCAAACACCAGGGAATGGGTAAACCCGCGTTGACCAACACCATTGCCGTAAGCAACGCCAAATTTACAGGACAATACGTCGGCGTCGGGCTGCATCGCCAg                                                                                     <  5:307885/149‑1 (MQ=255)
                                                           aaCCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGgat                                                                                  <  6:120998/149‑1 (MQ=255)
                                                           aaCCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGGCGGGCAGGATGGCCAGga                                                                                   <  6:107365/148‑1 (MQ=255)
                                                              cATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGgatgat                                                                               >  1:258026/1‑149 (MQ=255)
                                                                tCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGGAACCCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGGCGGGCAGCAGCGCCAGGATGATTc                                                                             <  1:87123/149‑1 (MQ=255)
                                                                   gCGTCCAATAAACACACATAACAGCGGGACAACAAACGAAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGCCGTCGGGCAGACTCGCCAGGATGATTCCgg                                                                          >  3:355957/1‑149 (MQ=255)
                                                                    cGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGc                                                                         <  4:180148/149‑1 (MQ=255)
                                                                          aTAAACACACATAACAGCGGGACAACAAAGGCAAACACCAGCGATTGGGTAAACCCGCGATGACCAAAAAGATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATaa                                                                   <  6:32083/149‑1 (MQ=255)
                                                                              acacacATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCGGGATGATTCCGGCAAATAACAGa                                                               >  4:262837/1‑149 (MQ=255)
                                                                                acaaaTAACAGCGGGACAACAAACGCAAACACCAGCGAATGTGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATAAGTCGGCGTCGGGCAGCATCGCCAGGATGAGTCCCGCAAAAAACAGCcg                                                             >  8:276898/1‑149 (MQ=255)
                                                                                 cacaTAACACCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGAcgc                                                            <  6:98633/149‑1 (MQ=255)
                                                                                  acaTAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGAcgcg                                                           >  1:381905/1‑149 (MQ=255)
                                                                                    aTAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGAtt                                                                              <  7:376619/128‑1 (MQ=255)
                                                                                    aTAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGAtt                                                                              >  8:376619/1‑128 (MQ=255)
                                                                                    aTAACAGCGGGACAACAAAAGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGGCGGCGTCGGGCAGCATCGCCAGGAAGATTCCGGCAAATAACAGACGAGGa                                                         >  2:252741/1‑149 (MQ=255)
                                                                                      aaCAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAgg                                                       >  6:44588/1‑149 (MQ=255)
                                                                                        cAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGa                                                     >  1:274013/1‑149 (MQ=255)
                                                                                        cAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGa                                                     <  2:258026/149‑1 (MQ=255)
                                                                                            gggACAACAATAGCAAACACCAGAGAATGGGTAAACCCGCGATGACCAATAACATTGCCGTATGCAAAGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAATTAACAGACGCGGAGGGATGAc                                                 <  2:36896/149‑1 (MQ=255)
                                                                                             ggACACCAAACGCAAACACCAGCGAATGGGTCAACCCGCGATGACCAAAAACATTGCCGTAAGCACCGCCAACTTTAAACGACAATACGTCGGCGGCGGGAAGCATCGCCAGGATGAGTCCGGCAAATAACAGACGCGGAGGGATGCCt                                                >  5:200553/1‑149 (MQ=255)
                                                                                              gACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATa                                                                    >  6:368720/1‑128 (MQ=255)
                                                                                              gACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATa                                                                    <  5:368720/128‑1 (MQ=255)
                                                                                                caacaaACGCAAACACCAGCGAATGGGTAAACCAGCGATGACCATAAACATTGCCGTAATCAAAGCCAAATTTAAACGACAATACGGCGGCTTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGTATGACTTTc                                             <  1:252741/149‑1 (MQ=255)
                                                                                                aaacacACGCAAACACCCGCGAATGGGTACTCCCTCGCTGACCAAATACCTTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGAGTCCGGCAAATAACAGACGCGGAGGGATGACTTTc                                             <  6:200553/148‑1 (MQ=255)
                                                                                                  acaaACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCCTCATCGCCAGGATGATGCCGGCAAATAACAGACGCTGAGGGATGACTTTCg                                            >  2:36856/1‑148 (MQ=255)
                                                                                                      aCGCAAACACCAGCGAATGGGTAAACCCGCGATGACCGAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaaccc                                       >  2:18830/1‑149 (MQ=255)
                                                                                                      aCGCAAACACCAGCGAATGGGTAAACCCGCGATGACCGAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaacc                                        >  2:12777/1‑148 (MQ=255)
                                                                                                       cGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaaccc                                       >  6:44803/1‑148 (MQ=255)
                                                                                                        gCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGga                                                                                   >  4:201302/1‑103 (MQ=255)
                                                                                                        gCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGga                                                                                   <  3:201302/103‑1 (MQ=255)
                                                                                                        gCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGGCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaacccaa                                     <  3:202916/149‑1 (MQ=255)
                                                                                                          aaaCACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaacccaacc                                   >  1:304408/1‑149 (MQ=255)
                                                                                                          aaaCACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaacccaacc                                   >  1:297704/1‑149 (MQ=255)
                                                                                                           aaCACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaacccaaccc                                  >  3:330719/1‑149 (MQ=255)
                                                                                                               ccAGCGAATGGGTAAACCCGCGATGACCGAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTaaa                              <  1:18830/149‑1 (MQ=255)
                                                                                                               ccAGCGAATGGGTAAACCCGCGATGACCGAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTaaa                              <  1:12777/149‑1 (MQ=255)
                                                                                                                cAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAAc                             >  1:382601/1‑149 (MQ=255)
                                                                                                                cAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAAc                             <  5:44803/149‑1 (MQ=255)
                                                                                                                cAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAAc                             <  2:381905/149‑1 (MQ=255)
                                                                                                                      aTGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCGGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGc                       <  3:262837/149‑1 (MQ=255)
                                                                                                                          gTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAg                   >  1:274531/1‑149 (MQ=255)
                                                                                                                          gTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAg                   >  1:11348/1‑149 (MQ=255)
                                                                                                                          gTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAg                   >  1:82651/1‑149 (MQ=255)
                                                                                                                          gTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAg                   <  2:274531/149‑1 (MQ=255)
                                                                                                                                ccGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAAc             >  8:305571/1‑149 (MQ=255)
                                                                                                                                ccGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGCATGATTCCGGCAAATAACAGACGCGGAGGGATGACATACGAACCCAACCCTAAACCAATGCATAGGTGAAc             >  6:239520/1‑149 (MQ=255)
                                                                                                                                        aCCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGgcgt     <  2:274013/149‑1 (MQ=255)
                                                                                                                                        aCCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGAATCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGgcgt     <  1:36856/149‑1 (MQ=255)
                                                                                                                                         ccAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGCACCCAACCCTAACCCAATGCATAGGGGAAAGgcggcgag    >  2:165081/1‑147 (MQ=255)
                                                                                                                                         ccAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGgcgtg    >  4:293274/1‑149 (MQ=255)
                                                                                                                                          cAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGCACCCAACCCTAAACCAATGCATAGGGGAACGGCGgcgtgc   >  3:113882/1‑149 (MQ=255)
                                                                                                                                          cAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGCCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTACGACCCCAACCCTAAACCAAAGCATAGGGGAACGGCGgcgtgc   >  6:241226/1‑149 (MQ=255)
                                                                                                                                           aaaaaCATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGgcgtgcg  >  7:317580/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
AGCTGTGCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCG  >  NC_000913/556250‑556537

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