Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,781,474 A→T 6.6% L1947Q (CTG→CAG)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,781,4740AT 6.6% 231.2 / 4.8 76L1947Q (CTG→CAG) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base T (3/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

TCCGGGGCACTGCCAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTC  >  NC_002947/4781327‑4781619
                                                                                                                                                   |                                                                                                                                                 
tCCGGGGCACTGCCAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAg                                                                                                                                                  >  5:337526/1‑149 (MQ=255)
       cACTGCCAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGc                                                                                                                                           <  6:81681/149‑1 (MQ=255)
            ccAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCt                                                                                                                                       >  1:206580/1‑148 (MQ=255)
            ccAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCt                                                                                                                                       >  2:224166/1‑148 (MQ=255)
            ccAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTc                                                                                                                                      <  8:25559/149‑1 (MQ=255)
             cAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCTGGGCCTGCACCtgt                                                                                                                                     <  2:373011/149‑3 (MQ=255)
                gAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCa                                                                                                                                          <  8:251258/141‑1 (MQ=255)
                gAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCa                                                                                                                                          >  7:251258/1‑141 (MQ=255)
                     gcggcggcCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTa                                                                                                                             <  6:337526/149‑1 (MQ=255)
                     gcggcggcCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTa                                                                                                                             <  7:160592/149‑1 (MQ=255)
                               caGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAgg                                                                                                                                                 <  7:207281/119‑1 (MQ=255)
                               caGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAgg                                                                                                                                                 >  8:207281/1‑119 (MQ=255)
                                 gCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGt                                                                                                                 >  3:348382/1‑149 (MQ=255)
                                     ggAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGc                                                                                                             >  7:552/1‑149 (MQ=255)
                                         gTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCggg                                                                                                          >  2:273984/1‑148 (MQ=255)
                                          tCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTa                                                                                                        <  4:348382/149‑1 (MQ=255)
                                            cgcgTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCTGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGt                                                                                                      >  6:54023/1‑149 (MQ=255)
                                            cgcgTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGt                                                                                                      >  2:285436/1‑149 (MQ=255)
                                              cgTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTg                                                                                                    <  2:2361/149‑1 (MQ=255)
                                                tCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCCGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGgt                                                                                                  <  6:191834/149‑1 (MQ=255)
                                                tCCAGCTGCTCGCTCAAGGCCAGCAACAAACCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCGGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGgt                                                                                                  <  3:25315/149‑1 (MQ=255)
                                                  cAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGgtgt                                                                                                <  6:377935/149‑1 (MQ=255)
                                                   aGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGtgtg                                                                                               >  5:354752/1‑149 (MQ=255)
                                                   aGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGACGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGtgtg                                                                                               >  7:187219/1‑149 (MQ=255)
                                                    gctgctCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGACGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTgtgc                                                                                              <  8:187219/149‑1 (MQ=255)
                                                           gctcAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCgg                                                                                                                                >  4:18298/1‑108 (MQ=255)
                                                           gctcAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCgg                                                                                                                                <  3:18298/108‑1 (MQ=255)
                                                                 ggCCAGCAACAAAACCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTcc                                                                                 >  4:235780/1‑149 (MQ=255)
                                                                   ccAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCa                                                                               <  4:17105/149‑1 (MQ=255)
                                                                    cAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAg                                                                              >  1:251781/1‑149 (MQ=255)
                                                                    cAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCACCACCACACGCGGGTACCCCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAg                                                                              >  1:333925/1‑149 (MQ=255)
                                                                           aaaGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGAc                                                                       >  6:192516/1‑149 (MQ=255)
                                                                                ccGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCa                                                                  >  5:356584/1‑149 (MQ=255)
                                                                                  ggAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGTGTGCGCGAACTCCCAGCCTTGACCTTCATg                                                                >  1:33517/1‑149 (MQ=255)
                                                                                  ggAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATg                                                                <  6:354752/149‑1 (MQ=255)
                                                                                    aaGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGt                                                              >  8:91357/1‑149 (MQ=255)
                                                                                     aGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTc                                                             <  3:175901/149‑1 (MQ=255)
                                                                                         agccagcAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGaa                                                          <  4:6331/148‑1 (MQ=255)
                                                                                          gccagcAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGc                                                        <  5:192516/149‑1 (MQ=255)
                                                                                                agcaCGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGgc                                                  <  6:356584/149‑1 (MQ=255)
                                                                                                agcaCGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGgc                                                  <  2:251781/149‑1 (MQ=255)
                                                                                                agcaCGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGgc                                                  <  1:273984/149‑1 (MQ=255)
                                                                                                     gCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGACGAAGCCTGGGCGCCaa                                             >  8:263632/1‑149 (MQ=255)
                                                                                                       tGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACg                                           >  5:357188/1‑149 (MQ=255)
                                                                                                         aTGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGcc                                                                            >  4:297626/1‑114 (MQ=255)
                                                                                                         aTGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGcc                                                                            <  3:297626/114‑1 (MQ=255)
                                                                                                          tGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTg                                        <  7:29415/149‑1 (MQ=255)
                                                                                                            gTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGa                                      >  6:15947/1‑149 (MQ=255)
                                                                                                                 gcTGTCGAACAGGCGTTGGTCATACGACAGGTGCTGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCa                                 <  5:54023/149‑1 (MQ=255)
                                                                                                                 gcTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCa                                 >  1:386232/1‑149 (MQ=255)
                                                                                                                  cTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCAc                                <  2:386232/149‑1 (MQ=255)
                                                                                                                   tGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACg                               >  5:159953/1‑149 (MQ=255)
                                                                                                                   tGTCGAACAGGCGTTGGTCATACGACAGCTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACg                               >  1:229184/1‑149 (MQ=255)
                                                                                                                     tCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACggg                             <  2:206580/149‑1 (MQ=255)
                                                                                                                       gAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGgcg                           >  3:265193/1‑149 (MQ=255)
                                                                                                                         aCAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGGGAGCGGGGAGTTGGTGTGGGTGCGCGAACTCCCCGCCTTGACCTTCAAGGTCGAAGCCTGGGCGCCAACGCGGGCCGCCACGGGcgcg                         >  6:63484/1‑149 (MQ=255)
                                                                                                                           aGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGTGCGCCAACGCTGGACTCCACGGGCGCGtt                       >  6:44504/1‑149 (MQ=255)
                                                                                                                           aGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTg                                        <  1:316718/132‑1 (MQ=255)
                                                                                                                           aGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTg                                        >  2:316718/1‑132 (MQ=255)
                                                                                                                            ggCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACACGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTtc                      <  6:6339/149‑1 (MQ=255)
                                                                                                                             gCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTg                                        >  5:61438/1‑130 (MQ=255)
                                                                                                                             gCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTg                                        <  6:61438/130‑1 (MQ=255)
                                                                                                                             gCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTtct                     >  8:252894/1‑149 (MQ=255)
                                                                                                                                  ggTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaa                 <  5:15947/148‑1 (MQ=255)
                                                                                                                                  ggTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaa                 <  8:552/148‑1 (MQ=255)
                                                                                                                                   gTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaaca               >  1:121742/1‑149 (MQ=255)
                                                                                                                                   gTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaaca               <  7:91357/149‑1 (MQ=255)
                                                                                                                                    tCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacac              <  1:224166/149‑1 (MQ=255)
                                                                                                                                      aTACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCTTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacacga            <  3:235780/149‑1 (MQ=255)
                                                                                                                                        aCGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacacgaa           >  2:219524/1‑148 (MQ=255)
                                                                                                                                        aCGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacacgaa           >  4:108770/1‑148 (MQ=255)
                                                                                                                                            gaGGTGCTGGGCCTGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGc      >  1:113957/2‑149 (MQ=255)
                                                                                                                                               gTGCTGGGCCTGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCTCGTTCTCGAACACGAACAGGCTGt   >  6:32924/1‑149 (MQ=255)
                                                                                                                                               gTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGt   >  6:358096/1‑149 (MQ=255)
                                                                                                                                               gTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGt   >  2:12117/1‑149 (MQ=255)
                                                                                                                                                tGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTc  <  3:108770/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCCGGGGCACTGCCAGGAAGGGCGGCGGCCACAGCATGGAAGTCCGCGTCCAGCTGCTCGCTCAAGGCCAGCAACAAAGCCCGGAAGTCAGCCAGCAGCACGCTGATGGTCGCGCTGTCGAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTC  >  NC_002947/4781327‑4781619

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: