Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,006,046 C→G 5.0% S32* (TCA→TGA)  zapC → FtsZ stabilizer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,006,0460CG 5.0% 269.4 / 5.6 80S32* (TCA→TGA) zapCFtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base C (34/42);  new base G (2/2);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

GAAGCGACAATGGATTGTGTTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATG  >  NC_000913/1005900‑1006192
                                                                                                                                                  |                                                                                                                                                  
gaagcgaCAATGGATTGTGTTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCAc                                                                                                                                                  <  3:304815/149‑1 (MQ=255)
gaagcgaCAATGGATTGTGGTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCAc                                                                                                                                                  <  4:425506/149‑1 (MQ=255)
       cAATGGATTGTGTTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTgc                                                                                                                                           <  3:303676/149‑1 (MQ=255)
              ttGTGTTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTgcgcgc                                                                                                                                       >  2:66777/1‑146 (MQ=255)
                gtgtTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGa                                                                                                                                   <  7:503692/148‑1 (MQ=255)
                 tgtTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATg                                                                                                                                 <  1:112524/149‑1 (MQ=255)
                   tttctgcGGTTTTATCGGGTAGGGGAGAGGCATCTGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGtt                                                                                                                               <  6:324815/146‑1 (MQ=255)
                    tgctgcGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGtt                                                                                                                               <  1:193726/148‑1 (MQ=255)
                    tgctgaGGTTTTATAGGGTAGGGGAGAGGCTGATGCGAATTAAACCAGACGATATCTGGCGTTGGTTTTACGATGAAGAGCTCGATCGAATGATGGTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGGGCAAGATGtt                                                                                                                               <  4:186266/148‑1 (MQ=255)
                       tgcGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACt                                                                                                                           >  6:111887/1‑149 (MQ=255)
                          ggTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCa                                                                                                                        <  5:434091/149‑1 (MQ=255)
                            ttttATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGa                                                                                                                      <  3:423567/149‑1 (MQ=255)
                                   gggtagggGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCttt                                                                                                               >  1:185547/1‑149 (MQ=255)
                                   gggtagggGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCtt                                                                                                                <  4:461283/148‑1 (MQ=255)
                                    gggagggGCGCGGCCGCTTCGTATTAAACCCGCCCCAAACCGGCGGTGGGATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCtttt                                                                                                              <  6:192817/146‑1 (MQ=255)
                                     gtagggGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTGACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTtc                                                                                                             >  4:262614/1‑149 (MQ=255)
                                     gtagggGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTGACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTtc                                                                                                             >  4:268799/1‑149 (MQ=255)
                                     gtagggGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTtc                                                                                                             <  5:111887/149‑1 (MQ=255)
                                       agggGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTtctc                                                                                                           >  8:262791/1‑149 (MQ=255)
                                        ggggAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTtctc                                                                                                           <  2:95048/148‑1 (MQ=255)
                                          ggAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCg                                                                                                        >  2:405806/1‑149 (MQ=255)
                                          ggAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCg                                                                                                        <  6:284975/149‑1 (MQ=255)
                                            agagGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCCCGATCGAATGATGCTCGATTTATCCAATGGTATGCTATTTCGCTCACGATTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCa                                                                                                      <  6:466840/149‑1 (MQ=255)
                                            agagGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCTAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGc                                                                                                                   >  3:119560/1‑136 (MQ=255)
                                            agagGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGc                                                                                                                   <  4:119560/136‑1 (MQ=255)
                                                gCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCt                                                                                                  <  1:405806/149‑1 (MQ=255)
                                                  aGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCttt                                                                                                <  8:361457/149‑1 (MQ=255)
                                                   gATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTGACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCtttt                                                                                               <  3:268799/149‑1 (MQ=255)
                                                   gATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTGACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCtttt                                                                                               <  3:262614/149‑1 (MQ=255)
                                                   gATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCtttt                                                                                               >  5:59745/1‑149 (MQ=255)
                                                        gAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGtt                                                                                          >  7:454679/1‑149 (MQ=255)
                                                              aaCCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTgacgac                                                                                    >  7:475443/1‑149 (MQ=255)
                                                               aCCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCaa                                                                                                                                     >  6:172382/1‑99 (MQ=255)
                                                               aCCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCaa                                                                                                                                     <  5:172382/99‑1 (MQ=255)
                                                                   gACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGATGAcgccgc                                                                               >  2:476016/1‑149 (MQ=255)
                                                                       aTAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGctc                                                                           >  4:29207/1‑149 (MQ=255)
                                                                            tGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGCCGCCGCGCTCTAttt                                                                      >  2:131374/1‑149 (MQ=255)
                                                                                    gTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGtt                                                                                                                               <  1:27466/84‑1 (MQ=255)
                                                                                    gTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGtt                                                                                                                               >  2:27466/1‑84 (MQ=255)
                                                                                     tATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTg                                                              >  7:408374/1‑148 (MQ=255)
                                                                                      aTTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTg                                                              >  3:363121/1‑147 (MQ=255)
                                                                                       ttACGATTAAGAGCACGATCGAATGATGCTCGATTTAGCGAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTgaa                                                            <  8:454679/148‑1 (MQ=255)
                                                                                       ttACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTgaa                                                            <  8:408374/148‑1 (MQ=255)
                                                                                         aCGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTgaaga                                                          >  6:381294/1‑148 (MQ=255)
                                                                                           gATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTgaagaa                                                         >  4:244640/1‑147 (MQ=255)
                                                                                             tGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGt                                                     <  8:282581/149‑1 (MQ=255)
                                                                                             tGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGt                                                     <  8:275175/149‑1 (MQ=255)
                                                                                                 gagCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGctct                                                                          >  7:442864/1‑124 (MQ=255)
                                                                                                 gagCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGctct                                                                          <  8:442864/124‑1 (MQ=255)
                                                                                                  agCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCcgc                                                <  3:29207/149‑1 (MQ=255)
                                                                                                    cACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTg                                                              <  7:81601/133‑1 (MQ=255)
                                                                                                    cACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTg                                                              >  8:81601/1‑133 (MQ=255)
                                                                                                    cACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGa                                              >  6:277097/1‑149 (MQ=255)
                                                                                                     aCGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGCTGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGAt                                             >  3:34286/1‑149 (MQ=255)
                                                                                                          cGAATGATACTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTaa                                        >  7:262880/1‑149 (MQ=255)
                                                                                                            aaTGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAAtt                                      >  8:113032/1‑149 (MQ=255)
                                                                                                               gatgCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTAt                                   >  3:189663/1‑149 (MQ=255)
                                                                                                                atgCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTAt                                   <  4:197104/148‑1 (MQ=255)
                                                                                                                atgCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATc                                  <  1:66777/149‑1 (MQ=255)
                                                                                                                       aTTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGAtt                                            <  7:74773/132‑1 (MQ=255)
                                                                                                                       aTTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGAtt                                            >  8:74773/1‑132 (MQ=255)
                                                                                                                             ccAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACaaaaa                     <  5:381294/149‑1 (MQ=255)
                                                                                                                               aaTGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGATGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGc                   <  1:476016/149‑1 (MQ=255)
                                                                                                                               aaTGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTAt                                                                        <  2:489874/96‑1 (MQ=255)
                                                                                                                               aaTGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTAt                                                                        >  1:489874/1‑96 (MQ=255)
                                                                                                                               aaTGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGaaaa                                                       >  4:296245/1‑113 (MQ=255)
                                                                                                                               aaTGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGaaaa                                                       <  3:296245/113‑1 (MQ=255)
                                                                                                                                   gTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGaa               <  8:262880/149‑1 (MQ=255)
                                                                                                                                        cTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTgg           <  5:277097/148‑1 (MQ=255)
                                                                                                                                        cTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGACTTTTCTCCCGCAGGCTTTTGCTTTGACGACGCCGCGCTCTATTTATCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGt          >  4:300413/1‑149 (MQ=255)
                                                                                                                                         tATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCATATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCtt                                                                 <  3:174389/93‑1 (MQ=255)
                                                                                                                                         tATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCtt                                                                 >  4:174389/1‑93 (MQ=255)
                                                                                                                                         tATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTg         >  6:190342/1‑149 (MQ=255)
                                                                                                                                          aTTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGc        >  4:429158/1‑149 (MQ=255)
                                                                                                                                           tttCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCt       >  1:213354/1‑149 (MQ=255)
                                                                                                                                           tttCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCt       <  3:244640/149‑1 (MQ=255)
                                                                                                                                           tttCGCTCACGGTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAACAGCCGAACTGGTg         <  8:475443/147‑1 (MQ=255)
                                                                                                                                               gCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGaa    <  6:59745/148‑1 (MQ=255)
                                                                                                                                                cTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATg  <  1:19978/149‑1 (MQ=255)
                                                                                                                                                 tCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATg  >  5:165655/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAAGCGACAATGGATTGTGTTGCTGCGGTTTTATAGGGTAGGGGAGAGGCAGATGCGAATTAAACCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATG  >  NC_000913/1005900‑1006192

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