Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,830,790 C→G 9.7% E134Q (GAA→CAA) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,830,7900CG 9.7% 185.1 / 10.2 62E134Q (GAA→CAA) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base G (2/4);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

CCCTTGGGCAGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCGCATACA  >  NC_002947/4830642‑4830936
                                                                                                                                                     |                                                                                                                                                  
cccTTGGGCAGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTc                                                                                                                                                    <  1:75569/149‑1 (MQ=255)
 ccTTGGGCAGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTcc                                                                                                                                                   <  5:15934/149‑1 (MQ=255)
      ggCAGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGcc                                                                                                                                              <  8:6630/149‑1 (MQ=255)
      ggCAGCGCGGGGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGcc                                                                                                                                              <  8:157776/149‑1 (MQ=255)
         aGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCCTGGAAGCCCTg                                                                                                                                           <  2:117886/149‑1 (MQ=255)
            gcgGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCCTGGAAGCCCTggcg                                                                                                                                        <  7:250211/149‑4 (MQ=255)
                 ggAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTACGCAGACAGGGCCGGGACCAAGCCAACCACAGCCAGCAGTAGGGGGGCGTTGACCGCCGCCAGGGCAACAAGCCCCAGCAggaag                                                                                                                                    >  2:295335/1‑145 (MQ=255)
                            tgtACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGcgcc                                                                                                                        <  2:101863/149‑1 (MQ=255)
                               aCTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGcgccgc                                                                                                                      >  3:228053/1‑148 (MQ=255)
                                      gAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTc                                                                                                               <  6:151264/149‑1 (MQ=255)
                                        aGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTcag                                                                                                             >  7:274547/1‑148 (MQ=255)
                                          gCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTcagcag                                                                                                          <  6:278468/149‑1 (MQ=255)
                                          gCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTcagcag                                                                                                          <  8:274547/149‑1 (MQ=255)
                                          gCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCCTGGAAGCCCTgg                                                                                                                                          <  7:75440/117‑2 (MQ=255)
                                          gCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCCTGGAAGCCCTgg                                                                                                                                          >  8:75440/1‑116 (MQ=255)
                                            gATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTcagcag                                                                                                          >  6:46480/1‑147 (MQ=255)
                                                  aTCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTcgcg                                                                                                  >  6:160016/1‑149 (MQ=255)
                                                   tCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCgg                                                                                                 >  6:199021/1‑149 (MQ=255)
                                                       ccTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCg                                                                                             >  5:115029/1‑149 (MQ=255)
                                                              cACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCCCGGTCAGCAGCAGTCGCGGCTCGGCGTCAt                                                                                      >  5:134517/1‑149 (MQ=255)
                                                                  ggTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCgat                                                                                  <  1:295335/149‑1 (MQ=255)
                                                                    ttCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCgatga                                                                                >  8:6078/1‑149 (MQ=255)
                                                                       aCGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATc                                                                              >  6:308732/1‑148 (MQ=255)
                                                                       aCGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATc                                                                              >  7:151637/1‑148 (MQ=255)
                                                                       aCGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATc                                                                              <  4:228053/148‑1 (MQ=255)
                                                                       aCGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCCTGGAAGCCCTgg                                                                                                                                          <  1:52353/88‑2 (MQ=255)
                                                                       aCGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCCTGGAAGCCCTgg                                                                                                                                          >  2:52353/1‑87 (MQ=255)
                                                                          ccAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACa                                                                          <  8:151637/149‑1 (MQ=255)
                                                                                cGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCaaa                                                                     <  3:127937/148‑1 (MQ=255)
                                                                                cGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAAc                                                                    >  4:122976/1‑149 (MQ=255)
                                                                                  cAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGc                                                                                                                                               <  4:106922/72‑1 (MQ=255)
                                                                                  cAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGc                                                                                                                                               >  3:106922/1‑72 (MQ=255)
                                                                                       aGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGACGCCACGGTCAGCAGCAGTCGCGGCTAGGCGTCATCGATGATCGACAGCAAACGCTCCTg                                                             >  8:39047/1‑149 (MQ=255)
                                                                                         ggCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGt                                                           >  1:266739/1‑149 (MQ=255)
                                                                                           ccGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGGTGATCGACAGCAAACGCTCCTGGtgc                                                         >  1:76567/1‑149 (MQ=255)
                                                                                             gggTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGACGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGtgctg                                                       >  2:251132/1‑149 (MQ=255)
                                                                                                 ccAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTgg                                                   >  3:240539/1‑149 (MQ=255)
                                                                                                  cAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGgc                                                  >  2:158518/1‑149 (MQ=255)
                                                                                                       cATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGcgcgcg                                             <  6:115029/149‑1 (MQ=255)
                                                                                                       cATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGcgcgcg                                             >  1:326988/1‑149 (MQ=255)
                                                                                                        aTCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCgg                                            <  2:326988/149‑1 (MQ=255)
                                                                                                         tCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGa                                           >  2:204590/1‑149 (MQ=255)
                                                                                                          ccACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGAc                                          >  7:111818/1‑149 (MQ=255)
                                                                                                           cACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGAc                                          <  4:240539/148‑1 (MQ=255)
                                                                                                             cTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTcc                                       <  3:122976/149‑1 (MQ=255)
                                                                                                                ccAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATc                                                                              <  7:108685/107‑1 (MQ=255)
                                                                                                                ccAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATc                                                                              >  8:108685/1‑107 (MQ=255)
                                                                                                                ccAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCcggc                                    >  3:300484/1‑149 (MQ=255)
                                                                                                                ccAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCcggc                                    >  8:261345/1‑149 (MQ=255)
                                                                                                                    cagTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCTGGt                                <  5:199021/149‑1 (MQ=255)
                                                                                                                    cagTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGt                                <  7:6078/149‑1 (MQ=255)
                                                                                                                       ttGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAgg                             <  5:46480/149‑1 (MQ=255)
                                                                                                                               gCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGCTCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCAc                      >  2:255706/1‑149 (MQ=255)
                                                                                                                                cGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGGTGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCAcc                     <  2:76567/149‑1 (MQ=255)
                                                                                                                                  ttGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACcgc                   <  2:266739/149‑1 (MQ=255)
                                                                                                                                     tCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGAt                <  8:111818/149‑1 (MQ=255)
                                                                                                                                         ccgccAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACg            >  1:170648/1‑149 (MQ=255)
                                                                                                                                             cAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGTAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCg        >  2:104528/1‑149 (MQ=255)
                                                                                                                                              aGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAg                              >  5:149922/1‑126 (MQ=255)
                                                                                                                                              aGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAg                              <  6:149922/126‑1 (MQ=255)
                                                                                                                                              aGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGAt                <  4:82520/140‑1 (MQ=255)
                                                                                                                                              aGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGAt                >  3:82520/1‑140 (MQ=255)
                                                                                                                                                   ttCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCGCATa    >  8:131264/1‑147 (MQ=255)
                                                                                                                                                   ttCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCGCATa    >  4:162725/1‑147 (MQ=255)
                                                                                                                                                     ccAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCGCATACa  <  7:39047/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCCTTGGGCAGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCGCATACA  >  NC_002947/4830642‑4830936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: