Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,921,695 A→G 5.9% I702V (ATC→GTC) ‡ PP_4330 → toxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,921,6950AG 5.9% 198.3 / 2.1 68I702V (ATC→GTC) ‡PP_4330toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAACGATCGATCAGGTCGAGCTGCG  >  NC_002947/4921550‑4921831
                                                                                                                                                 |                                                                                                                                         
cGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCg                                                                                                                                        >  8:151437/1‑149 (MQ=255)
cGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCg                                                                                                                                        >  3:242639/1‑149 (MQ=255)
     tcGATCAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGATGCAGCGCTGGATGGCAACGGTGGCATCGccac                                                                                                                                    <  4:55981/147‑1 (MQ=255)
             aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                                                                                                                           >  4:253648/1‑149 (MQ=255)
             aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                                                                                                                           <  7:151437/149‑1 (MQ=255)
             aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                                                                                                                           <  4:317242/149‑1 (MQ=255)
                gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                                                                                                                             <  5:10274/144‑1 (MQ=255)
                gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                                                                                                                             >  6:10274/1‑144 (MQ=255)
                gCAAGGCTATTGGCAGCCCGGGCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACACGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                                                                                                                       >  3:229079/1‑149 (MQ=255)
                    ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGc                                                                                                                     >  7:105799/1‑148 (MQ=255)
                      cTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGAc                                                                                                                  <  4:242639/149‑1 (MQ=255)
                       tATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTATGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgt                                                                                                                                 <  8:302009/133‑3 (MQ=255)
                       tATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTATGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgt                                                                                                                                 >  7:302009/1‑131 (MQ=255)
                           tgCAGCCCTGGGGCCCGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtg                                                                                                             >  2:244674/2‑149 (MQ=255)
                           ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt                                                                                                              <  8:76747/148‑1 (MQ=255)
                           ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt                                                                                                              <  3:253648/148‑1 (MQ=255)
                            gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt                                                                                                            >  6:142000/1‑149 (MQ=255)
                            gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt                                                                                                            >  7:260128/1‑149 (MQ=255)
                              aGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgc                                                                                                         >  5:2033/1‑149 (MQ=255)
                               gCCCGGGCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACACGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg                                                                                                       <  4:229079/149‑1 (MQ=255)
                                cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg                                                                                                       >  6:93769/1‑149 (MQ=255)
                                     gCCCCAGGGCTGCAGGTTTTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTa                                                                                                  >  2:241464/1‑149 (MQ=255)
                                       cccAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTAcc                                                                                                >  5:247566/1‑149 (MQ=255)
                                         cAGGGCTGCAGGTATTGCCGCTTGAAGAGCAATGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCgg                                                                                              <  6:247566/149‑1 (MQ=255)
                                           gggCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGAt                                                                                            <  5:142000/149‑1 (MQ=255)
                                           gggCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGAt                                                                                            <  7:213677/149‑1 (MQ=255)
                                                       ttGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCAt                                                                                  >  7:65633/1‑147 (MQ=255)
                                                       ttGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGa                                                                                <  5:93769/149‑1 (MQ=255)
                                                        tGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCTCCACCAATTACGAGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGa                                                                                <  5:252944/148‑1 (MQ=255)
                                                         gCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTCGCTGCGCTACCGGATGCTGCACCATGa                                                                                <  1:14948/148‑1 (MQ=255)
                                                          ccGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAggg                                                                             <  8:260128/149‑1 (MQ=255)
                                                            gCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTg                                                                           >  2:260340/1‑149 (MQ=255)
                                                              ttGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGcca                                                                                                                                     <  6:148695/90‑1 (MQ=255)
                                                              ttGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGcca                                                                                                                                     >  5:148695/1‑90 (MQ=255)
                                                              ttGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTggg                                                                         <  7:317195/149‑1 (MQ=255)
                                                                gAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                                                                                                                             >  6:220203/1‑96 (MQ=255)
                                                                gAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                                                                                                                             <  5:220203/96‑1 (MQ=255)
                                                                gAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTggg                                                                         >  8:317195/1‑147 (MQ=255)
                                                                gAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGc                                                                        <  8:105799/148‑1 (MQ=255)
                                                                  agagAAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGAGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTg                                                                      >  5:281041/1‑148 (MQ=255)
                                                                   gagCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTa                                                                    >  3:178250/1‑149 (MQ=255)
                                                                      cAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgtt                                                                                                                                >  5:156727/1‑84 (MQ=255)
                                                                      cAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgtt                                                                                                                                <  6:156727/87‑4 (MQ=255)
                                                                                      ccTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTc                                                 <  4:178250/149‑1 (MQ=255)
                                                                                          cTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCct                                               <  6:2033/147‑1 (MQ=255)
                                                                                           tGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCctgct                                            >  5:220003/1‑149 (MQ=255)
                                                                                             ccAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTAc                                                                   <  6:62748/124‑1 (MQ=255)
                                                                                             ccAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTAc                                                                   >  5:62748/1‑124 (MQ=255)
                                                                                             ccAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTAc                                                                   >  7:199367/1‑124 (MQ=255)
                                                                                             ccAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTAc                                                                   <  8:199367/124‑1 (MQ=255)
                                                                                                   ccGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACg                                    >  4:18503/1‑149 (MQ=255)
                                                                                                     gCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCct                                               <  4:87086/136‑1 (MQ=255)
                                                                                                     gCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCct                                               >  3:87086/1‑136 (MQ=255)
                                                                                                      cTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGACTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGc                                  >  1:8869/1‑148 (MQ=255)
                                                                                                              cgcAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACg                                    >  3:119482/1‑138 (MQ=255)
                                                                                                              cgcAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACg                                    <  4:119482/138‑1 (MQ=255)
                                                                                                               gcAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGaaa                        >  6:150712/1‑149 (MQ=255)
                                                                                                                cAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAAc                       <  1:260340/149‑1 (MQ=255)
                                                                                                                   tACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAAcgat                    <  8:65633/149‑1 (MQ=255)
                                                                                                                      gAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAAcgatcga                 >  5:108704/1‑149 (MQ=255)
                                                                                                                        ggCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAACgatcgatc               <  3:18503/149‑1 (MQ=255)
                                                                                                                                tGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAACGATCGATCAGGTCGAg       >  7:260138/1‑149 (MQ=255)
                                                                                                                                  gATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAACGATCGATCAGGTCGAGc      <  8:260138/148‑1 (MQ=255)
                                                                                                                                      gCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAACGATCGATCAGGTCGAGCTgcg  <  6:281041/148‑1 (MQ=255)
                                                                                                                                       cAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGgacga                                                     >  2:148432/1‑96 (MQ=255)
                                                                                                                                       cAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGgacga                                                     <  1:148432/96‑1 (MQ=255)
                                                                                                                                        aaCGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAggg                             <  2:141112/119‑1 (MQ=255)
                                                                                                                                        aaCGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAggg                             >  1:141112/1‑119 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
CGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGATGCTGCACCATGAGGGTGGGCTGTACATGGACGTGGACGACTTCCTGCTGCCGGACGGCAAGGGCGAAACGATCGATCAGGTCGAGCTGCG  >  NC_002947/4921550‑4921831

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: