Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,529,131 A→T 6.6% N52Y (AAC→TAC)  ftsZ → cell division protein FtsZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,529,1310AT 6.6% 183.7 / 3.1 61N52Y (AAC→TAC) ftsZcell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base T (2/2);  total (31/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 = 3.98e-01

AGCTCGTAGACAACGTCCCGCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGT  >  NC_002947/1528985‑1529275
                                                                                                                                                  |                                                                                                                                                
aGCTCGTAGACAACGTCCCGCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAAc                                                                                                                                                <  8:201559/149‑1 (MQ=255)
 gCTCGTAGACAACGTCCCGCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGCCGCTGAAAAACa                                                                                                                                               >  8:31991/1‑149 (MQ=255)
  cTCGTAGACAACGTCCCGCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACa                                                                                                                                               <  7:31991/148‑1 (MQ=255)
     gTAGACAACGTCCCGCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTCCGCCAACACCGATGCCCAGGCGCTGAAAAACATTgg                                                                                                                                           >  4:182989/1‑149 (MQ=255)
           aaCGTCCCGCAAAGTCCGGTCTATAAAGTGATCGGCGTTGGTGGTGGTTGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTGATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGgcgcgcg                                                                                                                                     <  4:21816/149‑1 (MQ=255)
           aaCGTCCCGCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGgcgcgcg                                                                                                                                     >  4:306233/1‑149 (MQ=255)
                   gCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATcc                                                                                                                             >  4:159190/1‑149 (MQ=255)
                   gCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATcc                                                                                                                             >  3:69213/1‑149 (MQ=255)
                          ccGGTCATTAAAGTGATCGGCGTTGGTAGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAAtt                                                                                                                      >  1:15308/1‑149 (MQ=255)
                              tCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCa                                                                                                                                   <  3:311060/132‑1 (MQ=255)
                              tCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCa                                                                                                                                   >  4:311060/1‑132 (MQ=255)
                              tCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTggg                                                                                                                   >  5:89373/1‑148 (MQ=255)
                                aTTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCAc                                                                                                                                  <  5:191172/131‑1 (MQ=255)
                                aTTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCAc                                                                                                                                  >  6:191172/1‑131 (MQ=255)
                                 ttAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAATACATTGGCGCGCGCACCATCCTGCAATTGGGCACg                                                                                                               >  3:115597/1‑149 (MQ=255)
                                 ttAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAATACATTGGCGCGCGCACCATCCTGCAATTGGGCACg                                                                                                               >  3:120132/1‑149 (MQ=255)
                                 ttAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTg                                                                                                                           >  5:18669/1‑137 (MQ=255)
                                 ttAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTg                                                                                                                           <  6:18669/137‑1 (MQ=255)
                                 ttAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACg                                                                                                               >  5:35287/1‑149 (MQ=255)
                                  tAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATc                                                                                                                              <  7:186489/133‑1 (MQ=255)
                                  tAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGACCAGGCGCTGAAAAACATTGGCGCGCGCACCATc                                                                                                                              >  8:186489/1‑133 (MQ=255)
                                    aaGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGgt                                                                                                            <  3:306233/149‑1 (MQ=255)
                                          tCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCGAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGAcc                                                                                                      <  3:182989/149‑1 (MQ=255)
                                          tCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGAc                                                                                                       >  4:246847/1‑148 (MQ=255)
                                                           ggcggcAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACgg                                                                                                              <  6:146358/124‑1 (MQ=255)
                                                                   cGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATc                                                                             <  5:43512/149‑1 (MQ=255)
                                                                       gTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACgg                                                                                                              >  4:175454/1‑112 (MQ=255)
                                                                       gTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACgg                                                                                                              <  3:175454/112‑1 (MQ=255)
                                                                            tCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCg                                                                    >  7:87324/1‑149 (MQ=255)
                                                                            tCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCg                                                                    <  2:15308/149‑1 (MQ=255)
                                                                                   gTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                             <  3:159190/149‑1 (MQ=255)
                                                                                     cAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAgg                                                           >  3:69245/1‑149 (MQ=255)
                                                                                        gagCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAgg                                                                                                  >  3:126406/1‑107 (MQ=255)
                                                                                        gagCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAgg                                                                                                  <  4:126406/107‑1 (MQ=255)
                                                                                                tCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAATACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGaa                                                <  4:115597/149‑1 (MQ=255)
                                                                                                tCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAATACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGaa                                                <  4:120132/149‑1 (MQ=255)
                                                                                                 cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAAtt                                                                                                                      <  4:26968/78‑1 (MQ=255)
                                                                                                 cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAAtt                                                                                                                      >  3:26968/1‑78 (MQ=255)
                                                                                                    aGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTgg                                                  <  7:38203/143‑1 (MQ=255)
                                                                                                    aGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTgg                                                  >  8:38203/1‑143 (MQ=255)
                                                                                                              ttCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgccggtg                                                                                   <  5:19304/100‑1 (MQ=255)
                                                                                                              ttCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgccggtg                                                                                   >  6:19304/1‑100 (MQ=255)
                                                                                                                 aTCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAgcg                                     >  3:32966/1‑143 (MQ=255)
                                                                                                                 aTCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAgcg                                     <  4:32966/143‑1 (MQ=255)
                                                                                                                    tGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGa                                              >  5:293140/1‑131 (MQ=255)
                                                                                                                    tGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGa                                              <  6:293140/131‑1 (MQ=255)
                                                                                                                    tGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCg                             >  5:31426/1‑148 (MQ=255)
                                                                                                                         cAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGtgc                       <  6:89373/149‑1 (MQ=255)
                                                                                                                         cAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGtgc                       >  7:33016/1‑149 (MQ=255)
                                                                                                                             acCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGAACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCAGCGCTGGAACACCCTCACCGCATTGCCGAAGtgctg                     >  7:291416/1‑147 (MQ=255)
                                                                                                                               cGATGCCCAGGCGCTGAAAAACATTGGCGCGCTATCCATCCTGAAATTGGGCACGGGTGTGAGCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAgg                 <  4:69245/149‑1 (MQ=255)
                                                                                                                               cGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAgg                 >  7:269070/1‑149 (MQ=255)
                                                                                                                               cGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAgg                 <  4:69213/149‑1 (MQ=255)
                                                                                                                                  tGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCa              <  8:87324/149‑1 (MQ=255)
                                                                                                                                     ccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTg                                <  4:47964/128‑1 (MQ=255)
                                                                                                                                     ccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTg                                >  3:47964/1‑128 (MQ=255)
                                                                                                                                       aGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAAc         <  3:257645/149‑1 (MQ=255)
                                                                                                                                          cgcTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGGGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATg      <  4:46163/149‑1 (MQ=255)
                                                                                                                                           gcTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATgg     >  4:238237/1‑149 (MQ=255)
                                                                                                                                             tGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATg      <  8:95503/146‑1 (MQ=255)
                                                                                                                                             tGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATg      >  7:95503/1‑146 (MQ=255)
                                                                                                                                               aaaaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                             <  2:218051/89‑1 (MQ=255)
                                                                                                                                               aaaaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCa                                                             >  1:218051/1‑89 (MQ=255)
                                                                                                                                               aaaaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGgtgt  >  3:45735/1‑148 (MQ=255)
                                                                                                                                               aaaaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGgtg   <  3:246847/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AGCTCGTAGACAACGTCCCGCAAAGTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGT  >  NC_002947/1528985‑1529275

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