Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 739,598 G→C 5.2% G344G (GGC→GGG) ‡ pilC ← type IV fimbrial assembly protein PilC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947739,5980GC 5.2% 707.0 / 20.6 210G344G (GGC→GGG) ‡pilCtype IV fimbrial assembly protein PilC
Reads supporting (aligned to +/- strand):  ref base G (109/90);  new base C (7/4);  total (116/94)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 210 total aligned reads displayed.

GCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCACCCGTTAC  >  NC_002947/739450‑739744
                                                                                                                                                    |                                                                                                                                                  
gCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGc                                                                                                                                                   <  1:4929803/150‑1 (MQ=255)
gCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGc                                                                                                                                                   >  2:3821331/1‑150 (MQ=255)
    gTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCt                                                                                                                                               >  1:2657275/1‑150 (MQ=255)
    gTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCt                                                                                                                                               >  2:1500677/1‑150 (MQ=255)
       cATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGa                                                                                                                                            <  2:4392864/150‑1 (MQ=255)
        aTCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGGCACGCTGGac                                                                                                                                           <  1:3449281/150‑2 (MQ=255)
         tCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGAtt                                                                                                                                          >  2:1384744/1‑150 (MQ=255)
           gCCACCACCAGGCCACACACCAGCAAGCCCAGGATCAGAAAGATGGACGGCTCCATCAGGCTGGTCAACTGATACAGAGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCTGCTATTTAAAGCATCTGGTCCAGAGTGCCGCTGGATTCa                                                                                                                                        >  2:2508148/1‑150 (MQ=255)
            ccaccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGc                                                                                                                                                <  1:2487852/141‑1 (MQ=255)
            ccaccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCAc                                                                                                                                       >  2:3051545/1‑150 (MQ=255)
            ccaccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCAc                                                                                                                                       <  2:4656864/150‑1 (MQ=255)
            ccaccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCAc                                                                                                                                       >  2:3051922/1‑150 (MQ=255)
            ccaccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGGCACGCTGGaccaga                                                                                                                                       <  1:4659090/150‑6 (MQ=255)
                caccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGAt                                                                                                                                   <  1:4434792/150‑1 (MQ=255)
                   caGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGc                                                                                                                                <  1:2757472/150‑1 (MQ=255)
                         acccaccAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCaa                                                                                                                          >  1:1331460/1‑150 (MQ=255)
                           ccaccaGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACt                                                                                                                        <  2:2592365/150‑1 (MQ=255)
                           ccaccaGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACt                                                                                                                        >  1:1767372/1‑150 (MQ=255)
                             accaGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGc                                                                                                                      >  2:3092602/1‑150 (MQ=255)
                                 gCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCa                                                                                                                  <  2:4741909/150‑1 (MQ=255)
                                   aaGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGc                                                                                                                >  1:1127521/1‑150 (MQ=255)
                                   aaGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGc                                                                                                                >  1:775086/1‑150 (MQ=255)
                                      cccAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGACCCAGCAAg                                                                                                             >  1:1662764/1‑150 (MQ=255)
                                       ccAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgg                                                                                                            >  2:2304970/1‑150 (MQ=255)
                                         aGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgggg                                                                                                          <  2:5396238/150‑1 (MQ=255)
                                           gATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgggggg                                                                                                        <  2:676478/150‑1 (MQ=255)
                                            aTCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAggggggg                                                                                                       <  2:5414284/150‑1 (MQ=255)
                                              cAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGaa                                                                                                     <  2:4237764/150‑1 (MQ=255)
                                              cAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGaa                                                                                                     <  2:2226327/150‑1 (MQ=255)
                                                gCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCACGGGGGGGAACa                                                                                                   >  2:145310/1‑150 (MQ=255)
                                                gCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACa                                                                                                   >  2:1617973/1‑150 (MQ=255)
                                                gCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACa                                                                                                   >  2:1210932/1‑150 (MQ=255)
                                                 cACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAg                                                                                                  >  2:2664627/1‑150 (MQ=255)
                                                 cACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAg                                                                                                  <  2:4980557/150‑1 (MQ=255)
                                                   cGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCg                                                                                                <  1:2901823/150‑1 (MQ=255)
                                                         ccGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCAc                                                                                          <  2:2468308/150‑1 (MQ=255)
                                                         ccGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGGCACGCTGGac                                                                                                                                           <  1:2997837/101‑2 (MQ=255)
                                                           ggCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCg                                                                                        <  1:201622/150‑1 (MQ=255)
                                                           ggCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCg                                                                                        <  1:3518824/150‑1 (MQ=255)
                                                            gCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCgg                                                                                       >  2:2095845/1‑150 (MQ=255)
                                                                 agcagGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATc                                                                                  >  1:2575562/1‑150 (MQ=255)
                                                                         ggTCAACTGATCCAGCGCCTGGCTGCCCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGAGTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGtt                                                                          >  2:3261873/1‑150 (MQ=255)
                                                                           tCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCa                                                                        >  1:1052705/1‑150 (MQ=255)
                                                                           tCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCa                                                                        <  1:37303/150‑1 (MQ=255)
                                                                              aCTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGGAAGGGGGGGACCAGCGGTTCGCCGGCCATCGCCGGGGTCAGcc                                                                     >  2:2921150/1‑150 (MQ=255)
                                                                               cTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGccc                                                                    <  1:2920373/150‑1 (MQ=255)
                                                                                 gATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCtt                                                                  >  2:2532624/1‑150 (MQ=255)
                                                                                  aTCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCgg                                                                                               >  1:4461895/1‑120 (MQ=255)
                                                                                  aTCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGGCACGCTGGac                                                                                                                                           >  1:1371904/1‑75 (MQ=255)
                                                                                   tCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGACCAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGc                                                                >  1:4116409/1‑150 (MQ=255)
                                                                                   tCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGc                                                                <  1:3858925/150‑1 (MQ=255)
                                                                                   tCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGc                                                                >  1:4116028/1‑150 (MQ=255)
                                                                                     cAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGccc                                                              >  1:417841/1‑150 (MQ=255)
                                                                                       gcgcCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCAt                                                            <  1:2378846/150‑1 (MQ=255)
                                                                                         gcCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTg                                                          <  2:3368329/150‑1 (MQ=255)
                                                                                                tGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGc                                                   >  1:2390594/1‑150 (MQ=255)
                                                                                                tGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGc                                                   <  2:5345636/150‑1 (MQ=255)
                                                                                                    cTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTg                                               <  1:2200974/150‑1 (MQ=255)
                                                                                                    cTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTg                                               >  2:451137/1‑150 (MQ=255)
                                                                                                    cTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTg                                               >  1:5180430/1‑150 (MQ=255)
                                                                                                     tGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTgg                                              >  2:1383821/1‑150 (MQ=255)
                                                                                                     tGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTgg                                              <  1:4321916/150‑1 (MQ=255)
                                                                                                     tGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTgg                                              <  1:3584030/150‑1 (MQ=255)
                                                                                                     tGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTgg                                              >  2:107181/1‑150 (MQ=255)
                                                                                                      ggtCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgc                                             <  1:2497737/148‑1 (MQ=255)
                                                                                                      gTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCGACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgc                                             <  1:125194/150‑1 (MQ=255)
                                                                                                      gTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgc                                             >  2:4314550/1‑150 (MQ=255)
                                                                                                      gTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgc                                             <  2:2374008/150‑1 (MQ=255)
                                                                                                        tCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgcgc                                           <  2:5029782/150‑1 (MQ=255)
                                                                                                         ccTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCTCATTGGCCATGCCTTGGCGCa                                          >  1:2066998/1‑150 (MQ=255)
                                                                                                          cTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAg                                         >  1:3698200/1‑150 (MQ=255)
                                                                                                           tCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGc                                        >  2:1344619/1‑150 (MQ=255)
                                                                                                           tCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGc                                        >  2:3171865/1‑150 (MQ=255)
                                                                                                              tAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGc                                     >  2:1391144/1‑150 (MQ=255)
                                                                                                                aTGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCaa                                   <  2:1030467/150‑1 (MQ=255)
                                                                                                                     tggcGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTc                                                                         >  2:2404607/1‑107 (MQ=255)
                                                                                                                     tggcGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCg                              >  2:4164594/1‑150 (MQ=255)
                                                                                                                      ggcggcTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGc                             >  2:2960542/1‑150 (MQ=255)
                                                                                                                        cggcTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCt                           <  1:3286220/150‑1 (MQ=255)
                                                                                                                         ggcTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTg                          >  1:4900846/1‑150 (MQ=255)
                                                                                                                         ggcTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTg                          >  1:5059168/1‑150 (MQ=255)
                                                                                                                         ggcTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTg                          >  1:105963/1‑150 (MQ=255)
                                                                                                                            tttttCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTc                       >  2:3690868/1‑150 (MQ=255)
                                                                                                                                tCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCt                                                                   <  1:1024916/102‑1 (MQ=255)
                                                                                                                                tCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGc                   >  1:3111118/1‑150 (MQ=255)
                                                                                                                                 ccAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCa                  <  1:1500605/150‑1 (MQ=255)
                                                                                                                                 ccAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCa                  >  2:2410954/1‑150 (MQ=255)
                                                                                                                                  cAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAg                 >  2:3779916/1‑150 (MQ=255)
                                                                                                                                     cATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTc              >  2:176061/1‑150 (MQ=255)
                                                                                                                                       tCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGc            <  2:253102/150‑1 (MQ=255)
                                                                                                                                          aatCCAGCGGCACGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAc         >  2:4955475/3‑150 (MQ=255)
                                                                                                                                           gTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTg                                                          >  2:3559824/1‑100 (MQ=255)
                                                                                                                                           gTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAcc        >  2:4096988/1‑150 (MQ=255)
                                                                                                                                           gTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAcc        >  2:5142591/1‑150 (MQ=255)
                                                                                                                                           atCCAGCGGCACGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAcc        >  2:881104/2‑150 (MQ=255)
                                                                                                                                            tCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAccc       <  1:3051922/150‑1 (MQ=255)
                                                                                                                                            tCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAccc       <  1:3051545/150‑1 (MQ=255)
                                                                                                                                            tCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAccc       >  1:243299/1‑150 (MQ=255)
                                                                                                                                            tCCAGCGGCACGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCCTGCCTTGGCGCAGCCGCACCACCGCCTGCTCGTGCAGGGCGCCAccc       >  2:4330059/1‑150 (MQ=255)
                                                                                                                                                 cGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCACCCGTTAc  <  2:2657275/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCACCCGTTAC  >  NC_002947/739450‑739744

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