Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,921,588 C→G 7.2% A666G (GCC→GGC) ‡ PP_4330 → toxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,921,5880CG 7.2% 224.3 / 6.7 71A666G (GCC→GGC) ‡PP_4330toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (3/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTA  >  NC_002947/4921446‑4921735
                                                                                                                                              |                                                                                                                                                   
aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                               >  7:76747/1‑149 (MQ=255)
          gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTa                                                                                                                                     <  4:56372/149‑1 (MQ=255)
          gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTa                                                                                                                                     >  8:89364/1‑149 (MQ=255)
          gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGATTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGccaata                                                                                                                                     <  1:112042/149‑6 (MQ=255)
                    gCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                            <  5:313325/148‑1 (MQ=255)
                     cTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTg                                                                                                                           <  5:182275/148‑1 (MQ=255)
                      tGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGc                                                                                                                                          >  8:294110/1‑132 (MQ=255)
                          cAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAagag                                                                                                                      >  2:13357/1‑148 (MQ=255)
                                 gccgccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGc                                                                                                                <  6:76247/147‑1 (MQ=255)
                                  ccgccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGACTGTTGCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGctt                                                                                                              <  5:203523/148‑1 (MQ=255)
                                    gccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                           <  7:89364/149‑1 (MQ=255)
                                    gccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGctt                                                                                                              <  8:331364/146‑1 (MQ=255)
                                    gccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGctt                                                                                                              >  7:331364/1‑146 (MQ=255)
                                     ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttc                                                                                                             >  3:71189/1‑146 (MQ=255)
                                     ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttc                                                                                                             <  4:71189/146‑1 (MQ=255)
                                                cgcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGACGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTg                                                                                               <  1:227454/149‑1 (MQ=255)
                                                     gAGCGACAGTGAATATAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCa                                                                                            >  6:128698/1‑147 (MQ=255)
                                                           cAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCt                                                                                    <  8:20696/149‑1 (MQ=255)
                                                              tGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTAt                                                                                  >  7:231267/1‑148 (MQ=255)
                                                                  tACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTAcgc                                                                             <  1:13357/149‑1 (MQ=255)
                                                                        aTGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCag                                                                                           >  1:257716/1‑129 (MQ=255)
                                                                        aTGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCag                                                                                           <  2:257716/129‑1 (MQ=255)
                                                                            gcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGATCAGATGGCAAGGCTATTGGCAGCCCGGGCCCCATGGCTGCAGGTATTGCCGCTTGAAGAGCaa                                                                                                                   <  5:104753/101‑1 (MQ=255)
                                                                            gcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                   <  8:231267/149‑1 (MQ=255)
                                                                            gcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                   >  5:172418/1‑149 (MQ=255)
                                                                              cTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCAGCAGGTATTGCCGCTAGAAGCGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGACCCTATTAAGCGCAGTACGAgg                                                                  >  3:55981/1‑148 (MQ=255)
                                                                                  ttCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAgcg                                                             >  3:317242/1‑149 (MQ=255)
                                                                                         gAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTAcgcg                                                                            <  6:188838/127‑1 (MQ=255)
                                                                                         gAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTAcgcg                                                                            >  5:188838/1‑127 (MQ=255)
                                                                                         gAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATg                                                      <  6:172418/149‑1 (MQ=255)
                                                                                         gAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATg                                                      <  5:128698/149‑1 (MQ=255)
                                                                                          aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTg                                                                                                                                           <  8:189371/63‑1 (MQ=37)
                                                                                          aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTg                                                                                                                                           >  7:189371/1‑63 (MQ=37)
                                                                                             gCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCaa                                                  >  6:252944/1‑149 (MQ=255)
                                                                                              cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGa                                                        >  5:194364/1‑142 (MQ=255)
                                                                                              cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGa                                                        <  6:194364/142‑1 (MQ=255)
                                                                                               ttCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACg                                                >  8:257981/1‑149 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                   <  7:328637/127‑1 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                   >  8:328637/1‑127 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                   >  4:32579/1‑127 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                   <  3:32579/127‑1 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                             <  7:257981/149‑1 (MQ=255)
                                                                                                    cTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTgg                                            >  8:213677/1‑148 (MQ=255)
                                                                                                        cGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCg                                       >  8:151437/1‑149 (MQ=255)
                                                                                                        cGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCg                                       >  3:242639/1‑149 (MQ=255)
                                                                                                             tcGATCAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGATGCAGCGCTGGATGGCAACGGTGGCATCGccac                                   <  4:55981/147‑1 (MQ=255)
                                                                                                                     aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                             >  1:303864/1‑130 (MQ=255)
                                                                                                                     aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                             <  2:303864/130‑1 (MQ=255)
                                                                                                                     aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                          <  4:317242/149‑1 (MQ=255)
                                                                                                                     aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                          >  4:253648/1‑149 (MQ=255)
                                                                                                                     aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                          <  7:151437/149‑1 (MQ=255)
                                                                                                                        gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                            <  5:10274/144‑1 (MQ=255)
                                                                                                                        gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                            >  6:10274/1‑144 (MQ=255)
                                                                                                                        gCAAGGCTATTGGCAGCCCGGG‑CCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACACGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                      >  3:229079/1‑149 (MQ=255)
                                                                                                                            ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGc                    >  7:105799/1‑148 (MQ=255)
                                                                                                                              cTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGAc                 <  4:242639/149‑1 (MQ=255)
                                                                                                                               tATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGc                                                                                              <  6:106307/71‑1 (MQ=255)
                                                                                                                               tATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGc                                                                                              >  5:106307/1‑71 (MQ=255)
                                                                                                                               tATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTATGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgt                                >  7:302009/1‑131 (MQ=255)
                                                                                                                               tATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTATGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgt                                <  8:302009/133‑3 (MQ=255)
                                                                                                                                   tgCAGCCCTGGGGCCCGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtg            >  2:244674/2‑149 (MQ=255)
                                                                                                                                   ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt             <  3:253648/148‑1 (MQ=255)
                                                                                                                                   ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt             <  8:76747/148‑1 (MQ=255)
                                                                                                                                    gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt           >  7:260128/1‑149 (MQ=255)
                                                                                                                                    gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt           >  6:142000/1‑149 (MQ=255)
                                                                                                                                      aGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgc         >  5:2033/1‑149 (MQ=255)
                                                                                                                                       gCCCGGG‑CCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACACGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg       <  4:229079/149‑1 (MQ=255)
                                                                                                                                        cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAg                                                                         >  3:302158/1‑83 (MQ=255)
                                                                                                                                        cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAg                                                                         <  4:302158/83‑1 (MQ=255)
                                                                                                                                        cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg       >  6:93769/1‑149 (MQ=255)
                                                                                                                                             gCCCCAGGGCTGCAGGTTTTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTa  >  2:241464/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTA  >  NC_002947/4921446‑4921735

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