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

GACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACT  >  NC_000913/4305791‑4306272
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ggcgatatgtttggggggTGCAATAAGGTTGGGTTCGCCCGTTGCCAGCGTCAC‑TCCAAGACCATTCCCCCCTTTGCCTGGCGGTTTGGGGGGCGGAAACAGACCGTAAGAGTGGCGCGCAGTGCGCCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGtcac                                                                                                                                                                                                                                             <  2:658137/232‑1 (MQ=255)
   gATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAgg                                                                                                                                                                                                                                       >  1:661379/1‑250 (MQ=255)
       gTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTc                                                                                                                                                                                                                                   >  1:500333/1‑251 (MQ=255)
        ttttCGACGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGCTAATGCCTTTCACGTCACCATAGTTGTCg                                                                                                                                                                                                                                  >  2:283209/1‑251 (MQ=255)
            cgtcgATGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCTCATTGCCTTGCACGTTGCGCTGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACATTTTCGCTGATGGCTTCGGACCTCAAGCCTGCAGGGGAAATCCCCTGAACTTTGCCGGATTTACCCTCTGCTTCACGGTGAATGCCTTTCACGTCACCATaagggacggggt                                                                                                                                                                                                                              >  2:514004/1‑240 (MQ=255)
               cgGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGc                                                                                                                                                                                                                           >  1:162613/1‑251 (MQ=255)
                ggTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑TCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCt                                                                                                                                                                                                                          <  1:213522/251‑1 (MQ=255)
                   gaaataatGCTCGGGCCGCCCTTTGCCAGCGCCAC‑GCCAAGGCCTTTCCCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGTTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCtgtg                                                                                                                                                                                                                       <  1:233732/249‑1 (MQ=255)
                            cTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTgtgtgtgt                                                                                                                                                                                                                  >  1:381815/1‑247 (MQ=255)
                                ggTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACg                                                                                                                                                                                                          >  1:678145/1‑251 (MQ=255)
                                ggTCGCCCGTTGCCAGCGTCAA‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGCTGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATACAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTACACGTCACCATAGTGGTCGGTGTGGCTGTGAgttgagtttgcg                                                                                                                                                                                                          >  2:295604/1‑241 (MQ=255)
                                    gCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGACGGTGTGGCtggggggggagttaacggca                                                                                                                                                                                                       >  2:55213/1‑232 (MQ=255)
                                     cccGTTGCCAG‑GTCCCCGCCAGGGCCATTACCCACATTGCCTGGCGCGTTGGGCGGCAGAAACAGACCGTAAGAGTAGAGCGCACGGGGGCTCT‑TGATGTT‑C‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGcaac                                                                                                                                                                                                     <  1:463288/251‑1 (MQ=255)
                                         ttGCCAGAGCCAC‑GTCAAGGCCATTACCCACATTGCCTTGCGCATTGTGCGGCAGTAACCGACCGTAAGAGCAGATCGCACGGCGCCTCA‑TGATGTT‑C‑CCCGACAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGACTGCAGGGGCAATCACCTGAACTTTGCAGGATTTAACATCGGCCTCACAGATAATGCCTTTACCGTCACCATCGTGGTCGGTGTGGCTGTGAGAGTAGTTAACGGcacccca                                                                                                                                                                                                  >  2:81237/1‑245 (MQ=255)
                                             cAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTACCAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGGGGGAGTAGGTAACCGCAAAAAACGGtt                                                                                                                                                                                             >  2:332364/1‑251 (MQ=255)
                                             cAGCGACAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCGGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGAGAAGGCCTTTCACGTAACCATAGTGGTCGGTGTGGCTGGGGATGTAAATAACGGCAAAAATCGGtt                                                                                                                                                                                             >  2:152768/1‑251 (MQ=255)
                                               gCGTCAC‑TCCAAGGCCTTTACCCACATTGCCGTGCGCGTTGGGCGGCAGTAACAGACCGTAAGGGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCAGAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGtttt                                                                                                                                                                                           <  1:567662/251‑1 (MQ=255)
                                                   cAC‑GCCAATGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTCACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGTGGCAATCACCTGCACTTTGCCGGATTTAACATCGGCTTCTGGGGTCATTCCTTTCACGTCACCATAGTGGTCGGTGAGGCTGGGTGTTTAGATAACGGCAACCATCCGTTTTTg                                                                                                                                                                                         >  2:20304/1‑249 (MQ=255)
                                                          aaGGCCATTACACACATTGCCTTGCGCGTTGTGCGGCAGTAACAGGACGAAAAAGTAGAGCGCAAGGCGGCTCA‑TGATGTT‑A‑CACGCAAGCACGTTTTCGATGATTGCTGCGTCCATAAAGCCTGCGGGGTCAATCACCTGTAATTTGCCGGATTTAACATCTGCTTCGGAGATAATGCCTTTCCCGTCGACATAGGGGTCGGTGTGGCTGAGAGTAGAGAAAGCGGCATCAATCGGgtttgtcaggcgat                                                                                                                                                                                 >  2:615995/1‑241 (MQ=255)
                                                           aGGCCATTACCCACATTGCCTTGCGCGTTGGGCGGCAGTAACAGACCGTAAGGGGAGAGCGCCCGGCGTCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATg                                                                                                                                                                                <  1:360453/251‑1 (MQ=255)
                                                           aGGCCATTACCAACGTGGCCTGGCGCGTGGTGCGGCAGTAACAGACCGTAAGAGTATATCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCGGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATg                                                                                                                                                                                <  1:615995/251‑1 (MQ=255)
                                                              ccATTACCCACATTGCCTTGCGCGTTGGGCGGAAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTg                                                                                                                                                                             <  1:168100/251‑1 (MQ=255)
                                                                 ttACCCACATTGCCTTGCGCGTTGAGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCCCCTGAAATTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCAAGTCACCATAGTGGTCGGCGTGGCTGTGAGTGTAGATAAAGGCACCAATCGGTTTTTGCGGACGATGCTGGa                                                                                                                                                                           >  2:323196/1‑250 (MQ=255)
                                                                   aCCCACATTGCCGTGCGGGTTGTGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGTGGGCTCA‑TGATGTT‑C‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAg                                                                                                                                                                         <  1:580231/250‑1 (MQ=255)
                                                                       ataTTGCCTTGCGCGTTGGGCGGCTTTAACAGCCCAAAAGAATAGAGCGCTCGGCGGTCCA‑TG‑GGTT‑ATCCCGCCAGCACGTTTTCGCTGATGGCTTCTTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGGCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAg                                                                                                                                                                    <  2:549239/249‑1 (MQ=255)
                                                                        caTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                                                                   >  1:239761/1‑251 (MQ=255)
                                                                             ccTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCa                                                                                                                                                                >  1:113493/1‑249 (MQ=255)
                                                                                             gCAGTAACAGACCGTAAGCGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATCAAGCCTGCAGGGGCAATCCCCTGCACTTTGCCGGTTTTAACATCGGCTTCAGAGATAATGCCCTTCACGTCACCATAGGGGTCGGTGTGGCTGTGAGTGTCGATAACGGCACCAATCGGTTTTTGCGGACGATGCTGGACGTACAGGTCAAGAGCGGCTTTTGCgtc                                                                                                                                               >  1:431047/1‑248 (MQ=255)
                                                                                               atgatcagacagaaggtggggagCGCCCGGTGGCTCG‑TG‑TGTTCA‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGACTTAACATCGGCTTCAGAGATAATGCCTTTCGCGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAgg                                                                                                                                            <  2:195812/231‑1 (MQ=255)
                                                                                                     aGACCGTCAGAGTAGAGCGCCCGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGATTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGTTAATGCCTTTAACGTCACCATCGTGGTCGGAGTGGCGGTGAGTGTAAGTAACGGGAACCATCGGTTTTTGCAGACAATACTGGCATGAAAGGGCACGTGCGGCTTTTGCGACAgggtgagt                                                                                                                                      >  2:680086/1‑245 (MQ=255)
                                                                                                      gACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATGATGCCTTTCACGTCACCATAGAGGTCGGTGTGGCTGTGAGGGGAGATACCGGCAACAATCGGTTTTTGCGGACGATGCTGGACGTAAAGGTCAAGTGCGGCTTTTGCGGccggggggggc                                                                                                                                     >  2:167992/1‑241 (MQ=255)
                                                                                                        ccgaaagaggggggcgccgggcggcTCAGTG‑TGTG‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGGCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGGGGCGTCAc                                                                                                                                   <  2:341804/232‑1 (MQ=255)
                                                                                                           tAAGAGTAGAGCGCCCGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGCGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTAAAGTGCGGCTTTTGCGGCAGGCGGCCTCAcc                                                                                                                                  >  1:517344/1‑249 (MQ=255)
                                                                                                                 tAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAt                                                                                                                                                                                                              >  2:419251/1‑167 (MQ=255)
                                                                                                                 aagagCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAt                                                                                                                                                                                                              <  1:419251/166‑1 (MQ=255)
                                                                                                                       gcACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAAc                                                                                                                   >  1:66573/1‑251 (MQ=255)
                                                                                                                           ggcggcTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCCGGTCG‑ATAACaat                                                                                                                >  1:253666/1‑250 (MQ=255)
                                                                                                                              ggcTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                                                                   <  1:500898/197‑1 (MQ=255)
                                                                                                                              ggcTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAgg                                                                                                                                                                   >  2:500898/1‑197 (MQ=255)
                                                                                                                                   a‑tgatgTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCttt                                                                                                       <  1:20304/251‑1 (MQ=255)
                                                                                                                                      gtttgt‑A‑CCCGCCCGCGCGTTGTCGCGGATGGCTCCTTCCATAACGCTGGCAGTGGCAATGCCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGTTACTGCCTTTCACGTCACCATAGGTGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTtc                                                                                                      <  2:321536/245‑1 (MQ=255)
                                                                                                                                          ttttccccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATAACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTtctc                                                                                                    <  2:500333/245‑1 (MQ=255)
                                                                                                                                            taccccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGc                                                                                                  <  2:162613/247‑1 (MQ=255)
                                                                                                                                                  gCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGGGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCG‑CTAACAATAATGCCCTTTTCGCCCTCaac                                                                                           >  1:601345/1‑251 (MQ=255)
                                                                                                                                                      gCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCaacgaac                                                                                       <  1:195630/251‑1 (MQ=255)
                                                                                                                                                          gTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCCACGAACGTa                                                                                    >  1:618395/1‑250 (MQ=255)
                                                                                                                                                               cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGAGCGGCTTTTGCGGCAGGGGGCGTCACCCGCGGGTCG‑GTAACAATAATACCTTTCTAGCCCCCAACGAACGTAAAGttt                                                                              >  1:315904/1‑250 (MQ=255)
                                                                                                                                                               cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAGAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTa                                                                              <  1:514004/251‑1 (MQ=255)
                                                                                                                                                                gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGACCTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAgg                                                                                                                                                                                                               >  2:288684/1‑121 (MQ=255)
                                                                                                                                                                gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCG‑GTAAAAATAATGCCTTTCTCGCCCTCCACGAACGTAATGTTa                                                                              >  1:663260/1‑250 (MQ=255)
                                                                                                                                                                gCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGgaga                                                                                                                                                                                                               <  1:288684/122‑1 (MQ=255)
                                                                                                                                                                    aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGat                                                                         <  1:283209/251‑1 (MQ=255)
                                                                                                                                                                     tGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAAt                                                                                                                                                                                                  >  1:420652/1‑130 (MQ=255)
                                                                                                                                                                     tGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCGGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAAt                                                                                                                                                                                                  <  2:420652/130‑1 (MQ=255)
                                                                                                                                                                      ggCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACGTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGGCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTGTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGatat                                                                       <  1:567699/251‑1 (MQ=255)
                                                                                                                                                                         ttCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGatatact                                                                    >  1:598803/1‑248 (MQ=255)
                                                                                                                                                                              ccttAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAGCTGCTTCAGAGATAATGCCTTTCACGTCCCCATAGTGGACGGTGTGGCTGGGAGTGTAGATAACGGAAACAACCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCGTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAAGGCCTTTCTCGCCCCCAACGAACGTAATGTTAGAGATATCCTGGc                                                                 <  2:455716/246‑1 (MQ=255)
                                                                                                                                                                                aTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTg                                                                                                                                                                                         >  1:177506/1‑128 (MQ=255)
                                                                                                                                                                                aTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTg                                                                                                                                                                                         <  2:177506/128‑1 (MQ=255)
                                                                                                                                                                                 tAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAACGGGTCGGATA‑CCATAATGCCTTTCTCGCCCTCAACGGACTTACTTTTAGAGATGTCCTGGCCgcgc                                                            >  1:651334/1‑251 (MQ=255)
                                                                                                                                                                                         gCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCGGGCCGCGCACCTGa                                                      <  2:661379/249‑1 (MQ=255)
                                                                                                                                                                                         gCAGGGGCAATAACCTGAACTTTGCCGGATTTAACAGCGGCGTCAGAGATAATGCCGTTCACGTCACCAGAGTGGGCGGTGGGGCTGTGAGTGTAGATAACGGCAACAATCGGTTGTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGa                                                      <  2:352896/249‑1 (MQ=255)
                                                                                                                                                                                             ggtcAATCATCGGAACTTTGCCGGATTTAACAGCGGCTTTAGAGAGAATGCCTTTCACGTCACCATAGTGGGCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGGCGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCG‑ATAACAATAATGCCTTTCTCTCCCGCAACGAACGTAATGTTAGAGATATCCTGCCCGCGCACCGGATACAtt                                                <  2:51847/248‑1 (MQ=255)
                                                                                                                                                                                                  aTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATc                                           <  1:185868/251‑1 (MQ=255)
                                                                                                                                                                                                     aCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGt                                        <  1:173452/251‑1 (MQ=255)
                                                                                                                                                                                                         gAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACt                                    <  2:663260/251‑1 (MQ=255)
                                                                                                                                                                                                                  cGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCCACGACCGTGATGTTAGAGATATCCTTGCCGCGCACCTGATACATTTTTTCGGTGACTTTTAACAgg                           >  1:185812/1‑251 (MQ=255)
                                                                                                                                                                                                                           cATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑CTAACAATAATGCCTTTCTCGCCCTCAACGCACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAt                                                 >  1:473402/1‑220 (MQ=255)
                                                                                                                                                                                                                           cATCGGCTCCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGGAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAt                                                 <  2:473402/220‑1 (MQ=255)
                                                                                                                                                                                                                                 cTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGCAATACCGTTg            >  1:30448/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGAt          <  1:268833/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACt  >  1:415481/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCGGGGCG‑GTAACAATAATGCCTTTCTCGCCCCCAACGAACGGAATGTTAGAGATATCCTGGCCgcgc                                                            >  1:256544/1‑188 (MQ=255)
                                                                                                                                                                                                                                                cTTTCACGGCACCATAGTGGTCGGTGTGGCTGGGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCgcgc                                                            <  2:256544/188‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
GACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCAC‑GCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCA‑TGATGTT‑A‑CCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCG‑ATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACT  >  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: ‑