Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,830,788 T→C 9.8% E134E (GAA→GAG) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,830,7880TC 9.8% 163.8 / 7.5 62E134E (GAA→GAG) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base C (2/4);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CCCTTGGGCAGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCGCATA  >  NC_002947/4830642‑4830934
                                                                                                                                                   |                                                                                                                                                  
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)
                                                                                                                                                   |                                                                                                                                                  
CCCTTGGGCAGCGCGGTGGAACCGGAGGTGTACTGCAGGAAGGCGATGTCATCGCCCTTGAGCACGGGTTCACGCCAGTTCGCAGCCAGGGCCGGGTCCAAGCCATCCACTGCCAGCAGTTGGGGGGCGTTGTCCGCCGCCAGGGCTTCCAGGCCCTGCAGGCTGTCGCGCAGCGCCGCCACGGTCAGCAGCAGTCGCGGCTCGGCGTCATCGATGATCGACAGCAAACGCTCCTGGTGCTGCTGGCGCGCGGACTCCGGCGGGTAGGCCGGCACCGCGATCACGCCCGCATA  >  NC_002947/4830642‑4830934

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: