Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,631:1 +G 6.2% coding (73/522 nt) hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,6311.G 6.2% 327.3 / 15.4 81coding (73/522 nt)hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base . (41/35);  new base G (2/3);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACC  >  NC_002947/980494‑980778
                                                                                                                                          |                                                                                                                                                     
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                                       >  3:242690/1‑145 (MQ=255)
                                                                                                       gCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                       <  8:220873/145‑1 (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)
                                                                                                                                          |                                                                                                                                                     
TTCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGC‑‑CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACC  >  NC_002947/980494‑980778

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: