Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,306,035 C→T G188D (GGC→GAC)  yjcS ← linear primary‑alkylsulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,0350CT94.3% 196.3 / 0.7 73G188D (GGC→GAC) yjcSlinear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (38/28);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACT  >  NC_000913/4305791‑4306272
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ggcgatatgtttggggggTGCAATAAGGTTGGGTTCGCCCGTTGCCAGCGTCACTCCAAGACCATTCCCCCCTTTGCCTGGCGGTTTGGGGGGCGGAAACAGACCGTAAGAGTGGCGCGCAGTGCGCCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGtcac                                                                                                                                                                                                                                             <  2:658137/232‑1 (MQ=255)
   gATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAgg                                                                                                                                                                                                                                       >  1:661379/1‑250 (MQ=255)
       gTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTc                                                                                                                                                                                                                                   >  1:500333/1‑251 (MQ=255)
        ttttCGACGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGCTAATGCCTTTCACGTCACCATAGTTGTCg                                                                                                                                                                                                                                  >  2:283209/1‑251 (MQ=255)
            cgtcgATGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCTCATTGCCTTGCACGTTGCGCTGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACATTTTCGCTGATGGCTTCGGACCTCAAGCCTGCAGGGGAAATCCCCTGAACTTTGCCGGATTTACCCTCTGCTTCACGGTGAATGCCTTTCACGTCACCATaagggacggggt                                                                                                                                                                                                                              >  2:514004/1‑240 (MQ=255)
               cgGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGc                                                                                                                                                                                                                           >  1:162613/1‑251 (MQ=255)
                ggTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTATCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCt                                                                                                                                                                                                                          <  1:213522/251‑1 (MQ=255)
                   gaaataatGCTCGGGCCGCCCTTTGCCAGCGCCACGCCAAGGCCTTTCCCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGTTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCtgtg                                                                                                                                                                                                                       <  1:233732/249‑1 (MQ=255)
                            cTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTgtgtgtgt                                                                                                                                                                                                                  >  1:381815/1‑247 (MQ=255)
                                ggTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACg                                                                                                                                                                                                          >  1:678145/1‑251 (MQ=255)
                                ggTCGCCCGTTGCCAGCGTCAAGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCATGCTGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATACAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTACACGTCACCATAGTGGTCGGTGTGGCTGTGAgttgagtttgcg                                                                                                                                                                                                          >  2:295604/1‑241 (MQ=255)
                                    gCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGACGGTGTGGCtggggggggagttaacggca                                                                                                                                                                                                       >  2:55213/1‑232 (MQ=255)
                                     cccGTTGCCAGGTCCCCGCCAGGGCCATTACCCACATTGCCTGGCGCGTTGGGCGGCAGAAACAGACCGTAAGAGTAGAGCGCACGGGGGCTCTTGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGcaac                                                                                                                                                                                                     <  1:463288/251‑1 (MQ=255)
                                         ttGCCAGAGCCACGTCAAGGCCATTACCCACATTGCCTTGCGCATTGTGCGGCAGTAACCGACCGTAAGAGCAGATCGCACGGCGCCTCATGATGTTCCCCGACAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGACTGCAGGGGCAATCACCTGAACTTTGCAGGATTTAACATCGGCCTCACAGATAATGCCTTTACCGTCACCATCGTGGTCGGTGTGGCTGTGAGAGTAGTTAACGGcacccca                                                                                                                                                                                                  >  2:81237/1‑245 (MQ=255)
                                             cAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTACCAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGGGGGAGTAGGTAACCGCAAAAAACGGtt                                                                                                                                                                                             >  2:332364/1‑251 (MQ=255)
                                             cAGCGACACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCGGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGAGAAGGCCTTTCACGTAACCATAGTGGTCGGTGTGGCTGGGGATGTAAATAACGGCAAAAATCGGtt                                                                                                                                                                                             >  2:152768/1‑251 (MQ=255)
                                               gCGTCACTCCAAGGCCTTTACCCACATTGCCGTGCGCGTTGGGCGGCAGTAACAGACCGTAAGGGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCAGAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGtttt                                                                                                                                                                                           <  1:567662/251‑1 (MQ=255)
                                                   cACGCCAATGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTCACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGTGGCAATCACCTGCACTTTGCCGGATTTAACATCGGCTTCTGGGGTCATTCCTTTCACGTCACCATAGTGGTCGGTGAGGCTGGGTGTTTAGATAACGGCAACCATCCGTTTTTg                                                                                                                                                                                         >  2:20304/1‑249 (MQ=255)
                                                          aaGGCCATTACACACATTGCCTTGCGCGTTGTGCGGCAGTAACAGGACGAAAAAGTAGAGCGCAAGGCGGCTCATGATGTTACACGCAAGCACGTTTTCGATGATTGCTGCGTCCATAAAGCCTGCGGGGTCAATCACCTGTAATTTGCCGGATTTAACATCTGCTTCGGAGATAATGCCTTTCCCGTCGACATAGGGGTCGGTGTGGCTGAGAGTAGAGAAAGCGGCATCAATCGGgtttgtcaggcgat                                                                                                                                                                                 >  2:615995/1‑241 (MQ=255)
                                                           aGGCCATTACCCACATTGCCTTGCGCGTTGGGCGGCAGTAACAGACCGTAAGGGGAGAGCGCCCGGCGTCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATg                                                                                                                                                                                <  1:360453/251‑1 (MQ=255)
                                                           aGGCCATTACCAACGTGGCCTGGCGCGTGGTGCGGCAGTAACAGACCGTAAGAGTATATCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCGGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATg                                                                                                                                                                                <  1:615995/251‑1 (MQ=255)
                                                              ccATTACCCACATTGCCTTGCGCGTTGGGCGGAAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTg                                                                                                                                                                             <  1:168100/251‑1 (MQ=255)
                                                                 ttACCCACATTGCCTTGCGCGTTGAGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCCCCTGAAATTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCAAGTCACCATAGTGGTCGGCGTGGCTGTGAGTGTAGATAAAGGCACCAATCGGTTTTTGCGGACGATGCTGGa                                                                                                                                                                           >  2:323196/1‑250 (MQ=255)
                                                                   aCCCACATTGCCGTGCGGGTTGTGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGTGGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAg                                                                                                                                                                         <  1:580231/250‑1 (MQ=255)
                                                                       ataTTGCCTTGCGCGTTGGGCGGCTTTAACAGCCCAAAAGAATAGAGCGCTCGGCGGTCCATGGGTTATCCCGCCAGCACGTTTTCGCTGATGGCTTCTTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGGCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAg                                                                                                                                                                    <  2:549239/249‑1 (MQ=255)
                                                                        caTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                                                                   >  1:239761/1‑251 (MQ=255)
                                                                             ccTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCa                                                                                                                                                                >  1:113493/1‑249 (MQ=255)
                                                                                             gCAGTAACAGACCGTAAGCGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATCAAGCCTGCAGGGGCAATCCCCTGCACTTTGCCGGTTTTAACATCGGCTTCAGAGATAATGCCCTTCACGTCACCATAGGGGTCGGTGTGGCTGTGAGTGTCGATAACGGCACCAATCGGTTTTTGCGGACGATGCTGGACGTACAGGTCAAGAGCGGCTTTTGCgtc                                                                                                                                               >  1:431047/1‑248 (MQ=255)
                                                                                               atgatcagacagaaggtggggagCGCCCGGTGGCTCGTGTGTTCACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGACTTAACATCGGCTTCAGAGATAATGCCTTTCGCGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAgg                                                                                                                                            <  2:195812/231‑1 (MQ=255)
                                                                                                     aGACCGTCAGAGTAGAGCGCCCGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGATTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGTTAATGCCTTTAACGTCACCATCGTGGTCGGAGTGGCGGTGAGTGTAAGTAACGGGAACCATCGGTTTTTGCAGACAATACTGGCATGAAAGGGCACGTGCGGCTTTTGCGACAgggtgagt                                                                                                                                      >  2:680086/1‑245 (MQ=255)
                                                                                                      gACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATGATGCCTTTCACGTCACCATAGAGGTCGGTGTGGCTGTGAGGGGAGATACCGGCAACAATCGGTTTTTGCGGACGATGCTGGACGTAAAGGTCAAGTGCGGCTTTTGCGGccggggggggc                                                                                                                                     >  2:167992/1‑241 (MQ=255)
                                                                                                        ccgaaagaggggggcgccgggcggcTCAGTGTGTGACCCGCCAGCACGTTTTCGCTGATGGCTTCGGCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGGGGCGTCAc                                                                                                                                   <  2:341804/232‑1 (MQ=255)
                                                                                                           tAAGAGTAGAGCGCCCGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGCGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTAAAGTGCGGCTTTTGCGGCAGGCGGCCTCAcc                                                                                                                                  >  1:517344/1‑249 (MQ=255)
                                                                                                                 tAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAt                                                                                                                                                                                                              >  2:419251/1‑167 (MQ=255)
                                                                                                                 aagagCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAt                                                                                                                                                                                                              <  1:419251/166‑1 (MQ=255)
                                                                                                                       gcACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAAc                                                                                                                   >  1:66573/1‑251 (MQ=255)
                                                                                                                           ggcggcTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCCGGTCGATAACaat                                                                                                                >  1:253666/1‑250 (MQ=255)
                                                                                                                              ggcTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                                                                   <  1:500898/197‑1 (MQ=255)
                                                                                                                              ggcTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                                                                   >  2:500898/1‑197 (MQ=255)
                                                                                                                                   atgatgTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCttt                                                                                                       <  1:20304/251‑1 (MQ=255)
                                                                                                                                      gtttgtACCCGCCCGCGCGTTGTCGCGGATGGCTCCTTCCATAACGCTGGCAGTGGCAATGCCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGTTACTGCCTTTCACGTCACCATAGGTGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTtc                                                                                                      <  2:321536/245‑1 (MQ=255)
                                                                                                                                          ttttccccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATAACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTtctc                                                                                                    <  2:500333/245‑1 (MQ=255)
                                                                                                                                            taccccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGc                                                                                                  <  2:162613/247‑1 (MQ=255)
                                                                                                                                                  gCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGGGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCGCTAACAATAATGCCCTTTTCGCCCTCaac                                                                                           >  1:601345/1‑251 (MQ=255)
                                                                                                                                                      gCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaacgaac                                                                                       <  1:195630/251‑1 (MQ=255)
                                                                                                                                                          gTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCCACGAACGTa                                                                                    >  1:618395/1‑250 (MQ=255)
                                                                                                                                                               cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGAGCGGCTTTTGCGGCAGGGGGCGTCACCCGCGGGTCGGTAACAATAATACCTTTCTAGCCCCCAACGAACGTAAAGttt                                                                              >  1:315904/1‑250 (MQ=255)
                                                                                                                                                               cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAGAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTa                                                                              <  1:514004/251‑1 (MQ=255)
                                                                                                                                                                gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGACCTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAgg                                                                                                                                                                                                               >  2:288684/1‑121 (MQ=255)
                                                                                                                                                                gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCGGTAAAAATAATGCCTTTCTCGCCCTCCACGAACGTAATGTTa                                                                              >  1:663260/1‑250 (MQ=255)
                                                                                                                                                                gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGgaga                                                                                                                                                                                                               <  1:288684/122‑1 (MQ=255)
                                                                                                                                                                    aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGat                                                                         <  1:283209/251‑1 (MQ=255)
                                                                                                                                                                     tGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAAt                                                                                                                                                                                                  >  1:420652/1‑130 (MQ=255)
                                                                                                                                                                     tGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCGGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAAt                                                                                                                                                                                                  <  2:420652/130‑1 (MQ=255)
                                                                                                                                                                      ggCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACGTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGGCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTGTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGatat                                                                       <  1:567699/251‑1 (MQ=255)
                                                                                                                                                                         ttCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGatatact                                                                    >  1:598803/1‑248 (MQ=255)
                                                                                                                                                                              ccttAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAGCTGCTTCAGAGATAATGCCTTTCACGTCCCCATAGTGGACGGTGTGGCTGGGAGTGTAGATAACGGAAACAACCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCGTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAAGGCCTTTCTCGCCCCCAACGAACGTAATGTTAGAGATATCCTGGc                                                                 <  2:455716/246‑1 (MQ=255)
                                                                                                                                                                                aTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTg                                                                                                                                                                                         >  1:177506/1‑128 (MQ=255)
                                                                                                                                                                                aTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTg                                                                                                                                                                                         <  2:177506/128‑1 (MQ=255)
                                                                                                                                                                                 tAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAACGGGTCGGATACCATAATGCCTTTCTCGCCCTCAACGGACTTACTTTTAGAGATGTCCTGGCCgcgc                                                            >  1:651334/1‑251 (MQ=255)
                                                                                                                                                                                         gCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCGGGCCGCGCACCTGa                                                      <  2:661379/249‑1 (MQ=255)
                                                                                                                                                                                         gCAGGGGCAATAACCTGAACTTTGCCGGATTTAACAGCGGCGTCAGAGATAATGCCGTTCACGTCACCAGAGTGGGCGGTGGGGCTGTGAGTGTAGATAACGGCAACAATCGGTTGTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGa                                                      <  2:352896/249‑1 (MQ=255)
                                                                                                                                                                                             ggtcAATCATCGGAACTTTGCCGGATTTAACAGCGGCTTTAGAGAGAATGCCTTTCACGTCACCATAGTGGGCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGGCGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCGATAACAATAATGCCTTTCTCTCCCGCAACGAACGTAATGTTAGAGATATCCTGCCCGCGCACCGGATACAtt                                                <  2:51847/248‑1 (MQ=255)
                                                                                                                                                                                                  aTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATc                                           <  1:185868/251‑1 (MQ=255)
                                                                                                                                                                                                     aCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGt                                        <  1:173452/251‑1 (MQ=255)
                                                                                                                                                                                                         gAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACt                                    <  2:663260/251‑1 (MQ=255)
                                                                                                                                                                                                                  cGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCCACGACCGTGATGTTAGAGATATCCTTGCCGCGCACCTGATACATTTTTTCGGTGACTTTTAACAgg                           >  1:185812/1‑251 (MQ=255)
                                                                                                                                                                                                                           cATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGCTAACAATAATGCCTTTCTCGCCCTCAACGCACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAt                                                 >  1:473402/1‑220 (MQ=255)
                                                                                                                                                                                                                           cATCGGCTCCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGGAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAt                                                 <  2:473402/220‑1 (MQ=255)
                                                                                                                                                                                                                                 cTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGCAATACCGTTg            >  1:30448/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGAt          <  1:268833/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACt  >  1:415481/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCGGGGCGGTAACAATAATGCCTTTCTCGCCCCCAACGAACGGAATGTTAGAGATATCCTGGCCgcgc                                                            >  1:256544/1‑188 (MQ=255)
                                                                                                                                                                                                                                                cTTTCACGGCACCATAGTGGTCGGTGTGGCTGGGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCgcgc                                                            <  2:256544/188‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
GACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACT  >  NC_000913/4305791‑4306272

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