Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,669,734 G→C 8.3% V22V (GTG→GTC) ‡ nuoN → NADH‑quinone oxidoreductase subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,669,7340GC 8.3% 175.7 / 6.1 60V22V (GTG→GTC) ‡nuoNNADH‑quinone oxidoreductase subunit N
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base C (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAA  >  NC_002947/4669588‑4669880
                                                                                                                                                   |                                                                                                                                                   
ccTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggt                                                                                                                                                   <  6:61994/149‑1 (MQ=255)
 ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                                  <  6:110768/149‑1 (MQ=255)
 ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                                  <  5:351809/149‑1 (MQ=255)
 ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                                  <  1:336027/149‑1 (MQ=255)
 ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                                  <  1:314244/149‑1 (MQ=255)
    tGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtggt                                                                                                                                                >  1:294750/1‑148 (MQ=255)
      ccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGAt                                                                                                                                             >  5:332193/1‑149 (MQ=255)
      ccgcaACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCATCTCTCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtggt                                                                                                                                                <  8:338343/146‑1 (MQ=255)
        gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATgg                                                                                                                                           <  5:286576/149‑2 (MQ=255)
           accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCGCCAGCATCACCACCACCGTGGTGAtggtgg                                                                                                                                        <  1:262124/149‑1 (MQ=39)
                          cagcagGGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACCCTTCATCGCCTTTGGCGCCGATGCTGATCACCACCATCACCATGGTGGTGGTGATGCTGGCGATCGCCTGGAAg                                                                                                                          <  1:302102/149‑1 (MQ=255)
                               gTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAAc                                                                                                                    <  3:86387/149‑1 (MQ=255)
                               gTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAAc                                                                                                                    <  3:341943/149‑1 (MQ=255)
                                  gCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCAc                                                                                                                 >  1:72340/1‑149 (MQ=255)
                                    tCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTc                                                                                                               >  2:239703/1‑149 (MQ=255)
                                      ggATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTTGAAGCGCAACCACTc                                                                                                               <  6:332193/147‑1 (MQ=255)
                                      ggATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                              >  8:294105/1‑148 (MQ=255)
                                      ggATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                              <  8:250181/148‑1 (MQ=255)
                                                tCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTc                                                                                                     <  3:283094/147‑1 (MQ=255)
                                                  aCTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgc                                                                                                 >  6:93781/1‑149 (MQ=255)
                                                    tCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgctg                                                                                               >  7:71398/1‑149 (MQ=255)
                                                           gCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCg                                                                                        >  6:247994/1‑149 (MQ=255)
                                                               cGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTg                                                                                     <  7:294105/147‑1 (MQ=255)
                                                                      gTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAAc                                                                            >  2:329767/1‑149 (MQ=255)
                                                                          gagCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                              >  5:321344/1‑112 (MQ=255)
                                                                          gagCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                              <  6:321344/112‑1 (MQ=255)
                                                                          gagCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTgg                                                                        <  2:72340/149‑1 (MQ=255)
                                                                           agCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGc                                                                       >  2:296469/1‑149 (MQ=255)
                                                                             cgcTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGccc                                                                     <  8:71398/149‑1 (MQ=255)
                                                                               cTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctg                                                                   <  1:296469/149‑1 (MQ=255)
                                                                               cTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATACTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctg                                                                   >  1:262489/1‑149 (MQ=255)
                                                                                tATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgc                                                                  >  7:119507/1‑149 (MQ=255)
                                                                                    gAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGt                                                               <  1:329767/148‑1 (MQ=255)
                                                                                        tCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTc                                                          >  1:257457/1‑149 (MQ=255)
                                                                                         caccacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTcc                                                         >  6:229148/1‑149 (MQ=255)
                                                                                         caccacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTc                                                          <  5:275767/148‑1 (MQ=255)
                                                                                          accacTCAATACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCt                                                        >  3:287949/1‑149 (MQ=255)
                                                                                            cacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGc                                                      >  7:203076/1‑149 (MQ=255)
                                                                                               tCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTcc                                                                                                    >  3:30329/1‑101 (MQ=255)
                                                                                               tCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTcc                                                                                                    <  4:30329/101‑1 (MQ=255)
                                                                                               tCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCgg                                                   >  8:235109/1‑149 (MQ=255)
                                                                                                 aaCACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgcg                                                 >  1:263225/1‑149 (MQ=255)
                                                                                                      ttCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTg                                              >  3:254915/1‑147 (MQ=255)
                                                                                                        cATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAgg                                          >  4:331529/1‑149 (MQ=255)
                                                                                                         aTCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAgg                                          <  8:203076/148‑1 (MQ=255)
                                                                                                               aTTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGc                                                      <  2:139741/131‑1 (MQ=255)
                                                                                                               aTTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGc                                                      >  1:139741/1‑131 (MQ=255)
                                                                                                                    cgcCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTg                               <  1:239703/149‑1 (MQ=255)
                                                                                                                    cgcCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTg                               <  2:263225/149‑1 (MQ=255)
                                                                                                                       cGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCg                            <  3:331529/149‑1 (MQ=255)
                                                                                                                        gATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCgg                           >  5:254790/1‑149 (MQ=255)
                                                                                                                          tGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTc                         <  2:294750/149‑1 (MQ=255)
                                                                                                                            cTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTAGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTca                        <  1:302146/149‑1 (MQ=255)
                                                                                                                              gATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACAt                     >  6:320158/1‑149 (MQ=255)
                                                                                                                                  gccagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCt                 >  6:15198/2‑149 (MQ=255)
                                                                                                                                  accaccATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGc                  <  4:287949/148‑1 (MQ=255)
                                                                                                                                  accaccATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCt                 <  5:76741/149‑1 (MQ=255)
                                                                                                                                   ccagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGc                             <  2:35434/136‑1 (MQ=255)
                                                                                                                                   ccagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGc                             >  1:35434/1‑136 (MQ=255)
                                                                                                                                    caccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACCTCGCTGa               >  4:204310/1‑149 (MQ=255)
                                                                                                                                    caccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGa               >  3:275842/1‑149 (MQ=255)
                                                                                                                                       caTCACCGCGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCa            >  8:106119/1‑149 (MQ=255)
                                                                                                                                           accacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCAt        <  3:69940/149‑1 (MQ=255)
                                                                                                                                             cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCAt        <  4:301273/147‑1 (MQ=255)
                                                                                                                                             cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCAt        >  3:301273/1‑147 (MQ=255)
                                                                                                                                             cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATc       >  4:168518/1‑148 (MQ=255)
                                                                                                                                             cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATc       >  3:278415/1‑148 (MQ=255)
                                                                                                                                                 gtggtggtgATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACaa  <  5:93781/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAA  >  NC_002947/4669588‑4669880

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: