Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,864,709 T→A 5.7% E225V (GAA→GTA)  gllP ← porin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,864,7090TA 5.7% 218.8 / 2.4 70E225V (GAA→GTA) gllPporin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base A (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

CAGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTGTTCT  >  NC_002947/2864561‑2864856
                                                                                                                                                     |                                                                                                                                                    
cAGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCtt                                                                                                                                                      <  8:12226/149‑1 (MQ=255)
 aGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTc                                                                                                                                                     >  4:162223/1‑149 (MQ=255)
        aaGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTgg                                                                                                                                              <  5:367822/149‑1 (MQ=255)
         aGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGACCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTggg                                                                                                                                             <  8:55720/149‑1 (MQ=255)
         aGCGGTGCGGTTGTTAAGTTCGCCTGCCTTGGGGGCACCGTCGTCCTTGCCCCACTAGTATCCCAGGTTTGCGCCCAGTACCCAGTTGCCCTTGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGAAGACATCTACCAGcagc                                                                                                                                              <  1:106744/149‑4 (MQ=255)
              tgcggtTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATAc                                                                                                                                        <  2:313360/149‑1 (MQ=255)
                  gttgttCAGTTCGCCTGCCTTGGCTCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTGGGCGCCCAGTACCCAGTTGCACAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATAc                                                                                                                                        <  7:181160/145‑1 (MQ=255)
                  gttgttCAGTTCGCCTGCCTTGGCGCCACAGTCGTCCTTGCCCCAGAAGTAGACCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCAGGTAGATGTATTCCAGCTGGGCATaa                                                                                                                                        >  8:181160/1‑144 (MQ=255)
                      ttCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGcc                                                                                                                                >  6:319093/1‑149 (MQ=255)
                      ttCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGGAGATGTCTTCCAGCTGGGCATACCAGGCGcc                                                                                                                                >  6:313831/1‑149 (MQ=255)
                       tCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGAAGACATCTACCAGcagcgct                                                                                                                                           >  3:351781/1‑131 (MQ=39)
                       tCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGAAGACATCTACCAGcagcgct                                                                                                                                           <  4:351781/137‑7 (MQ=39)
                         aGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGAc                                                                                                                             >  7:35678/1‑149 (MQ=255)
                          gTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTTCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGAcc                                                                                                                            <  4:209803/149‑1 (MQ=255)
                             cGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATCCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCCGGCGCCGACCAgg                                                                                                                         >  6:384982/1‑149 (MQ=255)
                              gcctgcctTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAgg                                                                                                                         <  5:282816/148‑1 (MQ=255)
                                      tGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTcc                                                                                                                >  4:360964/1‑149 (MQ=255)
                                            cACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGtt                                                                                                           >  8:14722/1‑148 (MQ=255)
                                             aCCGTCGTCCTTGCCCCAGAAGTAGCCCATGTTGGCGCCCAGTACCCAGTTTCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGa                                                                                                         <  4:37995/149‑1 (MQ=255)
                                                   gtcCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCCGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCgg                                                                                                    <  6:56710/148‑1 (MQ=255)
                                                           cccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccg                                                                                           >  5:44767/1‑149 (MQ=255)
                                                            ccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccgc                                                                                          >  7:362635/1‑149 (MQ=255)
                                                            ccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccgc                                                                                          >  7:272076/1‑149 (MQ=255)
                                                            ccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACATGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccgc                                                                                          >  5:200490/1‑149 (MQ=255)
                                                              aGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGccgccg                                                                                        <  1:108046/149‑1 (MQ=255)
                                                                 aGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGc                                                                                      >  1:204698/1‑148 (MQ=255)
                                                                 aGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCg                                                                                     <  8:35678/149‑1 (MQ=255)
                                                                    aGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAg                                                                                 >  4:301212/1‑149 (MQ=255)
                                                                          ggTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGaaa                                                                           <  6:200490/149‑1 (MQ=255)
                                                                               gcgcCCAGTACCCAGTTTCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGTTc                                                                      <  6:217565/149‑1 (MQ=255)
                                                                                cgcCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTcc                                                                     >  8:82594/1‑149 (MQ=255)
                                                                                    cAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAgg                                                                 >  2:300672/1‑149 (MQ=255)
                                                                                      gTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTg                                                               <  8:32968/149‑1 (MQ=255)
                                                                                          ccAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAgg                                                           >  8:173372/1‑149 (MQ=255)
                                                                                              ttGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCtgt                                                       <  8:272076/149‑1 (MQ=255)
                                                                                                gCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGCTCCGAGGTGAAGGCTGTgc                                                     >  8:271650/1‑149 (MQ=255)
                                                                                                 cccAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAg                                                                  <  1:69511/135‑1 (MQ=255)
                                                                                                 cccAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAg                                                                  >  2:69511/1‑135 (MQ=255)
                                                                                                  ccAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTgcgc                                                   <  1:236968/149‑1 (MQ=255)
                                                                                                   cAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGcgcg                                                  <  7:173372/149‑1 (MQ=255)
                                                                                                     aGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGt                                                 >  2:211884/1‑149 (MQ=255)
                                                                                                               gATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCt                                        >  4:276738/1‑148 (MQ=255)
                                                                                                                 tgctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCAt                                     >  3:82567/1‑149 (MQ=255)
                                                                                                                  gctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTcc                                                                     <  4:111059/116‑1 (MQ=255)
                                                                                                                  gctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTcc                                                                     >  3:111059/1‑116 (MQ=255)
                                                                                                                  gctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATg                                    >  2:204346/1‑149 (MQ=255)
                                                                                                                    tgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTc                                  >  4:336996/1‑149 (MQ=255)
                                                                                                                     gcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGTTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTcc                                 >  1:265645/1‑149 (MQ=255)
                                                                                                                     gcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTcc                                 >  6:172493/1‑149 (MQ=255)
                                                                                                                      cACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCt                                <  2:204698/149‑1 (MQ=255)
                                                                                                                       aCCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTg                               <  2:222992/149‑1 (MQ=255)
                                                                                                                         cTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCa                             >  3:9118/1‑149 (MQ=255)
                                                                                                                          tGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCAt                            >  5:277180/1‑149 (MQ=255)
                                                                                                                           ggAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATg                           >  3:6399/1‑149 (MQ=255)
                                                                                                                             aaGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCt                         <  3:276738/149‑1 (MQ=255)
                                                                                                                               gTAGCGCGGCTGGTAGTTGTCTTCCTGGTGGGCATACCAGGGGCCGACCAGGGTGGGCTCCTGGTTGAAGGGGTATTCGCCGCCGGCGAATCTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTg                        <  1:300672/149‑1 (MQ=255)
                                                                                                                                tAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGc                      <  7:14722/149‑1 (MQ=255)
                                                                                                                                tAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGc                      <  4:6399/149‑1 (MQ=255)
                                                                                                                                   cgcTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtc                   <  6:44767/149‑1 (MQ=255)
                                                                                                                                    gctgctGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtcg                  >  6:87729/1‑149 (MQ=255)
                                                                                                                                     ctgctgGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtcgt                 <  3:360964/149‑1 (MQ=255)
                                                                                                                                     cagctgGAAGACATCTACCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtcgt                 >  3:208543/3‑149 (MQ=255)
                                                                                                                                          ggTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTgcgc            >  3:314846/1‑149 (MQ=255)
                                                                                                                                           gTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGcgcg           <  2:265645/149‑1 (MQ=255)
                                                                                                                                              gATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGATGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGc        <  5:384982/149‑1 (MQ=255)
                                                                                                                                              gATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGc        <  1:211884/149‑1 (MQ=255)
                                                                                                                                                tGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTg      <  6:277180/149‑1 (MQ=255)
                                                                                                                                                 gTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTGt     >  8:61992/1‑149 (MQ=255)
                                                                                                                                                   cTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGCGGTGAAGGCTGTGCGCGTGAGCAGGCTCCTGTCCTGCATGCTGGCGGCGTTGCGGGGGcgggtc   >  4:260346/1‑144 (MQ=255)
                                                                                                                                                   cTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTg               >  7:300659/1‑137 (MQ=255)
                                                                                                                                                   cTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTg               <  8:300659/137‑1 (MQ=255)
                                                                                                                                                    ttCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTGTTCt  >  1:376880/1‑149 (MQ=255)
                                                                                                                                                    ttCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTGTTCt  >  7:274760/1‑149 (MQ=255)
                                                                                                                                                    ttCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTGTTCt  >  5:36634/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCCAGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAAGTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTGTTCT  >  NC_002947/2864561‑2864856

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: