Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,829,347 T→C 6.0% Y81C (TAC→TGC)  PP_3382 ← gluconate 2‑dehydrogenase cytochrome c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,829,3470TC 6.0% 192.4 / 2.1 66Y81C (TAC→TGC) PP_3382gluconate 2‑dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base C (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCCACGCCCTTCATGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACT  >  NC_002947/3829200‑3829491
                                                                                                                                                   |                                                                                                                                                
gcCACGCCCTTCATGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTa                                                                                                                                                 <  6:313541/149‑1 (MQ=255)
   aCGCCCTTCATGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTg                                                                                                                                              <  2:113843/149‑1 (MQ=255)
           cATGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGc                                                                                                                                      >  2:67117/1‑149 (MQ=255)
              gAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCAc                                                                                                                                   >  3:18285/1‑149 (MQ=255)
                       gTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACCGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGccgccg                                                                                                                          >  5:147589/1‑149 (MQ=255)
                          gagCGCCTGCATGTCCTGCTCGCTGACCGGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCgggg                                                                                                                       <  4:65780/149‑1 (MQ=255)
                          gagCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCgggg                                                                                                                       >  3:259840/1‑149 (MQ=255)
                                cTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGt                                                                                                                 >  4:50817/1‑149 (MQ=255)
                                cTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGt                                                                                                                 >  8:310982/1‑149 (MQ=255)
                                cTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGCACTGGCCAATGCCACTGGCCGCCGGGGTGATGt                                                                                                                 >  8:337024/1‑149 (MQ=255)
                                cTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGCACTGGCCAATGCCACTGGCCGCCGGGGTGATGt                                                                                                                 >  4:129153/1‑149 (MQ=255)
                                   cATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtgg                                                                                                              <  4:130078/149‑1 (MQ=255)
                                   cATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtg                                                                                                               <  4:259840/148‑1 (MQ=255)
                                     tGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggtg                                                                                                            <  2:200480/149‑1 (MQ=255)
                                     tGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggtg                                                                                                            <  3:293817/149‑1 (MQ=255)
                                     tGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggtg                                                                                                            <  7:49107/149‑1 (MQ=255)
                                      gTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggtg                                                                                                            <  2:309168/148‑1 (MQ=255)
                                      gTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggt                                                                                                             >  1:309168/1‑147 (MQ=255)
                                       tCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggtg                                                                                                            >  3:242738/1‑147 (MQ=255)
                                         cTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAg                                                                                                         <  7:10642/148‑1 (MQ=255)
                                            cTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTa                                                                                                       <  1:3944/147‑1 (MQ=255)
                                            cTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGCACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTac                                                                                                      <  3:129153/148‑1 (MQ=255)
                                            cTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGCACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTac                                                                                                      <  7:337024/148‑1 (MQ=255)
                                             tCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTAGCCGCCGGGGTGATGTTGGTGGAGTacac                                                                                                    >  5:352235/1‑149 (MQ=255)
                                               gCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTacac                                                                                                    <  6:229234/147‑1 (MQ=255)
                                                 tGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTg                                                                                                >  7:186502/1‑149 (MQ=255)
                                                  gACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggt                                                                                                             >  7:44738/1‑135 (MQ=255)
                                                  gACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggt                                                                                                             <  8:44738/135‑1 (MQ=255)
                                                       gcgcATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATc                                                                                          <  4:242738/149‑1 (MQ=255)
                                                              cgacgGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGttt                                                                                   <  3:50817/149‑1 (MQ=255)
                                                               gacgGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGa                                                                                                                    <  3:168459/115‑1 (MQ=255)
                                                               gacgGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGa                                                                                                                    >  4:168459/1‑115 (MQ=255)
                                                                         gcatggcCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAAc                                                                        >  7:295588/1‑149 (MQ=255)
                                                                           atggcCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAAccc                                                                      <  8:295588/149‑1 (MQ=255)
                                                                                 ggATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTg                                                                                                <  5:206763/117‑1 (MQ=255)
                                                                                 ggATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTg                                                                                                >  6:206763/1‑117 (MQ=255)
                                                                                      gagCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAgg                                                           <  8:102066/149‑1 (MQ=255)
                                                                                      gagCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAgg                                                           <  1:67117/149‑1 (MQ=255)
                                                                                      gagCGTACCGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAgg                                                           <  6:147589/149‑1 (MQ=255)
                                                                                       agCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAggg                                                          >  5:217133/1‑149 (MQ=255)
                                                                                                ccGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATTTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAATGGGTTTgccgcc                                                 >  5:118858/1‑149 (MQ=255)
                                                                                                ccGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTgccgc                                                  <  8:186502/148‑1 (MQ=255)
                                                                                                     cTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGc                                                                                                                                            <  5:22450/53‑1 (MQ=255)
                                                                                                     cTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGc                                                                                                                                            >  6:22450/1‑53 (MQ=255)
                                                                                                      ttGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGc                                           >  4:39588/1‑149 (MQ=255)
                                                                                                       tGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCg                                          >  5:298632/1‑149 (MQ=255)
                                                                                                                 ccTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATc                                                                              >  3:317566/1‑103 (MQ=255)
                                                                                                                 ccTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATc                                                                              <  4:317566/103‑1 (MQ=255)
                                                                                                                    gCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCTAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGGCa                             >  8:315425/1‑149 (MQ=255)
                                                                                                                      gTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACg                           <  3:39588/149‑1 (MQ=255)
                                                                                                                        aCGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAAGCCACg                           >  1:124439/1‑147 (MQ=255)
                                                                                                                         cGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCa                         <  4:18285/148‑1 (MQ=255)
                                                                                                                          ggCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTAGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGt                       <  6:352235/149‑1 (MQ=255)
                                                                                                                            ccTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCaa                                                                                            <  6:160575/78‑1 (MQ=255)
                                                                                                                            ccTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCaa                                                                                            >  5:160575/1‑78 (MQ=255)
                                                                                                                                gTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCgg                 >  1:66564/1‑149 (MQ=255)
                                                                                                                                 tCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGgc                >  4:126193/1‑149 (MQ=255)
                                                                                                                                  cGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTgc                                                     >  7:115114/1‑111 (MQ=255)
                                                                                                                                  cGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTgc                                                     <  8:115114/111‑1 (MQ=255)
                                                                                                                                      gTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCggg           <  2:124439/149‑1 (MQ=255)
                                                                                                                                        cTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGc          >  7:107259/1‑148 (MQ=255)
                                                                                                                                         ttCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTTGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGc          >  1:255665/1‑147 (MQ=255)
                                                                                                                                         ttCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCg                  >  3:325474/1‑139 (MQ=255)
                                                                                                                                         ttCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCg                  <  4:325474/139‑1 (MQ=255)
                                                                                                                                              aGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTAc   <  2:66564/149‑1 (MQ=255)
                                                                                                                                                cTGTACTGGCGTTTGCGATTGGCGGTGGGGGTGTTGTTTGTGGATTACACATTGACTTTCGGCGTTTCCATCGGGATCCCGCGTGCATAGGGTTTGCCGCCTTTGGCGTTATTGCAGTCCACGCTGTCTCCGTCGCGGGCCAGGTACt  <  6:298632/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCCACGCCCTTCATGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACT  >  NC_002947/3829200‑3829491

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: