Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,921,589 C→G 7.0% A666A (GCC→GCG) ‡ PP_4330 → toxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,921,5890CG 7.0% 222.2 / 6.5 71A666A (GCC→GCG) ‡PP_4330toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base C (33/33);  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 = 5.79e-01

AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACC  >  NC_002947/4921446‑4921737
                                                                                                                                               |                                                                                                                                                    
aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                                 >  7:76747/1‑149 (MQ=255)
          gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTa                                                                                                                                       >  8:89364/1‑149 (MQ=255)
          gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTa                                                                                                                                       <  4:56372/149‑1 (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                                                                     >  5:172418/1‑149 (MQ=255)
                                                                            gcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                     <  8:231267/149‑1 (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                                                                                                                                             >  7:189371/1‑63 (MQ=37)
                                                                                          aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTg                                                                                                                                             <  8:189371/63‑1 (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                                                                     <  3:32579/127‑1 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                     >  4:32579/1‑127 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                     >  8:328637/1‑127 (MQ=255)
                                                                                                  gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAg                                                                     <  7:328637/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                            <  7:151437/149‑1 (MQ=255)
                                                                                                                     aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                            <  4:317242/149‑1 (MQ=255)
                                                                                                                     aaCGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                            >  4:253648/1‑149 (MQ=255)
                                                                                                                        gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                              >  6:10274/1‑144 (MQ=255)
                                                                                                                        gCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAc                              <  5:10274/144‑1 (MQ=255)
                                                                                                                        gCAAGGCTATTGGCAGCCCGGGCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACACGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                        >  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               <  8:76747/148‑1 (MQ=255)
                                                                                                                                   ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt               <  3:253648/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)
                                                                                                                                       gCCCGGGCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACACGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg         <  4:229079/149‑1 (MQ=255)
                                                                                                                                        cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAg                                                                           <  4:302158/83‑1 (MQ=255)
                                                                                                                                        cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAg                                                                           >  3:302158/1‑83 (MQ=255)
                                                                                                                                        cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg         >  6:93769/1‑149 (MQ=255)
                                                                                                                                             gCCCCAGGGCTGCAGGTTTTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTa    >  2:241464/1‑149 (MQ=255)
                                                                                                                                               cccAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTAcc  >  5:247566/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACC  >  NC_002947/4921446‑4921737

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: