Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,640 C→A 5.4% R28S (CGC→AGC)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,6400CA 5.4% 271.5 / 3.0 84R28S (CGC→AGC) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (43/36);  new base A (2/3);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCA  >  NC_002947/980494‑980787
                                                                                                                                                   |                                                                                                                                                     
ttCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATcgc                                                                                                                                                     <  1:139389/149‑1 (MQ=255)
   aacaacCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAg                                                                                                                                                  <  5:8709/149‑1 (MQ=255)
      aacCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTg                                                                                                                                               <  4:218306/149‑1 (MQ=255)
       acCCGAACGTGCGCTGTGTGTGTGGCTGCTGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTATGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                              <  7:75735/149‑1 (MQ=255)
        ccccAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGc                                                                                                                                             >  6:309144/1‑149 (MQ=255)
         cccAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGcc                                                                                                                                            >  2:253855/1‑149 (MQ=255)
         cccAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCGGGCCAGCTCGCGATAGCGAgtggc                                                                                                                                              >  8:159984/1‑144 (MQ=255)
           cAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTg                                                                                                                                               <  1:60207/144‑1 (MQ=255)
           cAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTg                                                                                                                                               >  2:60207/1‑144 (MQ=255)
            atcGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAGTGGCCagc                                                                                                                                         <  8:261597/147‑1 (MQ=255)
            aaCGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCcgc                                                                                                                                         >  4:92528/1‑149 (MQ=255)
            aaCGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCcgc                                                                                                                                         <  8:28522/149‑1 (MQ=255)
             aCGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCcgcg                                                                                                                                        <  4:153485/149‑1 (MQ=255)
                tGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAGTGGCCagcttgt                                                                                                                                     <  3:105734/149‑5 (MQ=255)
                    cgGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCa                                                                                                                                 <  3:92528/149‑1 (MQ=255)
                        gAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                              >  7:25992/1‑132 (MQ=255)
                        gAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                              <  8:25992/132‑1 (MQ=255)
                                cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCtt                                                                                                                     >  5:160568/1‑149 (MQ=255)
                                    ggcgAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGc‑‑c                                                                                                               <  6:17684/149‑1 (MQ=255)
                                       gAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAGTGGCCagcttg                                                                                                                                      <  8:196119/125‑4 (MQ=38)
                                       gAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAGTGGCCagcttg                                                                                                                                      >  7:196119/1‑122 (MQ=38)
                                       gAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGAc                                                                                                            >  8:194130/1‑149 (MQ=255)
                                              tCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTcc                                                                                                      >  2:114303/1‑148 (MQ=255)
                                                aaCGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGa                                                                                                    <  5:309144/148‑1 (MQ=255)
                                                 aCGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                              <  3:295680/107‑1 (MQ=255)
                                                 aCGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTgg                                                                                                                                              >  4:295680/1‑107 (MQ=255)
                                                   gTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGt                                                                                                <  6:288146/149‑1 (MQ=255)
                                                    tttGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTTGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGt                                                                                                >  2:250747/1‑148 (MQ=255)
                                                    tttGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGt                                                                                                >  2:85743/1‑148 (MQ=255)
                                                      tGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGa                                                                                              <  7:194130/148‑1 (MQ=255)
                                                              gctgTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTa                                                                                      <  3:128731/148‑1 (MQ=255)
                                                               ctgTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTAtt                                                                                    <  6:160568/149‑1 (MQ=255)
                                                                      tACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGTTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTg                                                                              >  1:99459/1‑148 (MQ=255)
                                                                      tACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTgg                                                                             <  3:258766/149‑1 (MQ=255)
                                                                         tCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                               <  6:155940/144‑1 (MQ=255)
                                                                         tCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                               >  5:155940/1‑144 (MQ=255)
                                                                          ccTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGaaa                                                                          >  1:10739/1‑148 (MQ=255)
                                                                           cTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGaaaa                                                                         >  7:109224/1‑148 (MQ=255)
                                                                               tCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCcg                                                                    >  1:15501/1‑149 (MQ=255)
                                                                                cATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGAc                                                                                                            >  6:265719/1‑108 (MQ=255)
                                                                                cATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGAc                                                                                                            <  5:265719/108‑1 (MQ=255)
                                                                                  ttACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCgg                                                                 >  5:102851/1‑149 (MQ=255)
                                                                                    aCGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGcc                                                               >  6:196302/1‑149 (MQ=255)
                                                                                      gCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                               <  7:63945/131‑1 (MQ=255)
                                                                                      gCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                               >  8:63945/1‑131 (MQ=255)
                                                                                            tttGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacga                                                       >  3:141949/1‑149 (MQ=255)
                                                                                            tttGACCTCCAGACAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacga                                                       >  8:131582/1‑149 (MQ=255)
                                                                                                aCCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGa                                                                                                    <  6:86251/100‑1 (MQ=255)
                                                                                                aCCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGa                                                                                                    >  5:86251/1‑100 (MQ=255)
                                                                                                aCCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                    <  4:141949/148‑1 (MQ=255)
                                                                                                 ccTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCtt                                                  <  1:114303/149‑1 (MQ=255)
                                                                                                 ccTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCt                                                   >  3:6365/1‑148 (MQ=255)
                                                                                                   tCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCACCGACGCTTAc                                                >  8:286700/1‑149 (MQ=255)
                                                                                                   tCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                >  4:175940/1‑149 (MQ=255)
                                                                                                   tCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCCAAGACCATCCTGATCGCAAACC‑‑CGACGCCTCCCAGCGTGAGCAGCGAATTGCGCTCGAAAAATCCGCGGCCCTAAACCACGCATAc                                                >  5:212927/1‑149 (MQ=255)
                                                                                                   tCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTAGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                >  1:270167/1‑149 (MQ=255)
                                                                                                     cAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                <  5:196302/147‑1 (MQ=255)
                                                                                                       gCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                <  8:220873/145‑1 (MQ=255)
                                                                                                       gCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                >  3:242690/1‑145 (MQ=255)
                                                                                                       gCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                >  7:220873/1‑145 (MQ=255)
                                                                                                       gCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                <  4:242690/145‑1 (MQ=255)
                                                                                                       gCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGa                                            >  4:267163/1‑149 (MQ=255)
                                                                                                         cAAGATTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCAGCGAAGTCCATCCTGATCGCTTTGA‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACg                                          >  5:223505/1‑149 (MQ=255)
                                                                                                           aGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCt                                        <  4:102652/149‑1 (MQ=255)
                                                                                                             cTTCCGTCTGGATCTCGACAAGCTGGCCACTTGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgc                                      <  1:250747/149‑1 (MQ=255)
                                                                                                             cTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgc                                      >  5:215392/1‑149 (MQ=255)
                                                                                                              ttCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                >  7:13616/1‑138 (MQ=255)
                                                                                                              ttCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                <  8:13616/138‑1 (MQ=255)
                                                                                                               tCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgcg                                     >  4:237656/1‑148 (MQ=255)
                                                                                                                 cGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                    >  7:153581/1‑69 (MQ=255)
                                                                                                                 cGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                    <  8:153581/69‑1 (MQ=255)
                                                                                                                  gTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGt                                 <  1:85743/149‑1 (MQ=255)
                                                                                                                       gATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGcgcc                            >  7:292716/1‑149 (MQ=255)
                                                                                                                           tcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCGGGGCCGCCTCCGCGCGTGCGCAGCGTTTTGCGGTGGAAAAAACCGCGGCCCTCACCGCCGGTTACCACACGGTGCGCAGTGcgccgc                          >  1:93263/1‑149 (MQ=255)
                                                                                                                             gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcg                      >  5:9303/1‑149 (MQ=255)
                                                                                                                                  gCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                 <  3:175940/149‑1 (MQ=255)
                                                                                                                                  gCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgc                 <  2:168029/149‑1 (MQ=255)
                                                                                                                                   cTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcgcg                <  4:6365/149‑1 (MQ=255)
                                                                                                                                        cACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgc                                      <  6:85999/122‑1 (MQ=255)
                                                                                                                                        cACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgc                                      >  5:85999/1‑122 (MQ=255)
                                                                                                                                        cACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAcc           >  6:154782/1‑149 (MQ=255)
                                                                                                                                              gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg      >  8:127618/1‑149 (MQ=255)
                                                                                                                                              gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgct       <  6:215392/148‑1 (MQ=255)
                                                                                                                                              gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgct       <  8:292716/148‑1 (MQ=255)
                                                                                                                                                  tCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCa  >  4:38981/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCA  >  NC_002947/980494‑980787

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: