Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,041,198 G→T 5.4% A46D (GCT→GAT)  PP_0900 ← PAP2 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,041,1980GT 5.4% 238.5 / 2.5 74A46D (GCT→GAT) PP_0900PAP2 family protein
Reads supporting (aligned to +/- strand):  ref base G (31/39);  new base T (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

ACGGCGCCAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGACGGCTCCAGG  >  NC_002947/1041050‑1041346
                                                                                                                                                     |                                                                                                                                                    
aCGGCGCCAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAg                                                                                                                                                      >  1:158450/1‑149 (MQ=255)
aCGGCGCCAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAg                                                                                                                                                      >  3:272475/1‑149 (MQ=255)
aCGGCGCCAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAg                                                                                                                                                      >  4:232535/1‑149 (MQ=255)
  ggCGCTAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCg                                                                                                                                                    >  3:95645/1‑149 (MQ=255)
  ggCGCCAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCTGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCg                                                                                                                                                    <  3:391873/149‑1 (MQ=255)
    cgcCAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCggg                                                                                                                                                  <  8:32212/149‑1 (MQ=255)
       cagccaAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTc                                                                                                                                               <  5:295621/149‑1 (MQ=255)
       cagccaAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTc                                                                                                                                               >  5:65047/1‑149 (MQ=255)
        agccaAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTc                                                                                                                                               >  6:32266/1‑148 (MQ=255)
             aaCAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCgaagaa                                                                                                                                         <  4:272475/149‑1 (MQ=255)
                aGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAg                                                                                                                                      >  3:156472/1‑149 (MQ=255)
                        ccAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGc                                                                                                                                <  4:95645/147‑1 (MQ=255)
                           gAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGAGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACa                                                                                                                            >  6:116275/1‑148 (MQ=255)
                           gAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGAGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACa                                                                                                                            <  5:116275/148‑1 (MQ=255)
                            aaaaaGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATc                                                                                                                          >  2:177044/1‑149 (MQ=255)
                                   tCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTc                                                                                                                    >  6:96033/1‑148 (MQ=255)
                                    cgcgGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCg                                                                                                                                                    <  7:117314/115‑1 (MQ=255)
                                    cgcgGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCg                                                                                                                                                    >  8:117314/1‑115 (MQ=255)
                                    cgcgGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTTGa                                                                                                                  <  6:62718/149‑1 (MQ=255)
                                     gcgGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCAGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGAc                                                                                                                 >  1:104605/1‑149 (MQ=255)
                                       ggCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCgaag                                                                                                                                           <  1:45726/121‑1 (MQ=255)
                                       ggCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCgaag                                                                                                                                           >  2:45726/1‑121 (MQ=255)
                                         cAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGaa                                                                                                             <  5:96033/149‑1 (MQ=255)
                                         cAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGaa                                                                                                             >  3:56030/1‑149 (MQ=255)
                                           aaGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGt                                                                                                           <  8:16072/149‑1 (MQ=255)
                                             ggTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCa                                                                                                         <  5:91296/149‑1 (MQ=255)
                                                 agcagcGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCt                                                                                                      <  1:237790/148‑1 (MQ=255)
                                                           tCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGAt                                                                                            >  5:407368/1‑148 (MQ=255)
                                                             aCCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGACAGCAACAGAAGGaac                                                                                           >  8:57082/1‑145 (MQ=255)
                                                              cccccAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAg                                                                                        <  5:12765/149‑1 (MQ=255)
                                                              cccccAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAg                                                                                        >  4:29831/1‑149 (MQ=255)
                                                                 ccAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGc                                                                                      >  4:208253/1‑148 (MQ=255)
                                                                   aaCAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCAc                                                                                   <  2:104605/149‑1 (MQ=255)
                                                                    aCAGGATCACCACTTGCTTGACCCGGTCATGCAGGATTTTTTCCAACAGGTAACTGTGTCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCATGATCAAGGCCAcc                                                                                  <  1:276686/149‑1 (MQ=255)
                                                                            accacTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGc                                                                          <  2:184905/149‑1 (MQ=255)
                                                                             ccacTTGCTTGACCCGGTCATGCAGGATGTGTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGc                                                                          >  1:242161/1‑149 (MQ=255)
                                                                             ccacTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCATCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGc                                                                          >  1:265911/1‑148 (MQ=255)
                                                                             ccacTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCATCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGc                                                                          >  1:271584/1‑148 (MQ=255)
                                                                              cacTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGc                                                                          >  7:34502/1‑147 (MQ=255)
                                                                              cacTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGAAGTCga                                                                          >  4:220302/1‑146 (MQ=255)
                                                                                 ttgcttgACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAggg                                                                     >  1:326684/1‑149 (MQ=255)
                                                                                   gcttgACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAggg                                                                     <  5:32266/147‑1 (MQ=255)
                                                                                     ttgACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCACGGCCACCGCCGTCGCCGGGGGact                                                                 >  6:283789/1‑147 (MQ=255)
                                                                                      tgACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATc                                                                >  2:273342/1‑149 (MQ=255)
                                                                                       gACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATc                                                                <  3:208253/148‑1 (MQ=255)
                                                                                        aCCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATccc                                                              >  1:358548/1‑149 (MQ=255)
                                                                                           cGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAg                                                           <  4:125819/149‑1 (MQ=255)
                                                                                              tCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTa                                                        <  6:154293/149‑1 (MQ=255)
                                                                                                  gCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGc                                                    >  4:264206/1‑149 (MQ=255)
                                                                                                   cAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGc                                                    <  3:232535/148‑1 (MQ=255)
                                                                                                    aGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGc                                                    >  2:322500/1‑147 (MQ=255)
                                                                                                    aGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGc                                                    >  1:248288/1‑147 (MQ=255)
                                                                                                    aGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCa                                                  <  1:177044/149‑1 (MQ=255)
                                                                                                     ggATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCATCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCag                                                 <  2:265911/149‑1 (MQ=255)
                                                                                                     ggATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCATCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCag                                                 <  2:271584/149‑1 (MQ=255)
                                                                                                     ggATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGTGGAATCCCCAGGTACAGCCag                                                 <  3:29831/149‑1 (MQ=255)
                                                                                                     ggATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCag                                                 >  8:332209/1‑149 (MQ=255)
                                                                                                       aTGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCagta                                               <  7:57082/149‑1 (MQ=255)
                                                                                                       aTGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCagta                                               <  2:326684/149‑1 (MQ=255)
                                                                                                         gTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGt                                                                                                           >  5:129680/1‑87 (MQ=255)
                                                                                                         gTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGt                                                                                                           <  6:129680/87‑1 (MQ=255)
                                                                                                                     ggTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGa                                  <  6:417474/149‑1 (MQ=255)
                                                                                                                             tgGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGt                          <  6:407368/149‑1 (MQ=255)
                                                                                                                             tgGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGt                          <  4:56030/149‑1 (MQ=255)
                                                                                                                             gtgcggcCAATGAACTGCCCGGCAGCGGGGTCGAAGAACATGTTGGCGACATCCATGTCGACTGAAGGCAGCTCCAGCAGGATCGAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTTCAGCCAGTAGTTGATCGGGCGAGGGCgg                            <  6:65047/145‑1 (MQ=255)
                                                                                                                                    cAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACa                                                      <  3:219200/114‑1 (MQ=255)
                                                                                                                                    cAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACa                                                      >  4:219200/1‑114 (MQ=255)
                                                                                                                                     aaTGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTtgct                  <  1:273342/149‑1 (MQ=255)
                                                                                                                                     aaTGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTtgct                  <  2:248288/149‑1 (MQ=255)
                                                                                                                                     aaTGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTtgct                  <  4:156472/149‑1 (MQ=255)
                                                                                                                                         aaCTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCAt              <  8:34502/149‑1 (MQ=255)
                                                                                                                                            tGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGAc           >  7:1198/1‑149 (MQ=255)
                                                                                                                                               ccGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGACGGc        >  2:150841/1‑149 (MQ=255)
                                                                                                                                               ccGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGACGGc        <  7:332209/149‑1 (MQ=255)
                                                                                                                                                 ggCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGACGGCTc      <  1:322500/149‑1 (MQ=255)
                                                                                                                                                    aGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAg                                                           >  6:310107/1‑93 (MQ=255)
                                                                                                                                                    aGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAg                                                           <  5:310107/93‑1 (MQ=255)
                                                                                                                                                     gCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCTAGGCCACCGCAGTCGCCAGGGGATTGCCGAGGTACAGCGAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGACGGCTCCAgg  <  5:283789/149‑1 (MQ=255)
                                                                                                                                                     gCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGACGGCTCCAgg  <  3:264206/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGGCGCCAGCCAAACAGGCGCTTCCAGAAAAAGCTCGCGGCAAAGGTCAGCAGCGATATCACCCCCAACAGGATCACCACTTGCTTGACCCGGTCATGCAGGATGTTTTCCAACAGGTAACTGTGGCGGCCAATGAACTGCCCGGCAGCGGGGTCGAAGAACAGGTTGGCGACATCCATGTCGACTGAAGTCAGCTCCAGCAGGATCAAGGCCACCGCAGTCGCCAGGGGAATCCCCAGGTACAGCCAGTAGTTGATCGGGCGAGGGCGGGTGCGTTGCTGCATGACGGCTCCAGG  >  NC_002947/1041050‑1041346

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: