Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,882,857 A→G 6.6% intergenic (+56/‑169) oprN → / → PP_3428 multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,882,8570AG 6.6% 212.7 / 4.0 76intergenic (+56/‑169)oprN/PP_3428multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (2/3);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

GGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCC  >  NC_002947/3882712‑3883002
                                                                                                                                                 |                                                                                                                                                 
tgTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGcc                                                                                                                                               >  2:472985/2‑150 (MQ=255)
 gTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGCCGGCTGGCAACCGGGGAGGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCCACCGTTTGCCCCCCAGGCCa                                                                                                                                              >  2:737257/1‑150 (MQ=255)
  tGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCAt                                                                                                                                             <  1:1327809/150‑1 (MQ=255)
       tCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGccc                                                                                                                                        >  1:287221/1‑150 (MQ=255)
            cTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCggg                                                                                                                                   >  2:327477/1‑150 (MQ=255)
                 tCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCGGGGcaaa                                                                                                                              <  2:1495419/150‑4 (MQ=255)
                 tCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCttt                                                                                                                              >  2:565270/1‑150 (MQ=255)
                 tCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTCGCCCCGCGGGCCATGGCCTGCGGGGcaaa                                                                                                                              <  2:1345347/150‑4 (MQ=255)
                        aTCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTCGCCAGCGCGCACAGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGTGGGGCTTTTTTATGt                                                                                                                       >  2:1052460/1‑150 (MQ=255)
                          cGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGttt                                                                                                                     <  2:1591671/150‑1 (MQ=255)
                           gTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTg                                                                                                                    >  1:1899711/1‑150 (MQ=255)
                              tcGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCa                                                                                                                 >  2:1134155/1‑150 (MQ=255)
                                 aTCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCACAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCg                                                                                                              >  1:1535543/1‑150 (MQ=255)
                                  tCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTCTGTTTGCCATCGt                                                                                                             >  2:1029950/1‑150 (MQ=255)
                                  tCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGt                                                                                                             <  2:69821/150‑1 (MQ=255)
                                   cTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGTCGGTGGCTAGCGCGCACTGATTCCACTGCCTCAACGGCTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCTTCGTa                                                                                                            <  2:1832833/150‑1 (MQ=255)
                                    tACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTAt                                                                                                           <  1:1905588/150‑1 (MQ=255)
                                    tACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTAt                                                                                                           >  2:787187/1‑150 (MQ=255)
                                     aCAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACGGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGtttt                                                                                                          >  2:462993/1‑150 (MQ=255)
                                      cAAGGCCCTGGGTGGCGGCTGGCAACCGGCGACGGTGGCCAGCGCGCACTGATCCCACTGCCTCAACGACTCCTTTGGTTGGCCCCTCCCCCAACCGTTTGCCCCGCAGGGCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTc                                                                                                         <  2:389887/150‑1 (MQ=255)
                                      cAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTc                                                                                                         >  1:301006/1‑150 (MQ=255)
                                        aGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAc                                                                                                       <  1:900859/150‑1 (MQ=255)
                                          gCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                     <  1:325805/150‑1 (MQ=255)
                                          gCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                     <  2:578366/150‑1 (MQ=255)
                                                    gcggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAg                                                                                           <  1:1854903/150‑1 (MQ=255)
                                                      ggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                             >  2:1544426/1‑114 (MQ=255)
                                                      ggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                             <  1:1544426/114‑1 (MQ=255)
                                                         tGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCAt                                                                                      >  2:357613/1‑150 (MQ=255)
                                                         tGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCAt                                                                                      <  2:1899711/150‑1 (MQ=255)
                                                             aaCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGGCTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGaa                                                                                  <  2:287221/150‑1 (MQ=255)
                                                             aaCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGaa                                                                                  >  1:602260/1‑150 (MQ=255)
                                                             aaCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACCCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATCCGCAGTGCATGGaa                                                                                  >  2:785454/1‑150 (MQ=255)
                                                              aCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAAt                                                                                 <  1:327477/150‑1 (MQ=255)
                                                               ccGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCAATACGCAGTGCATGGAATGc                                                                               >  1:983980/1‑150 (MQ=255)
                                                               ccGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTACTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATg                                                                                >  1:1129036/1‑150 (MQ=255)
                                                                  gagaCGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGcc                                                                              <  2:1363629/149‑1 (MQ=255)
                                                                  gagaCGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGcc                                                                              >  1:1363629/1‑149 (MQ=255)
                                                                       ggTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCa                                                                        >  1:472824/1‑150 (MQ=255)
                                                                        gTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCATGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAg                                                                       <  1:1066447/150‑1 (MQ=255)
                                                                         tGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAgg                                                                      >  2:1598157/1‑150 (MQ=255)
                                                                                   gcACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggcgcct                                                            <  1:357613/150‑1 (MQ=255)
                                                                                       tGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAc                                                                                                       <  1:1362259/103‑1 (MQ=255)
                                                                                       tGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAc                                                                                                       >  2:1362259/1‑103 (MQ=255)
                                                                                           tCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACg                                                    <  2:1535543/150‑1 (MQ=255)
                                                                                                tGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTGATGTTTGCCATCGTATTCACCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGc                                              <  2:983980/150‑1 (MQ=255)
                                                                                                    tCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcg                                           <  2:602260/150‑1 (MQ=255)
                                                                                                     cAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcgg                                          >  2:570938/1‑150 (MQ=255)
                                                                                                       aCGACTCCTTTGGTTGGCTCCTCCCCCACCCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggtt                                        <  1:1598157/150‑1 (MQ=255)
                                                                                                            tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCaa                                                                                                  >  1:1011533/1‑87 (MQ=255)
                                                                                                            tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcggtt                                   >  2:1156934/1‑150 (MQ=255)
                                                                                                            tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCCGGGCTTTTTTATGTTTGCCATCGTATTCACCCCaa                                                                                                  <  2:1011533/87‑1 (MQ=255)
                                                                                                               tttGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCa                                                                                                        >  1:411486/1‑78 (MQ=255)
                                                                                                               tttGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCa                                                                                                        <  2:411486/78‑1 (MQ=255)
                                                                                                                   gttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttc                                       >  2:1230971/1‑139 (MQ=255)
                                                                                                                   gttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttc                                       <  1:1230971/139‑1 (MQ=255)
                                                                                                                   gttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTa                            >  1:538966/1‑150 (MQ=255)
                                                                                                                   gttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTAATCAGCCGAATACGCAATACCACACGTCCGAACTCCAGTCACGGTCCTTGCCAACGGCGTGCCCGGTTCGGTTGCATTTa                            >  1:1868523/1‑150 (MQ=255)
                                                                                                                    ttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTa                            >  2:1119951/1‑149 (MQ=255)
                                                                                                                    ttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTa                            <  1:1119951/149‑1 (MQ=255)
                                                                                                                           ctccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCgg                    >  1:506447/1‑150 (MQ=255)
                                                                                                                           ctccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTATTATGTTTGTTATAGCATTCAACCAATTTTGGAGTGCATGGAATGCCACCGCAGGCGTGGCGCCTGTCCAACGGCGTGCCCGGTTCGGTTGCATTTACACCCCgg                    >  1:425650/1‑150 (MQ=255)
                                                                                                                              ccccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTg                 >  1:1009238/1‑150 (MQ=255)
                                                                                                                              ccccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTg                 >  1:1135598/1‑150 (MQ=255)
                                                                                                                                ccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGc               >  1:551810/1‑150 (MQ=255)
                                                                                                                                 cAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGcc              <  2:538966/150‑1 (MQ=255)
                                                                                                                                    ccGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCaa           <  1:565270/150‑1 (MQ=255)
                                                                                                                                    ccGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCaa           >  1:1347993/1‑150 (MQ=255)
                                                                                                                                    ccGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCaa           <  1:787187/150‑1 (MQ=255)
                                                                                                                                     cGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAg          <  2:1009238/150‑1 (MQ=255)
                                                                                                                                     cGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAg          <  2:1135598/150‑1 (MQ=255)
                                                                                                                                          gCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCACTACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTg     <  1:1052460/150‑1 (MQ=255)
                                                                                                                                           ccccGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCAACTCCAGTCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGc    >  1:1277672/1‑150 (MQ=255)
                                                                                                                                           ccccGCAGGCCATGGCTCGCGGGGCTTTTTTATGTTTGCCTTCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGTCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGc    <  2:69366/150‑1 (MQ=255)
                                                                                                                                             ccGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcggtt                                   >  1:903333/1‑117 (MQ=255)
                                                                                                                                             ccGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcggtt                                   <  2:903333/117‑1 (MQ=255)
                                                                                                                                             ccGCAGGCCCTGGACAACGGGGCTGTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGccc  >  2:465997/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCC  >  NC_002947/3882712‑3883002

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