Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,148 A→T 6.9% L141Q (CTG→CAG)  PP_0221 ← methionine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,1480AT 6.9% 178.5 / 4.0 58L141Q (CTG→CAG) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base T (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGC  >  NC_002947/273999‑274292
                                                                                                                                                      |                                                                                                                                                
gCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagca                                                                                                                                                  <  2:763172/150‑1 (MQ=255)
       gTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGTATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagcagcag                                                                                                                                           >  2:71454/1‑150 (MQ=255)
                 tcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCg                                                                                                                                 <  1:910211/150‑1 (MQ=255)
                  cgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGc                                                                                                                                >  2:457921/1‑150 (MQ=255)
                   gtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCt                                                                                                                               <  1:457921/150‑1 (MQ=255)
                     aGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGt                                                                                                                             <  2:1207885/150‑1 (MQ=255)
                     aGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGt                                                                                                                             <  2:1186840/150‑1 (MQ=255)
                     aGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGt                                                                                                                             <  1:536427/150‑1 (MQ=255)
                       gcgcGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTg                                                                                                                           >  2:638877/1‑150 (MQ=255)
                          cgTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATc                                                                                                                        <  2:511672/150‑1 (MQ=255)
                               cAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTc                                                                                                                   <  2:56793/150‑1 (MQ=255)
                                   tGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTcctaa                                                                                                               <  1:948175/150‑5 (MQ=255)
                                   tGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTg                                                                                                               <  1:551585/150‑1 (MQ=255)
                                     cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACAGCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCgggt                                                                                                                              <  1:105484/134‑4 (MQ=25)
                                     cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACAGCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCgggt                                                                                                                              >  2:105484/1‑131 (MQ=25)
                                       ccGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGa                                                                                                           >  1:613824/1‑150 (MQ=255)
                                        cGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGAt                                                                                                          >  2:874499/1‑150 (MQ=255)
                                           cTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGc                                                                                                       >  1:513329/1‑150 (MQ=255)
                                                       gATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGcc                                                                                           >  2:256692/1‑150 (MQ=255)
                                                            tCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGCAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAg                                                                                      >  1:623399/1‑150 (MQ=255)
                                                             cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                                                                                     >  1:122502/1‑150 (MQ=255)
                                                             cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                                                                                     >  1:745912/1‑150 (MQ=255)
                                                                cTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAg                                                                                  >  2:206237/1‑150 (MQ=255)
                                                                 tttttAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGcc                                                                                           >  1:433106/1‑140 (MQ=255)
                                                                 tttttAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGcc                                                                                           <  2:433106/140‑1 (MQ=255)
                                                                    ttAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCg                                                                              <  1:256692/150‑1 (MQ=255)
                                                                               gAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCtttt                                                                   <  1:638877/150‑1 (MQ=255)
                                                                                          ccTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCCGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCGGGCAGCTCGGCAAAGCTTTTGTACTTTTTCg                                                        >  2:528096/1‑150 (MQ=255)
                                                                                           cTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGa                                                       <  2:613824/150‑1 (MQ=255)
                                                                                                tGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGACACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCGTTCAGGTAGCTAGGCAAGGATTTTGTACTTTTTCGCGTaga                                                  >  2:437181/1‑148 (MQ=255)
                                                                                                tGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATAAACTCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACa                                                  <  1:369042/150‑1 (MQ=255)
                                                                                                 ggCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAg                                                 <  1:71454/150‑1 (MQ=255)
                                                                                                   ccTTCTAGCCCACGCCCGCCTTGAGGGTGATCAACCCGGCCTTCGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGc                                               <  1:437181/150‑1 (MQ=255)
                                                                                                    cTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGcc                                              >  2:901334/1‑150 (MQ=255)
                                                                                                    cTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGcc                                              >  2:121891/1‑150 (MQ=255)
                                                                                                       gTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAAc                                           >  2:429920/1‑150 (MQ=255)
                                                                                                           cccACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCttt                                                                    >  2:1196140/1‑121 (MQ=255)
                                                                                                           cccACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCttt                                                                    <  1:1196140/121‑1 (MQ=255)
                                                                                                             cACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGa                                                                                                           <  2:1065056/80‑1 (MQ=255)
                                                                                                             cACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGa                                                                                                           >  1:1065056/1‑80 (MQ=255)
                                                                                                             cACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttgtt                                     <  2:745912/150‑1 (MQ=255)
                                                                                                               cGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtga                                   <  1:206237/150‑1 (MQ=255)
                                                                                                                  cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                    <  2:1150824/130‑1 (MQ=255)
                                                                                                                  cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                    >  1:1150824/1‑130 (MQ=255)
                                                                                                                  cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgatga                                >  1:742167/1‑150 (MQ=255)
                                                                                                                     gCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCaaa                                                                         <  2:221451/106‑1 (MQ=255)
                                                                                                                     gCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCaaa                                                                         >  1:221451/1‑106 (MQ=255)
                                                                                                                          aGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                        <  1:7430/150‑1 (MQ=255)
                                                                                                                          aGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                        <  1:121891/150‑1 (MQ=255)
                                                                                                                            ggTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGc                      >  1:329750/1‑150 (MQ=255)
                                                                                                                                  caacagcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCGGGTTGATCAGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAgg                >  1:1032000/7‑150 (MQ=25)
                                                                                                                                     cccGGCCTTGGCTAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGCAGTTCTGGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcac             <  1:74221/150‑1 (MQ=255)
                                                                                                                                       ccGGCCGTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGTGTCGTTGGCGATCCCGTCGATCGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcaca            >  2:644891/1‑150 (MQ=255)
                                                                                                                                        cGGCCTTGGCCAGCAGCAGCACGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACaa           >  1:755963/1‑150 (MQ=255)
                                                                                                                                           ccTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGt        >  2:263668/1‑150 (MQ=255)
                                                                                                                                            cTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGtt       <  2:513329/150‑1 (MQ=255)
                                                                                                                                             ttGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTg      >  1:622497/1‑150 (MQ=255)
                                                                                                                                               ggCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCTAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGaa    <  1:429920/150‑1 (MQ=255)
                                                                                                                                               ggCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGaa    <  1:874499/150‑1 (MQ=255)
                                                                                                                                                 ccAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGc  <  2:329750/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGC  >  NC_002947/273999‑274292

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