Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,776,087 T→C 5.9% Q1060Q (CAA→CAG) ‡ pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,776,0870TC 5.9% 203.3 / 2.2 68Q1060Q (CAA→CAG) ‡pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (39/25);  new base C (2/2);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCATACGACAGCGCCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGTCGGCGCCT  >  NC_002947/4775941‑4776235
                                                                                                                                                  |                                                                                                                                                    
gCATACGACAGCGCCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTg                                                                                                                                                    <  2:142362/149‑1 (MQ=255)
  aTACGACAGCGCCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCACCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCt                                                                                                                                                  >  4:185798/1‑149 (MQ=255)
     cGACAGCGCCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTgg                                                                                                                                                <  4:243405/148‑1 (MQ=255)
            gcCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGCAACTGgcgctggaa                                                                                                                                        <  2:158900/149‑6 (MQ=255)
            gcCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGCAACTGgcgct                                                                                                                                            >  5:275339/1‑144 (MQ=255)
            gcCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGCAACTGgcgct                                                                                                                                            <  6:275339/145‑2 (MQ=255)
                ggggCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCAcgc                                                                                                                                         >  5:330881/1‑144 (MQ=255)
                ggggCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCAcgc                                                                                                                                         <  6:330881/144‑1 (MQ=255)
                 gggCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCaa                                                                                                                                   >  6:336306/1‑149 (MQ=255)
                     tctcCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGgc                                                                                                                                               >  2:347174/1‑133 (MQ=255)
                     tctcCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGgc                                                                                                                                               <  1:347174/133‑1 (MQ=255)
                       tcCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAatcat                                                                                                                               <  8:242312/147‑1 (MQ=255)
                       tcCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAatcat                                                                                                                               <  5:336306/147‑1 (MQ=255)
                       tcCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCt                                                                                                                             >  2:184618/1‑149 (MQ=255)
                         cccGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCt                                                                                                                             <  1:184618/147‑1 (MQ=255)
                          ccGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCg                                                                                                                          >  3:124023/1‑149 (MQ=255)
                                caccacAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGAcgc                                                                                                                      >  7:199029/1‑147 (MQ=255)
                                   cacaGCGATGGCCGGCAGCTGCGCATCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGc                                                                                                                  >  3:240071/1‑148 (MQ=255)
                                   cacaGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGc                                                                                                                  >  1:93197/1‑148 (MQ=255)
                                                 gCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCAcc                                                                                                    >  1:297221/1‑148 (MQ=255)
                                                 gCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCAcc                                                                                                    >  4:164659/1‑148 (MQ=255)
                                                     cTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGtt                                                                                               <  2:297221/149‑1 (MQ=255)
                                                            cccGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCgaagaa                                                                                        <  5:297753/149‑1 (MQ=255)
                                                                  gtGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTc                                                                                  <  1:217342/149‑1 (MQ=255)
                                                                     gCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCt                                                                               >  8:34421/1‑149 (MQ=255)
                                                                           gccACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGc                                                                          >  6:208798/1‑148 (MQ=255)
                                                                             cACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCaa                                                                       >  6:38350/1‑149 (MQ=255)
                                                                              aCCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAAc                                                                      <  4:11641/149‑1 (MQ=255)
                                                                               cccGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAAcc                                                                     <  7:34421/149‑1 (MQ=255)
                                                                               cccGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCCCAAGCCAAcc                                                                     <  2:11417/149‑1 (MQ=255)
                                                                                  gCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGAcgc                                                                                                                      >  6:332290/1‑97 (MQ=255)
                                                                                  gCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGAcgc                                                                                                                      <  5:332290/97‑1 (MQ=255)
                                                                                   cTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGcc                                                                 <  3:164659/149‑1 (MQ=255)
                                                                                    tCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGccc                                                                >  7:328847/1‑149 (MQ=255)
                                                                                           aGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCtt                                                          <  8:206470/148‑1 (MQ=255)
                                                                                           aGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCtt                                                          >  7:206470/1‑148 (MQ=255)
                                                                                           aGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTc                                                         <  5:38350/149‑1 (MQ=255)
                                                                                               cGCCAGGTCCTGGGCTTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGc                                                      >  5:220706/1‑148 (MQ=255)
                                                                                               cGCCAGGTCCTGGGCTTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGc                                                      <  6:220706/148‑1 (MQ=255)
                                                                                                  cAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACAt                                                  <  2:58755/149‑1 (MQ=255)
                                                                                                    ggTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCg                                                >  4:296278/1‑149 (MQ=255)
                                                                                                          gggCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGCCCTTCAGCACATCGGcc                                             >  2:121708/1‑146 (MQ=255)
                                                                                                          gggCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGcc                                             >  6:51035/1‑146 (MQ=255)
                                                                                                           ggCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCa                                           <  3:185798/147‑1 (MQ=255)
                                                                                                                 gAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCAcgc                                                                                                                                         <  5:144542/47‑1 (MQ=255)
                                                                                                                 gAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCAcgc                                                                                                                                         <  5:15043/47‑1 (MQ=255)
                                                                                                                 gAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCAcgc                                                                                                                                         >  6:144542/1‑47 (MQ=255)
                                                                                                                 gAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCAcgc                                                                                                                                         >  6:15043/1‑47 (MQ=255)
                                                                                                                 gAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATcgccg                                   >  2:131325/1‑149 (MQ=255)
                                                                                                                      cacTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAgc                              <  8:199029/149‑1 (MQ=255)
                                                                                                                       acTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAgcg                             >  7:81483/1‑149 (MQ=255)
                                                                                                                         tGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCt                           >  4:53102/1‑149 (MQ=255)
                                                                                                                         tGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCt                           >  2:69002/1‑149 (MQ=255)
                                                                                                                            gcAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTtgc                        >  1:367707/1‑149 (MQ=255)
                                                                                                                              aGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTtgctg                      >  1:351051/1‑149 (MQ=255)
                                                                                                                                cGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCa                    <  4:124023/149‑1 (MQ=255)
                                                                                                                                  gCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGc                  >  8:264207/1‑149 (MQ=255)
                                                                                                                                      ttCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGc              >  2:363410/1‑149 (MQ=255)
                                                                                                                                      ttCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGc              >  3:282196/1‑149 (MQ=255)
                                                                                                                                       gcgtgcgcCAGCAACTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGGCCTTCAGCACATCGGCCCAGATCGCCGCCAgc                              >  6:238683/5‑132 (MQ=255)
                                                                                                                                        ccAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGCCCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTg            >  5:333615/1‑149 (MQ=255)
                                                                                                                                         cAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTgg           >  2:195317/1‑149 (MQ=255)
                                                                                                                                          aGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGt          <  2:93197/149‑1 (MQ=255)
                                                                                                                                          aGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGt          >  3:288243/1‑149 (MQ=255)
                                                                                                                                            cgcCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGTCg        >  2:175558/1‑149 (MQ=255)
                                                                                                                                            cgcCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGTCg        >  8:174709/1‑149 (MQ=255)
                                                                                                                                              ccAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGTCGgc      >  8:245666/1‑149 (MQ=255)
                                                                                                                                                 gTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGTCGGCGcc   >  4:318109/1‑149 (MQ=255)
                                                                                                                                                  ttGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGTCGGCGCCt  <  1:69002/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCATACGACAGCGCCAGGGGCTCTCCCCGTGGCACCACAGCGATGGCCGGCAGCTGCGCACCCGTCGTGGCAGCCGCCACCCGCTCGGCGAAGCCCGCCAGGTCCTGGGCCTCGAAGACACTGCGCAGCGGCACTTCCAGCGCCAGTTGCTGGCGCACGCGGCCAATCATCTGCGACGCGAGCAAGGAATGGCCACCCAGTTCGAAGAAGCTGTCGCTCAAGCCAACCTGCCCGACCTTCAGCACATCGGCCCAGATCGCCGCCAGCGCTTGCTGCAGCGGGCTGGTCGGCGCCT  >  NC_002947/4775941‑4776235

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: ‑