Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 739,594 T→G 5.2% T345P (ACG→CCG)  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,5940TG 5.2% 305.4 / 5.0 96T345P (ACG→CCG) pilCtype IV fimbrial assembly protein PilC
Reads supporting (aligned to +/- strand):  ref base T (42/49);  new base G (3/2);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

GCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCACCCGTTACC  >  NC_002947/739447‑739742
                                                                                                                                                     |                                                                                                                                                    
gCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTg                                                                                                                                                     <  1:163148/149‑1 (MQ=255)
gCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTg                                                                                                                                                     >  4:366689/1‑149 (MQ=255)
 cAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGc                                                                                                                                                    <  3:366689/149‑1 (MQ=255)
      aCATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTg                                                                                                                                               >  3:220370/1‑149 (MQ=255)
        aTCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGa                                                                                                                                             >  5:105863/1‑149 (MQ=255)
        aTCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTCCAGCATCTGGTCCAGCGGCACGCTGGac                                                                                                                                            <  2:368820/149‑2 (MQ=255)
         tCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGAt                                                                                                                                            <  8:50672/149‑1 (MQ=255)
             caccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCAc                                                                                                                                        <  3:148325/149‑1 (MQ=255)
             caccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCAc                                                                                                                                        <  1:102484/149‑1 (MQ=255)
             caccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCAc                                                                                                                                        >  3:334742/1‑149 (MQ=255)
               ccaccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGGCACGCTGGaccagatg                                                                                                                                      <  5:150917/149‑8 (MQ=255)
                caccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGa                                                                                                                                     >  3:317389/1‑149 (MQ=255)
                  ccaGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATg                                                                                                                                   <  6:67472/149‑1 (MQ=255)
                   caGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCt                                                                                                                                                >  3:373486/1‑135 (MQ=255)
                   caGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCt                                                                                                                                                <  4:373486/135‑1 (MQ=255)
                   caGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATgg                                                                                                                                  >  4:31835/1‑149 (MQ=255)
                    aGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGc                                                                                                                                 <  2:351436/149‑1 (MQ=255)
                     ggCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGcc                                                                                                                                <  2:15357/149‑1 (MQ=255)
                     ggCCACCCACCACCAATCCCAGGATCAGCACGATGGCCGGCTCCAGCTGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGCTCCAGCGTGCCGCTGGATTCACCGATGGc                                                                                                                                 <  1:219755/149‑1 (MQ=255)
                        cacccaccAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTc                                                                                                                             <  6:174574/149‑1 (MQ=255)
                         acccaccAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCAGTc                                                                                                                             >  2:24302/1‑148 (MQ=255)
                           ccaccaGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAAc                                                                                                                          >  7:84878/1‑149 (MQ=255)
                           ccaccaGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAAc                                                                                                                          >  5:12831/1‑149 (MQ=255)
                           ccaccaGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAAc                                                                                                                          >  3:294498/1‑149 (MQ=255)
                                aGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCAc                                                                                                                     <  4:317389/149‑1 (MQ=255)
                                       ccAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAg                                                                                                              >  1:379403/1‑149 (MQ=255)
                                        cAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgg                                                                                                             >  1:46923/1‑149 (MQ=255)
                                         aGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAggg                                                                                                            <  3:31835/149‑1 (MQ=255)
                                           gATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAggggg                                                                                                          <  1:191072/149‑1 (MQ=255)
                                               aGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGaa                                                                                                      >  6:48617/1‑149 (MQ=255)
                                                gCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGaa                                                                                                      >  3:289491/1‑147 (MQ=255)
                                                    gATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCg                                                                                                 >  8:328731/1‑149 (MQ=255)
                                                    gATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCg                                                                                                 >  4:195881/1‑149 (MQ=255)
                                                      tGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTTGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCt                                                                                                                                                <  6:89363/100‑1 (MQ=255)
                                                      tGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTCGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCt                                                                                                                                                >  5:89363/1‑100 (MQ=255)
                                                                        tGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgggg                                                                                                           >  4:230524/1‑119 (MQ=255)
                                                                        tGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgggg                                                                                                           <  3:230524/119‑1 (MQ=255)
                                                                        tGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTgg                                                                             >  5:204291/1‑149 (MQ=255)
                                                                          gTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGtt                                                                           >  8:213411/1‑149 (MQ=255)
                                                                           tCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGtt                                                                           >  6:16411/1‑148 (MQ=255)
                                                                           tCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGtt                                                                           >  1:244012/1‑148 (MQ=255)
                                                                             aaCTGATCCATCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCa                                                                         <  1:361073/148‑1 (MQ=255)
                                                                             aaCTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgggg                                                                                                           >  1:376736/1‑114 (MQ=255)
                                                                             aaCTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAgggg                                                                                                           <  2:376736/114‑1 (MQ=255)
                                                                              aCTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCTTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGc                                                                       <  1:105052/149‑1 (MQ=255)
                                                                              aCTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGc                                                                       >  2:186302/1‑149 (MQ=255)
                                                                                tGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGccc                                                                     >  1:277534/1‑149 (MQ=255)
                                                                                  aTCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGGGGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTCGTTCAGCCCt                                                                    >  1:302086/1‑149 (MQ=255)
                                                                                  aTCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGccc                                                                     <  6:63396/147‑1 (MQ=255)
                                                                                  aTCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCtt                                                                   <  4:294498/149‑1 (MQ=255)
                                                                                  aTCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCtt                                                                   <  2:379403/149‑1 (MQ=255)
                                                                                   tCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTg                                                                  <  1:167194/149‑1 (MQ=255)
                                                                                     cAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGcc                                                                >  2:143861/1‑149 (MQ=255)
                                                                                     cAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGcc                                                                <  7:328731/149‑1 (MQ=255)
                                                                                          ccTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTg                                                           <  2:170293/149‑1 (MQ=255)
                                                                                           cTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTgg                                                          >  3:324274/1‑149 (MQ=255)
                                                                                            tGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGc                                                         >  3:70657/1‑149 (MQ=255)
                                                                                             ggCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGc                                                         >  3:294283/1‑148 (MQ=255)
                                                                                              gCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGc                                                         <  4:70657/147‑1 (MQ=255)
                                                                                               cTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCAt                                                      <  5:191368/149‑1 (MQ=255)
                                                                                                       ttCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTgg                                               <  3:195881/148‑1 (MQ=255)
                                                                                                       ttCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgc                                              >  4:340453/1‑149 (MQ=255)
                                                                                                         ccTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgcg                                             >  5:137941/1‑148 (MQ=255)
                                                                                                          cTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGgcg                                             >  1:6400/1‑147 (MQ=255)
                                                                                                             gTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCa                                                                         <  4:132042/116‑1 (MQ=255)
                                                                                                             gTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAg                                          <  2:6400/147‑1 (MQ=255)
                                                                                                              tAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCg                                       >  2:267360/1‑149 (MQ=255)
                                                                                                                aTGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGc                                      <  8:84878/148‑1 (MQ=255)
                                                                                                                aTGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGc                                      <  2:46923/148‑1 (MQ=255)
                                                                                                                 tggctggcGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCaa                                    <  4:334742/149‑1 (MQ=255)
                                                                                                                   gctggcGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAaca                                  <  5:48617/149‑1 (MQ=255)
                                                                                                                       gcggcTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGc                              >  7:196247/1‑149 (MQ=255)
                                                                                                                        cggcTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGcc                             >  1:240401/1‑149 (MQ=255)
                                                                                                                           cTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGc                          <  2:244012/149‑1 (MQ=255)
                                                                                                                              tttCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTc                        <  6:137941/148‑1 (MQ=255)
                                                                                                                                  cAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTg                     >  1:173813/1‑147 (MQ=255)
                                                                                                                                  cAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGc                    <  2:302086/148‑1 (MQ=255)
                                                                                                                                  cAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGc                    <  8:196247/148‑1 (MQ=255)
                                                                                                                                  cAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGc                    <  5:307216/148‑1 (MQ=255)
                                                                                                                                  cAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCa                   <  4:289491/149‑1 (MQ=255)
                                                                                                                                    gCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAgg                 <  4:324274/149‑1 (MQ=255)
                                                                                                                                     cATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGt                <  2:277534/149‑1 (MQ=255)
                                                                                                                                     cATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTACAGGt                <  1:143861/149‑1 (MQ=255)
                                                                                                                                      aTCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTc               <  8:255311/149‑1 (MQ=255)
                                                                                                                                       tCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCg              >  5:59293/1‑149 (MQ=255)
                                                                                                                                        gtgaatCCAGCGGCACGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGAAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCg              >  4:162230/6‑149 (MQ=255)
                                                                                                                                         tGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGc             >  5:189720/1‑148 (MQ=255)
                                                                                                                                           gTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCg              <  7:213411/145‑1 (MQ=255)
                                                                                                                                           gTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCa           <  1:267360/148‑1 (MQ=255)
                                                                                                                                            atCCAGCGGCACGCTGGATTCACCGATGGCCGTCAACTGCTCCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAc          >  7:251197/2‑149 (MQ=255)
                                                                                                                                             tCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGAGCAGCCGAAACACCGCCTGCTCTGGCAGGTCGACAcc         >  7:337960/1‑149 (MQ=255)
                                                                                                                                               cAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAccc        <  4:294283/148‑1 (MQ=255)
                                                                                                                                               cAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCAccc        <  6:204291/148‑1 (MQ=255)
                                                                                                                                               cAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCACCCg       >  2:254710/1‑149 (MQ=255)
                                                                                                                                                  cGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCa                                                                         >  5:141310/1‑80 (MQ=255)
                                                                                                                                                  cGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCa                                                                         <  6:141310/80‑1 (MQ=255)
                                                                                                                                                     tGCCGCTGTATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCACCCGTTAcc  <  4:220370/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCCTGGCTGACCTGTTCCTCGTAATGGCTGGCGGCTTTTTCCAGCATCTGGTCCAGCGTGCCGCTGGATTCACCGATGGCCGTCAACTGCACCAGCAAGGGGGGGAACAGCGGTTCACCGGCCATCGCCTGGTTCAGCCCTTGCCCATTGGCCATGCCTTGGCGCAGCCGCAACACCGCCTGCTCGTGCAGGTCGCCACCCGTTACC  >  NC_002947/739447‑739742

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