Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 132,467 G→C 5.3% intergenic (+41/‑111) cc → / → dsbA cytochrome c4/thiol:disulfide interchange protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947132,4670GC 5.3% 243.9 / 3.7 76intergenic (+41/‑111)cc/dsbAcytochrome c4/thiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base C (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

GGCAACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGCCGC  >  NC_002947/132319‑132615
                                                                                                                                                    |                                                                                                                                                    
ggCAACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATGCAGGGCCTGAATTAACTTTCGGTTATTGTTGTAGGTCAATGATGAAAGGGCGGCCg                                                                                                                                                      >  8:328969/1‑149 (MQ=255)
 gCAACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGTAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCgg                                                                                                                                                     >  2:183523/1‑149 (MQ=255)
   aaCCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGGTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGc                                                                                                                                                   >  8:168248/1‑149 (MQ=255)
          ccAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCtt                                                                                                                                            <  4:330216/149‑1 (MQ=255)
           cAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCttt                                                                                                                                           <  1:187500/149‑1 (MQ=255)
           cAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCttt                                                                                                                                           >  1:292513/1‑149 (MQ=255)
           cAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCttt                                                                                                                                           >  6:319338/1‑149 (MQ=255)
            aaCGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttc                                                                                                                                          <  4:188693/149‑2 (MQ=255)
              cGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttcat                                                                                                                                        <  2:197934/149‑4 (MQ=255)
                   gCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTg                                                                                                                                                 >  1:159552/1‑135 (MQ=255)
                   gCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTg                                                                                                                                                 <  2:159552/135‑1 (MQ=255)
                   gCGATGCCATGACCATGCGCACCATCGCTGGAAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc                                                                                                                                   <  1:183523/149‑1 (MQ=255)
                       tGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCa                                                                                                                               <  2:292513/149‑1 (MQ=255)
                         ccATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGc                                                                                                                             <  3:262169/149‑1 (MQ=255)
                          cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGGTGCCCTTTTTTCCGTCCGCAGcc                                                                                                                            >  3:298674/1‑149 (MQ=255)
                          cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGcc                                                                                                                            >  3:88738/1‑149 (MQ=255)
                           aTGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCg                                                                                                                           >  5:132858/1‑149 (MQ=255)
                                cATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTaca                                                                                                                      >  7:2827/1‑149 (MQ=255)
                                cATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTaca                                                                                                                      <  8:24552/149‑1 (MQ=255)
                                cATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTa                                                                                                                        <  8:2827/147‑1 (MQ=255)
                                    cgcACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTa                                                                                                                   <  3:132360/148‑1 (MQ=255)
                                    cgcACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTa                                                                                                                   >  5:236829/1‑147 (MQ=255)
                                     gcACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACa                                                                                                                 <  8:123839/149‑1 (MQ=255)
                                     gcACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGATGTAGGTCATTGATGAAAGGGCGGCCGGGCTGCCCTGTTTTCCGTCCGCAGCCGTTACCCTa                                                                                                                   >  8:91689/1‑147 (MQ=255)
                                         cATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAAc                                                                                                             >  5:231369/1‑149 (MQ=255)
                                          aTCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACt                                                                                                            >  8:233673/1‑149 (MQ=255)
                                              cTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACa                                                                                                                 >  1:341742/1‑140 (MQ=255)
                                              cTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACa                                                                                                                 <  2:341742/140‑1 (MQ=255)
                                                 gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGAGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGGTACACTACAGAACTCTAGcac                                                                                                     >  8:27329/1‑147 (MQ=255)
                                                   aaGCTGAGCAACCATGACATCGAAGCGTTGGCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTc                                                                                                  <  6:236829/149‑1 (MQ=255)
                                                    aGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTc                                                                                                  <  7:77117/149‑1 (MQ=255)
                                                         aGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGc                                                                                             <  2:76352/149‑1 (MQ=255)
                                                            ggccATGGCATCGAAGCGTTGGCCAGCTACATCCAGGGGCTGGATTAACTTTCGGGTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGGTGGGGTGTTTTGCGGCCGCAGCGGGTGCCGTACAGAACTCGGGCCGGGCCCGGCCTg                                                                                          >  4:296919/3‑149 (MQ=255)
                                                               cATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGTCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtct                                                                                       <  4:88738/149‑1 (MQ=255)
                                                               cATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtct                                                                                       >  8:339183/1‑149 (MQ=255)
                                                               cATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtct                                                                                       <  6:132858/149‑1 (MQ=255)
                                                                 tGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtct                                                                                       >  1:286233/1‑147 (MQ=255)
                                                                 tGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCa                                                                                     <  5:11069/149‑1 (MQ=255)
                                                                    cATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGt                                                                                   <  3:239718/148‑1 (MQ=255)
                                                                                   cAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTgg                                                                   >  8:261282/1‑149 (MQ=255)
                                                                                    aGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTGGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGc                                                                  <  5:319338/149‑1 (MQ=255)
                                                                                    aGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCAGCATCAGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTATCAGTGCAGGTCGCGTCGTGGc                                                                  >  7:33940/1‑149 (MQ=255)
                                                                                      cTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGa                                                                >  2:169785/1‑149 (MQ=255)
                                                                                         cATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAcc                                                              >  1:95906/1‑148 (MQ=255)
                                                                                         cATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAc                                                               <  4:118251/147‑1 (MQ=255)
                                                                                          aTCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCa                                                            >  7:70267/1‑149 (MQ=255)
                                                                                           tCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCCACGACGCGAcctgc                                                           <  6:231369/149‑4 (MQ=255)
                                                                                             cAGGTCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGa                                                         <  4:361131/149‑1 (MQ=255)
                                                                                             cAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAACGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCCACGACGCGAcctgcac                                                         <  7:328969/149‑6 (MQ=255)
                                                                                              aGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCATTGATGAATGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGAc                                                        <  4:93269/149‑1 (MQ=255)
                                                                                              aGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGAc                                                        <  7:339183/149‑1 (MQ=255)
                                                                                               gggCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCTGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACt                                                       <  7:91689/149‑1 (MQ=255)
                                                                                               gggCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACt                                                       >  4:185227/1‑149 (MQ=255)
                                                                                                ggCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTg                                                      <  7:233673/149‑1 (MQ=255)
                                                                                                 gCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCTGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGt                                                                             <  8:34800/125‑1 (MQ=255)
                                                                                                 gCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGt                                                                             >  7:34800/1‑125 (MQ=255)
                                                                                                 gCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGc                                                     >  7:337541/1‑149 (MQ=255)
                                                                                                       aTTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg                                                 >  7:159101/1‑147 (MQ=255)
                                                                                                             tttCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGCGCCCTTCCCGGCCGGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAg                                         >  4:83235/1‑149 (MQ=255)
                                                                                                                 gttttaTGTTGTAGGTCAATGATGTAAGGGCGGGCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTTCACTACAGATCTCGAGCCCTTCGCGTCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATg                                     <  4:23842/144‑1 (MQ=255)
                                                                                                                  gTTAATGTTGTAGGTCAATGATGACAGCGCGGCCGGGCTGCCCTTTTTTCCCTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                    <  4:298674/149‑1 (MQ=255)
                                                                                                                  gTTAATGTTGTAGGTCAATGATGAAAGGGCCGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGCGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                    >  3:101842/1‑149 (MQ=255)
                                                                                                                       tgttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTaa                                >  3:59250/1‑148 (MQ=255)
                                                                                                                            tAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGAt                          >  4:12768/1‑149 (MQ=255)
                                                                                                                              ggTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATtc                        >  7:214239/1‑149 (MQ=255)
                                                                                                                               gTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGAt                          >  6:133706/1‑146 (MQ=255)
                                                                                                                               gTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGAt                          <  5:133706/146‑1 (MQ=255)
                                                                                                                               gTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATtct                       >  4:125285/1‑149 (MQ=255)
                                                                                                                                  aaTGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCa                     <  8:337541/148‑1 (MQ=255)
                                                                                                                                  aaTGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAg                    <  8:70267/149‑1 (MQ=255)
                                                                                                                                      atgaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACGGCTCAGGAGTAAAGCATGCGTAAACTGCTTCTCAgcgct                >  5:157047/1‑149 (MQ=255)
                                                                                                                                      aaaaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgct                >  4:28125/4‑149 (MQ=255)
                                                                                                                                      aaaaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgct                >  8:306241/4‑149 (MQ=255)
                                                                                                                                       tgaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgct                <  8:159101/148‑1 (MQ=255)
                                                                                                                                        gaAAGGGCGGCCGGGCTGCCGTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctgc              >  7:275845/1‑149 (MQ=255)
                                                                                                                                        gaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctgc              >  7:74974/1‑149 (MQ=255)
                                                                                                                                              gcggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGa                                                         <  5:354292/100‑1 (MQ=255)
                                                                                                                                              gcggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGa                                                         >  6:354292/1‑100 (MQ=255)
                                                                                                                                                    gggCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgccgc  <  2:286233/149‑1 (MQ=255)
                                                                                                                                                    gggCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGGCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgccgc  <  7:27329/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCAACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGCCGC  >  NC_002947/132319‑132615

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: