Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 125,468 C→A 6.6% A96D (GCT→GAT)  PP_0120 → zinc ABC transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947125,4680CA 6.6% 226.6 / 3.4 70A96D (GCT→GAT) PP_0120zinc ABC transporter
Reads supporting (aligned to +/- strand):  ref base C (36/29);  new base A (3/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCCGGCGCATCCCCCCACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATGCGCGAGT  >  NC_002947/125323‑125615
                                                                                                                                                 |                                                                                                                                                   
ccgccCGGCGCATCCCCCCACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCt                                                                                                                                                    <  2:1625/147‑1 (MQ=255)
 cgccCGGCGCATCCCCCCACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCagtg                                                                                                                                                 >  7:319923/1‑149 (MQ=255)
        cgcATCCCCCCACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGc                                                                                                                                          <  2:75849/149‑1 (MQ=255)
          cATCCCCCCACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTg                                                                                                                                        >  7:356685/1‑149 (MQ=255)
              cccccACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCTAGGCCACAGTGGCTGTGCAGTCGCTGCCgg                                                                                                                                    <  1:30256/149‑1 (MQ=255)
               ccccACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCggg                                                                                                                                   >  3:372772/1‑149 (MQ=255)
                 ccaccaCTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCa                                                                                                                                 <  2:99058/149‑1 (MQ=255)
                        tATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCt                                                                                                                          >  5:86075/1‑149 (MQ=255)
                        tATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCt                                                                                                                          >  4:138909/1‑149 (MQ=255)
                           gcgctgcggCCCTCCGTCGTGCGGCGGGTAGGCGACGCCGAGCTGTTGTACTGGTTCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCACTGCCGGGCATGAAGCTgc                                                                                                                        <  8:319923/148‑1 (MQ=255)
                           gcgctgcgcCCCTCCGACGTGCGGCGGGTAGGCGATGCCGAACTCTTGTACTGGATCGGCCCGGACATGGAAAAGTTCCTTGCGCGTGTGCTGGGTAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGACGCTgag                                                                                                                       >  2:286616/1‑147 (MQ=255)
                                   ccccTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCCCAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTg                                                                                                               >  3:264430/1‑149 (MQ=255)
                                    cccTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCTTGAAGCTGCGCCACTTTgg                                                                                                              >  3:231334/1‑149 (MQ=255)
                                    cccTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTgg                                                                                                              >  1:211830/1‑149 (MQ=255)
                                             gTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGa                                                                                                           >  3:16209/1‑143 (MQ=255)
                                             gTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGa                                                                                                           <  4:16209/143‑1 (MQ=255)
                                             gTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAg                                                                                                     >  3:182153/1‑149 (MQ=255)
                                               gcggagGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAACGCCACAGTGGATGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGc                                                                                                    >  7:357299/1‑148 (MQ=255)
                                                  gcggGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCAtt                                                                                                <  7:241950/149‑1 (MQ=255)
                                                      gTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTgcg                                                                                                                       >  7:353560/1‑122 (MQ=255)
                                                      gTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTgcg                                                                                                                       <  8:353560/122‑1 (MQ=255)
                                                      gTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCa                                                                                            >  1:259863/1‑149 (MQ=255)
                                                       tAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCGGTAGCAAGGCCACAGAGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCAt                                                                                           <  6:84901/149‑1 (MQ=255)
                                                       tAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCTTGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCAt                                                                                           <  4:231334/149‑1 (MQ=255)
                                                         ggCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGa                                                                                         <  4:182153/149‑1 (MQ=255)
                                                         ggCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGa                                                                                         <  2:121619/149‑1 (MQ=255)
                                                                ccGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGgaagaa                                                                                  <  6:86075/149‑1 (MQ=255)
                                                                  gACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGa                                                                                                           >  3:170224/1‑122 (MQ=255)
                                                                  gACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGa                                                                                                           <  4:170224/122‑1 (MQ=255)
                                                                       gttgtACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCAt                                                                            >  2:54069/1‑148 (MQ=255)
                                                                          gtACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGAc                                                                         >  2:252324/1‑148 (MQ=255)
                                                                           gaCTGGAGCCGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGa                                                                       <  4:308211/148‑1 (MQ=255)
                                                                            aCTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGAt                                                                      >  3:374192/1‑149 (MQ=255)
                                                                                gATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATg                                                                  <  4:374192/149‑1 (MQ=255)
                                                                                    ggCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGAcc                                                               <  5:352896/148‑1 (MQ=255)
                                                                                      cccGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCAc                                                             >  3:122825/1‑148 (MQ=255)
                                                                                               ggAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATg                                                                  >  4:271529/1‑134 (MQ=255)
                                                                                               ggAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATg                                                                  <  3:271529/134‑1 (MQ=255)
                                                                                               ggAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGcc                                                   >  6:77803/1‑149 (MQ=255)
                                                                                                gAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGga                                                                                      <  4:173628/113‑1 (MQ=255)
                                                                                                gAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGga                                                                                      >  3:173628/1‑113 (MQ=255)
                                                                                                         cTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCt                                         <  3:138909/149‑1 (MQ=255)
                                                                                                               cgTGTGCTGGGAAGCCGTAGCAAGGGCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAATCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCGTGGATgc                                   <  8:357299/149‑1 (MQ=255)
                                                                                                                gtgtGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCATGGATgat                                  >  6:245731/1‑147 (MQ=255)
                                                                                                                 tgtgCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATgcgc                                 <  4:372772/149‑1 (MQ=255)
                                                                                                                 tgtgCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATgcgc                                 <  1:54069/149‑1 (MQ=255)
                                                                                                                      tGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCCTTCCCATGAGGAAGAAGACCATGACGAGCATGACCACGAGCACCGCACAGGTCGCCTGGATGCGCATCtg                            >  8:63558/1‑149 (MQ=255)
                                                                                                                        ggAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCtgtg                          >  1:301293/1‑149 (MQ=255)
                                                                                                                         gAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACATCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTgg                         >  4:285642/1‑149 (MQ=255)
                                                                                                                            gCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTa                      >  3:253329/1‑149 (MQ=255)
                                                                                                                            gCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTa                      >  8:80228/1‑149 (MQ=255)
                                                                                                                            gCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTa                      <  1:252324/149‑1 (MQ=255)
                                                                                                                              cGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGa                                      <  7:155894/131‑1 (MQ=255)
                                                                                                                              cGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGa                                      >  8:155894/1‑131 (MQ=255)
                                                                                                                               gTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCCCCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTAtcg                   >  5:166350/1‑149 (MQ=255)
                                                                                                                                   cAAGGCCACAGTGGCTGTCCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCgg               >  2:195026/1‑149 (MQ=255)
                                                                                                                                    aaGGCCACAGTGGATGTGCCGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAg                                                                                        >  1:251114/1‑75 (MQ=255)
                                                                                                                                    aaGGCCACAGTGGATGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAg                                                                                        >  1:245695/1‑75 (MQ=255)
                                                                                                                                    aaGGCCACAGTGGATGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAg                                                                                        <  2:245695/75‑1 (MQ=255)
                                                                                                                                    aaGGCCACAGGGGATGTGCAGTCGCTGCCGGGCATGAAGCGGCGCCACTTTGGCGAGGACAGCCATTCCCATGAg                                                                                        <  2:251114/75‑1 (MQ=255)
                                                                                                                                     aGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTc             <  4:122825/149‑1 (MQ=255)
                                                                                                                                         cacaGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATg         >  5:380901/1‑149 (MQ=255)
                                                                                                                                         cacaGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATg         >  7:360795/1‑149 (MQ=255)
                                                                                                                                          acaGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTc             <  2:60223/144‑1 (MQ=255)
                                                                                                                                          acaGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTc             >  1:60223/1‑144 (MQ=255)
                                                                                                                                            aGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATgcgc      >  3:134312/1‑149 (MQ=255)
                                                                                                                                              tGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGTCTGGATGCGCATCTGTGGCTATCGTCGGTCAATgcgc      >  3:360439/1‑147 (MQ=255)
                                                                                                                                              tGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATGCGCGa    <  4:264430/149‑1 (MQ=255)
                                                                                                                                              tGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATGCGCGa    <  7:80228/149‑1 (MQ=255)
                                                                                                                                                gCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATGCGCGAGt  >  2:231869/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCGCCCGGCGCATCCCCCCACCACTATGCGCTGCGCCCCTCCGACGTGCGGCGGGTAGGCGATGCCGACCTGTTGTACTGGATCGGCCCGGACATGGAAAACTTCCTTGCGCGTGTGCTGGGAAGCCGTAGCAAGGCCACAGTGGCTGTGCAGTCGCTGCCGGGCATGAAGCTGCGCCACTTTGGCGAGGACAGCCATTCCCATGAGGAAGAAGACCATGACGATCATGACCACGATCACCGCCCAGGTAGCCTGGATGCGCATCTGTGGCTATCGTCGGTCAATGCGCGAGT  >  NC_002947/125323‑125615

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: ‑