Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,783,744 A→G 10.2% G1205G (GGT→GGC pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,783,7440AG10.2% 160.3 / 7.6 59G1205G (GGT→GGCpvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base G (3/3);  total (32/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 = 6.81e-01

CATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGGTGTCGTACGCG  >  NC_002947/4783595‑4783892
                                                                                                                                                     |                                                                                                                                                    
cATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGa                                                                                                                                                      <  1:350966/150‑1 (MQ=255)
cATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGa                                                                                                                                                      <  2:218344/150‑1 (MQ=255)
   tGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCa                                                                                                                                                   <  2:1111623/150‑1 (MQ=255)
       gcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCGGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACGGCCGTTGCTGACCATCgg                                                                                                                                               >  2:1002819/1‑150 (MQ=255)
          ggggTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGa                                                                                                                                            <  2:312123/150‑1 (MQ=255)
              tGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTc                                                                                                                                        <  1:167505/150‑1 (MQ=255)
                ccGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGCTGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGa                                                                                                                                      <  1:249198/150‑1 (MQ=255)
                 cGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGaa                                                                                                                                     <  1:872824/150‑1 (MQ=255)
                 cGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGaa                                                                                                                                     <  1:361289/150‑1 (MQ=255)
                     cTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAAt                                                                                                                                 >  1:145479/1‑150 (MQ=255)
                           ggCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGcc                                                                                                                           >  1:962795/1‑150 (MQ=255)
                             cAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCt                                                                                                                         <  2:223405/150‑1 (MQ=255)
                             cAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCt                                                                                                                         >  1:28476/1‑150 (MQ=255)
                                       gCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCg                                                                                                               >  2:108403/1‑150 (MQ=255)
                                             ccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACgg                                                                                                         <  1:108403/150‑1 (MQ=255)
                                              cTAGGCTAGCCTGCTGGTCGGCCGTCTGCAGGTCTCCGAACAGGATTCGCCAAGACACCCCATCCACCGCCAGGTGATGCAACGCCAGCAACAGCCGTTGTTGACCATCGGCGAGATCGAACAATGACGCCCTGAACAGCTCAACACggg                                                                                                        >  2:260806/1‑150 (MQ=255)
                                               tAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAg                                                                                                                   <  1:833753/138‑1 (MQ=255)
                                               tAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAg                                                                                                                   >  2:833753/1‑138 (MQ=255)
                                                  gCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg                                                                                                    >  2:355890/1‑150 (MQ=255)
                                                  gCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCACAGCCTATTTTTTTTTTACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg                                                                                                    >  2:98511/1‑150 (MQ=255)
                                                    tAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGt                                                                                                  >  2:390738/1‑150 (MQ=255)
                                                        ctgctgGTAGGCCGTCTGCAGGTCTTCGAATAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGAAAGGTCGat                                                                                              >  2:357060/1‑149 (MQ=255)
                                                        ctgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAAc                                                                                                                                    >  2:249787/1‑112 (MQ=255)
                                                        ctgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAAc                                                                                                                                    <  1:249787/112‑1 (MQ=255)
                                                                          cAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGtt                                                                            <  1:355890/150‑1 (MQ=255)
                                                                             gTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCAc                                                                                  >  1:621245/1‑141 (MQ=255)
                                                                             gTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCAc                                                                                  <  2:621245/141‑1 (MQ=255)
                                                                                ttCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGCGGCTGGCCTGCACCTCGTTggcggc                                                                      >  2:635638/1‑150 (MQ=255)
                                                                                ttCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTggcggc                                                                      >  2:213080/1‑150 (MQ=255)
                                                                                    aaCAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGtt                                                                  >  1:534060/1‑150 (MQ=255)
                                                                                      cAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTc                                                                <  1:215943/150‑1 (MQ=255)
                                                                                             cgcCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCt                                                         <  1:213080/150‑1 (MQ=255)
                                                                                                cAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCa                                                      >  1:1009226/1‑150 (MQ=255)
                                                                                                cAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCa                                                      >  1:616703/1‑150 (MQ=255)
                                                                                                    acacCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGc                                                                                    <  1:279340/116‑1 (MQ=255)
                                                                                                    acacCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGc                                                                                    >  2:279340/1‑116 (MQ=255)
                                                                                                        cccATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTc                                              >  1:94542/1‑150 (MQ=255)
                                                                                                        cccATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGCTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTc                                              <  1:357060/150‑1 (MQ=255)
                                                                                                          cATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTCGGCTTGATGTTCTAGCGCTGCAAGGCTGTCCa                                            >  1:840237/1‑150 (MQ=255)
                                                                                                           ctCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCag                                           <  2:842572/149‑1 (MQ=255)
                                                                                                            tCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGCGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCaga                                          <  2:616703/150‑1 (MQ=255)
                                                                                                              caccaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGt                                        <  2:1016839/150‑1 (MQ=255)
                                                                                                               accaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt                                                                                           <  2:568266/98‑1 (MQ=255)
                                                                                                               accaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt                                                                                           >  1:568266/1‑98 (MQ=255)
                                                                                                               accaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGt                                        <  2:621263/149‑1 (MQ=255)
                                                                                                               accaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGt                                        >  1:621263/1‑149 (MQ=255)
                                                                                                                         gATACACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCa                             >  2:218781/1‑150 (MQ=255)
                                                                                                                            gCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCagga                          >  1:1039490/1‑150 (MQ=255)
                                                                                                                              aCCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTGCGGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg                        >  1:57783/1‑150 (MQ=255)
                                                                                                                              aCCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg                        <  2:28476/150‑1 (MQ=255)
                                                                                                                                cGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACATCGAGGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTGGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGc                      <  1:62630/150‑1 (MQ=255)
                                                                                                                                  tcAGCAACGGCTGTTGCTGGCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTa                                                              <  1:695506/107‑1 (MQ=255)
                                                                                                                                  tcAGCAACGGCTGTTGCTGGCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTa                                                              >  2:695506/2‑108 (MQ=255)
                                                                                                                                  tcAGCAACGGCTGTTGCTGGCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTc                                              >  1:1055314/2‑124 (MQ=255)
                                                                                                                                  tcAGCAACGGCTGTTGCTGGCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTc                                              >  1:1055234/2‑124 (MQ=255)
                                                                                                                                  tcAGCAACGGCTGTTGCTGGCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTc                                              <  2:1055314/123‑1 (MQ=255)
                                                                                                                                  tcAGCAACGGCTGTTGCTGGCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTc                                              <  2:1055234/123‑1 (MQ=255)
                                                                                                                                  tcAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCATGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCt                    >  2:606672/2‑150 (MQ=255)
                                                                                                                                         cAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCgg             >  1:221517/1‑150 (MQ=255)
                                                                                                                                                   tGACCATCGGCGAGTTCGAAGAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCAATATTGACTACCTGCAGGAGGGCCTGCACCGGTGTCGTAcgc   >  1:561454/1‑150 (MQ=255)
                                                                                                                                                    gACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGGTGTCGTAcgcg  <  1:390738/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGGTGTCGTACGCG  >  NC_002947/4783595‑4783892

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