Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 922,773 G→C 6.7% R85R (CGC→CGG)  PP_0803 ← protein secretion ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947922,7730GC 6.7% 210.5 / 3.5 71R85R (CGC→CGG) PP_0803protein secretion ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base C (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAGGCCGCCCAGGCCA  >  NC_002947/922626‑922918
                                                                                                                                                   |                                                                                                                                                     
gcgcAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCcgc                                                                                                                                                      >  7:150918/1‑149 (MQ=255)
    aGGCTCTCCGGGAAGGGCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCGCGGCCAGGTGCTGGCCCAGCTc                                                                                                                                                  <  7:122847/149‑1 (MQ=255)
       ctctCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCGCGGCCAGGTGCTGGCCCAGCTCGAc                                                                                                                                               <  6:208345/149‑1 (MQ=255)
          tcCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCACTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGATCGACCAg                                                                                                                                            >  7:27523/1‑149 (MQ=255)
            cGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAg                                                                                                                                            <  2:115848/147‑1 (MQ=255)
                     cAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCgg                                                                                                                                >  7:177670/1‑149 (MQ=255)
                      aGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCcgc                                                                                                                                                      >  6:88366/1‑127 (MQ=255)
                        tGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCa                                                                                                                             >  8:46629/1‑149 (MQ=255)
                        tGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCa                                                                                                                             >  8:228807/1‑149 (MQ=255)
                           cTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGc                                                                                                                          >  6:63966/1‑149 (MQ=255)
                           cTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGc                                                                                                                          >  6:258135/1‑149 (MQ=255)
                           cTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGc                                                                                                                          <  3:168767/149‑1 (MQ=255)
                              gCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTc                                                                                                                        <  8:150918/148‑1 (MQ=255)
                                 ggTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCt                                                                                                                    >  4:50671/1‑149 (MQ=255)
                                  gTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCCCGCTCATCTg                                                                                                                   >  1:156093/1‑149 (MQ=255)
                                   tCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTgc                                                                                                                  <  7:46629/149‑1 (MQ=255)
                                     aCCTCTTCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGTTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTgc‑‑‑gc                                                                                                             <  8:276229/149‑1 (MQ=255)
                                     aCCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTgc‑‑‑gc                                                                                                             <  6:153632/149‑1 (MQ=255)
                                      ccTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTgc‑‑‑gc                                                                                                             >  5:256831/1‑148 (MQ=255)
                                       ctttgcctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑Gcga                                                                                                           >  3:99436/4‑149 (MQ=255)
                                       ctctGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑Gcga                                                                                                           <  7:30235/149‑1 (MQ=255)
                                       ctctGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCGCGGCCAGGTGCTGGCCCAGCTCGAcccg                                                                                                                                            <  3:175294/120‑3 (MQ=37)
                                       ctctGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCGCGGCCAGGTGCTGGCCCAGCTCGAcccg                                                                                                                                            >  4:175294/1‑118 (MQ=37)
                                        tctGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑Gcgac                                                                                                          >  3:72749/1‑149 (MQ=255)
                                          tgcctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑Gcgacga                                                                                                        >  1:2505/1‑149 (MQ=255)
                                            cctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGAt                                                                                                       <  5:9758/148‑1 (MQ=255)
                                               gcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCg                                                                                                   >  1:139189/1‑149 (MQ=255)
                                                cAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGTGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGc                                                                                                  >  5:94377/1‑149 (MQ=255)
                                                 aGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGcc                                                                                                 >  8:136175/1‑149 (MQ=255)
                                                  ggCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGccc                                                                                                <  4:99436/149‑1 (MQ=255)
                                                   gCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCt                                                                                               <  8:177670/149‑1 (MQ=255)
                                                   gCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCt                                                                                               >  6:106515/1‑149 (MQ=255)
                                                   gCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCt                                                                                               <  7:228807/149‑1 (MQ=255)
                                                       gCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCCGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGCTACCGCCcgcgca                                                                                           >  7:224500/1‑147 (MQ=255)
                                                        ccTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAAc                                                                                          <  2:139189/149‑1 (MQ=255)
                                                             ctggcCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTg                                                                                     <  6:256831/149‑1 (MQ=255)
                                                                 cccTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAAt                                                                                  >  4:88296/1‑148 (MQ=255)
                                                                      gCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCt                                                                             <  6:94377/148‑1 (MQ=255)
                                                                         gcgcGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGctc                                                                         <  3:288611/149‑1 (MQ=255)
                                                                           gcgGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCCGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑CCGACGAAACCGCCCTCCAACGACTGAATGACCTGCtctc                                                                       >  4:168355/1‑149 (MQ=255)
                                                                            cgGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTcg                                                                      <  7:204735/149‑1 (MQ=255)
                                                                             ggTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGCCTGAATGACCTGCTCTcgc                                                                     >  3:253166/1‑149 (MQ=255)
                                                                               tACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGa                                                                   <  5:106515/149‑1 (MQ=255)
                                                                                aCTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGaa                                                                  >  5:189238/1‑149 (MQ=255)
                                                                                     gCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTg                                                              <  8:27523/148‑1 (MQ=255)
                                                                                     gCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTg                                                              <  2:2505/148‑1 (MQ=255)
                                                                                      ccTCGCTTTCACCTAGGCTGGAAGCGGTCTTGGTCGGGTCGTGGCGGGCCAGCAGCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGTATGACCTGCTCTCGCGAACTTGGa                                                            <  8:224500/149‑1 (MQ=255)
                                                                                       cTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGCCCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGaa                                                           >  5:196893/1‑149 (MQ=255)
                                                                                       cTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGaa                                                           <  6:189238/149‑1 (MQ=255)
                                                                                       cTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGaa                                                           <  4:72749/149‑1 (MQ=255)
                                                                                       cTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTACTATAAGACGATACCGCCCTCGAACGAAGGAATGACCTGCTATAGCGAACTTg                                                              >  1:268979/1‑149 (MQ=255)
                                                                                          gCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATg                                                         <  3:88296/148‑1 (MQ=255)
                                                                                                 cccACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGcc                                                 >  6:147625/1‑149 (MQ=255)
                                                                                                   cACGCTGGAAGCCGTCTTGGTCCGGTCGACCTCGGAAAGCACCTGGCCGCGCACCACCAGGTC‑ACCTTCCGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAAAGACTGAATGACCAGCTCTCGCGAACTAGGAATGACCTTACCtg                                               >  1:90776/1‑149 (MQ=255)
                                                                                                          gAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGt                                        >  2:294226/1‑149 (MQ=255)
                                                                                                           aaGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTg                                       <  7:136175/149‑1 (MQ=255)
                                                                                                           aaGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTg                                       <  5:147625/149‑1 (MQ=255)
                                                                                                              ccGTGTTGGTCGGGTCGTGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCAGGCTCTTCTGC‑‑‑GCGTCGATACCGCCCTCGAACGACTGAATGACCTGCTGTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGa                                     <  2:275517/148‑1 (MQ=255)
                                                                                                                     ggtcgggtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCg                              <  3:50671/148‑1 (MQ=255)
                                                                                                                      gtcgggtcgAGCTGGGCCAGCACCAGGCCCCGCACGACCAGGTC‑ACCTTAGGCCACGCTCATATGC‑‑‑GCGACGAAACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAAAGACCTTGCCTGTGCCGGTGGATACCTCGGt                            >  7:183853/1‑149 (MQ=255)
                                                                                                                         gggtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCAc                         <  1:294226/149‑1 (MQ=255)
                                                                                                                         ggggcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCAc                         <  3:168355/149‑1 (MQ=255)
                                                                                                                          ggtcgTGCTGGGCCAGCACGTGGCCGCGCTCGACCAGGTC‑ACGTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCAcc                        <  4:253166/149‑1 (MQ=255)
                                                                                                                            tcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTc                      >  5:86453/1‑149 (MQ=255)
                                                                                                                              gAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCTACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGACCGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTc                      >  1:217425/1‑148 (MQ=255)
                                                                                                                                gCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAAc                  >  2:201680/1‑149 (MQ=255)
                                                                                                                                      ccAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGt                            <  1:78442/133‑1 (MQ=255)
                                                                                                                                      ccAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGt                            >  2:78442/1‑133 (MQ=255)
                                                                                                                                       cAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAGgccg           >  4:290623/1‑149 (MQ=255)
                                                                                                                                         gCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAGgccgcc         >  6:266614/1‑149 (MQ=255)
                                                                                                                                                gCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAGGCCGCccaggcca  <  6:86927/149‑1 (MQ=255)
                                                                                                                                                gCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAGGCCGCccaggc    >  5:86927/1‑147 (MQ=255)
                                                                                                                                                gCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAGGCCGCccaggc    >  2:8191/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTC‑ACCTTCGGCCACGCTCATCTGC‑‑‑GCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAGGCCGCCCAGGCCA  >  NC_002947/922626‑922918

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