Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 619,395:1 +A 5.7% coding (606/1329 nt) pfeS‑I ← histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947619,3951.A 5.7% 276.7 / 10.4 70coding (606/1329 nt)pfeS‑Ihistidine kinase
Reads supporting (aligned to +/- strand):  ref base . (33/33);  new base A (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GCCGCGCCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATGGGTCCAGG  >  NC_002947/619250‑619543
                                                                                                                                                  |                                                                                                                                                    
gCCGCGCCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAGGGCTCTGGAGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCgg                                                                                                                                                  >  1:629478/1‑150 (MQ=255)
 ccGCGCCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCggg                                                                                                                                                 <  2:851839/150‑1 (MQ=255)
 ccGCGCCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCggg                                                                                                                                                 <  1:595063/150‑1 (MQ=255)
   gcgcCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCgggg                                                                                                                                                <  1:715457/149‑1 (MQ=255)
   gcgcCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCgggg                                                                                                                                                >  2:715457/1‑149 (MQ=255)
   gcgcCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGc                                                                                                                                               >  1:14678/1‑150 (MQ=255)
     gcCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTg                                                                                                                                             <  2:593332/150‑1 (MQ=255)
      ccAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGt                                                                                                                                            >  2:1113134/1‑150 (MQ=255)
       cAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTc                                                                                                                                           >  2:246428/1‑150 (MQ=255)
       cAGCGGCCTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAGGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAGGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGTAGGGGCAGG‑GCGGGGCTGTc                                                                                                                                           <  2:23899/150‑1 (MQ=255)
        aGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTcc                                                                                                                                          <  2:790103/150‑1 (MQ=255)
         gcggcgTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCa                                                                                                                                         <  1:1113134/150‑1 (MQ=255)
           ggcgTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGa                                                                                                                                       <  2:108525/150‑1 (MQ=255)
           ggcgTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGa                                                                                                                                       <  2:1123116/150‑1 (MQ=255)
               tGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGatcat                                                                                                                                   >  2:149116/1‑150 (MQ=255)
                  gcAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGataatctg                                                                                                                                >  2:479439/1‑148 (MQ=255)
                         gTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAg                                                                                                                         <  1:468425/150‑1 (MQ=255)
                         gTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAg                                                                                                                         >  2:814418/1‑150 (MQ=255)
                          tGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAgg                                                                                                                        <  1:282306/150‑1 (MQ=255)
                            agaacccGTGCGTAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGc                                                                                                                      >  2:599585/7‑150 (MQ=255)
                            agaCAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGc                                                                                                                      <  1:1093820/150‑1 (MQ=255)
                                     gcgcAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGcgcg                                                                                                             <  1:246428/150‑1 (MQ=255)
                                       gcAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCAACCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGACGCGGt                                                                                                           >  1:682676/1‑150 (MQ=255)
                                         aaCAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTcg                                                                                                         >  2:1122620/1‑150 (MQ=255)
                                                      gTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCa                                                                                            <  1:534816/150‑1 (MQ=255)
                                                       ttGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCNAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAg                                                                                           <  1:479439/150‑1 (MQ=255)
                                                       ttGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAg                                                                                           <  2:308776/150‑1 (MQ=255)
                                                        tGCGCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGc                                                                                          <  1:718441/150‑1 (MQ=255)
                                                          cTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCgg                                                                                        >  2:723677/1‑150 (MQ=255)
                                                           tCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCgg                                                                                                                                                  <  1:410831/91‑1 (MQ=255)
                                                           tCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCgg                                                                                                                                                  >  2:410831/1‑91 (MQ=255)
                                                           tCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGc                                                                                       >  1:176701/1‑150 (MQ=255)
                                                           tCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGc                                                                                       >  1:996608/1‑150 (MQ=255)
                                                             aaGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCggg                                                                                                                                                 >  1:1073373/1‑90 (MQ=255)
                                                             aaGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCggg                                                                                                                                                 <  2:1073373/90‑1 (MQ=255)
                                                                       cgcAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGc                                                                           <  1:1057339/150‑1 (MQ=255)
                                                                                ccGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAgc                                                                  >  1:1207491/1‑150 (MQ=255)
                                                                                ccGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAgc                                                                  <  2:14678/150‑1 (MQ=255)
                                                                                  ggCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCCGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAgcgc                                                                <  1:439795/150‑1 (MQ=255)
                                                                                      aTGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAgg                                                            >  1:1200476/1‑150 (MQ=255)
                                                                                              gAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTagcag                                                    >  1:468944/1‑150 (MQ=255)
                                                                                               aaGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTagcagc                                                   <  1:599585/150‑1 (MQ=255)
                                                                                                   ccTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGc                                                                                                  >  2:553215/1‑99 (MQ=255)
                                                                                                   ccTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGc                                                                                                  <  1:553215/99‑1 (MQ=255)
                                                                                                    cTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACaa                                                                                >  1:973868/1‑116 (MQ=255)
                                                                                                    cTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACaa                                                                                <  2:973868/116‑1 (MQ=255)
                                                                                                       aGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCcag                                           <  2:208942/150‑1 (MQ=255)
                                                                                                         ccAGCTAACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCcagca                                         <  2:248075/150‑1 (MQ=255)
                                                                                                            gCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCGAGATCATAGGCCCGCAGTGCGTCGGCGCGGTCGCGCAGGCGGTTGAGCGGCACGGCAAGATGTCGGTGTAGCGCAAGGCCTAGCAGCACGGCcagcgtta                                      >  2:935571/1‑146 (MQ=255)
                                                                                                               cACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCgg                                   <  1:814418/150‑1 (MQ=255)
                                                                                                               cACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCgg                                   <  2:993262/150‑1 (MQ=255)
                                                                                                                aCCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCcagca                                         >  2:358301/1‑143 (MQ=255)
                                                                                                                aCCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCcagca                                         <  1:358301/143‑1 (MQ=255)
                                                                                                                  gatctggacagccCCGC‑CCTGCCCCTGCAGGAGCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAgcagca                                                  <  1:180222‑M1/119‑1 (MQ=255)
                                                                                                                  gatctggacagccCCGC‑CCTGCCCCTGCAGGAGCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAgcagca                                                  >  2:180222‑M1/14‑132 (MQ=255)
                                                                                                                   cAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGc                               >  2:300124/1‑150 (MQ=255)
                                                                                                                      tggacagccCCGC‑CCTGCCCCTGCAGGAGCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGcagcag                                                 <  1:761380‑M1/120‑1 (MQ=255)
                                                                                                                      tggacagccCCGC‑CCTGCCCCTGCAGGAGCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGcagcag                                                 >  2:761380‑M1/10‑129 (MQ=255)
                                                                                                                      cTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCcagca                                         <  1:469202/137‑1 (MQ=255)
                                                                                                                      cTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCcagca                                         >  2:469202/1‑137 (MQ=255)
                                                                                                                      cTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAAc                            >  2:1150102/1‑150 (MQ=255)
                                                                                                                          tcaCGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCa                        >  1:39018/1‑150 (MQ=255)
                                                                                                                          tcaCGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCa                        <  1:300124/150‑1 (MQ=255)
                                                                                                                           caCGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCAt                       >  2:985292/1‑150 (MQ=255)
                                                                                                                              gccgcTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCg                    >  1:843635/1‑150 (MQ=255)
                                                                                                                              gccgcTCCTGCAGGGGCAGG‑GCGGGGCCGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAgc                                                                  >  1:823489/1‑104 (MQ=255)
                                                                                                                              gccccTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAgc                                                                  <  2:823489/104‑1 (MQ=255)
                                                                                                                                 GCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCAGCAGCGCGTCGGCGCGGTCGCGCAGGCGGTTCATCCGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGa                 >  2:739458‑M1/1‑150 (MQ=255)
                                                                                                                                  CTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGAc                >  1:135918‑M1/1‑150 (MQ=255)
                                                                                                                                   tCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGAcc               >  1:20614/1‑150 (MQ=255)
                                                                                                                                       gcaggggcagg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACAt           <  2:996608/150‑1 (MQ=255)
                                                                                                                                         aggggcagg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGt                                                                                                           >  1:1207970/1‑52 (MQ=255)
                                                                                                                                         aggggcagg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGt                                                                                                           <  2:1207970/52‑1 (MQ=255)
                                                                                                                                         aggggcagg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATgg         <  2:475978/150‑1 (MQ=255)
                                                                                                                                         aggggcagg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGACGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATgg         <  2:682676/150‑1 (MQ=255)
                                                                                                                                          gtttcagg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATggg        >  1:475978/5‑150 (MQ=255)
                                                                                                                                          ggggcagg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATggg        >  1:636200/1‑150 (MQ=255)
                                                                                                                                               agg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCTGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATGGGTCCAg   <  2:15893/150‑1 (MQ=255)
                                                                                                                                                gg‑gCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATGGGTCCAgg  <  1:985292/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCGCGCCAGCGGCGTGCGCAATTCGTGAGACAGCGTGCGCAACAGCAGGCGCTGTTGCTCAAGGCTCTGGCGCAGGCGCCCGGCCATGTGCTCGAAGGCCTGAGCCAGCTCACCCAGCTCATCACGCCGCTCCTGCAGGGGCAGG‑GCGGGGCTGTCCAGATCATCGGCCCGCAGGGCGTCGGCGCGGTCGCGCAGGCGGTTCAGCGGCACGACAAGATGCCGGTATAGCGCAAGGCCTAGCAGCAGGGCCAGCAGTGCGGGGGCAACGCCATGCGTGACCACATGGGTCCAGG  >  NC_002947/619250‑619543

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