Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,139 T→G 8.0% K144T (AAG→ACG)  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,1390TG 8.0% 163.6 / 2.4 58K144T (AAG→ACG) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base G (3/2);  total (30/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 = 1.00e+00

GTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTG  >  NC_002947/273991‑274288
                                                                                                                                                     |                                                                                                                                                     
gTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCNGTGNCCNTGTAGCCCACGCNGGGCTTGAGGGTGATCAACCCGGCCtt                                                                                                                                                      <  1:1000543/150‑1 (MQ=255)
    cccTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGcc                                                                                                                                                  >  1:56793/1‑150 (MQ=255)
      cTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcag                                                                                                                                                >  1:511672/1‑150 (MQ=255)
        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:121891/1‑150 (MQ=255)
                                                                                                            cTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGcc                                          >  2:901334/1‑150 (MQ=255)
                                                                                                               gTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAAc                                       >  2:429920/1‑150 (MQ=255)
                                                                                                                   cccACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCttt                                                                <  1:1196140/121‑1 (MQ=255)
                                                                                                                   cccACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCttt                                                                >  2:1196140/1‑121 (MQ=255)
                                                                                                                     cACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGa                                                                                                       >  1:1065056/1‑80 (MQ=255)
                                                                                                                     cACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGa                                                                                                       <  2:1065056/80‑1 (MQ=255)
                                                                                                                     cACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttgtt                                 <  2:745912/150‑1 (MQ=255)
                                                                                                                       cGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtga                               <  1:206237/150‑1 (MQ=255)
                                                                                                                          cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                >  1:1150824/1‑130 (MQ=255)
                                                                                                                          cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                <  2:1150824/130‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
GTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTG  >  NC_002947/273991‑274288

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: