Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,783,733 A→G 11.7% L1209P (CTG→CCG)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,783,7330AG11.7% 151.6 / 10.3 60L1209P (CTG→CCG) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base G (4/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

GCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG  >  NC_002947/4783584‑4783881
                                                                                                                                                     |                                                                                                                                                    
gCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCNACa                                                                                                                                                      <  2:696523/150‑1 (MQ=255)
        gcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGc                                                                                                                                              >  2:872824/1‑150 (MQ=255)
           cATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGa                                                                                                                                           <  2:218344/150‑1 (MQ=255)
           cATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGa                                                                                                                                           <  1:350966/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                                                                                                        >  2:833753/1‑138 (MQ=255)
                                                          tAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAg                                                                                                        <  1:833753/138‑1 (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                                                                                >  1:568266/1‑98 (MQ=255)
                                                                                                                          accaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt                                                                                <  2:568266/98‑1 (MQ=255)
                                                                                                                          accaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGt                             >  1:621263/1‑149 (MQ=255)
                                                                                                                          accaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGt                             <  2:621263/149‑1 (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:1055234/123‑1 (MQ=255)
                                                                                                                                             tcAGCAACGGCTGTTGCTGGCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTc                                   <  2:1055314/123‑1 (MQ=255)
                                                                                                                                             tcAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCATGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCt         >  2:606672/2‑150 (MQ=255)
                                                                                                                                                    cAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCgg  >  1:221517/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG  >  NC_002947/4783584‑4783881

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: