Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 134,397 A→C 7.9% D98A (GAT→GCT) ‡ PP_0129 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947134,3970AC 7.9% 172.4 / 3.6 60D98A (GAT→GCT) ‡PP_0129GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base C (3/2);  total (31/29)
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

AAGGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGCT  >  NC_002947/134249‑134543
                                                                                                                                                    |                                                                                                                                                  
aaGGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCACAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGa                                                                                                                                                    >  5:253263/1‑149 (MQ=255)
               tGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCAc                                                                                                                                      <  2:271123/148‑1 (MQ=255)
                 cATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTgg                                                                                                                                   <  7:148700/149‑1 (MQ=255)
                   tGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGtt                                                                                                                                 <  5:19522/149‑1 (MQ=255)
                     aaGGAAATGCGCGATGTCATTCGCAGCGTCAACATGGACGCCGGCCGCGCTGTCTTGATTCCGCGCCTTGAGAAAGGGGTACTGGATTGCGAACAGCGTCGCGAAACCCGCATGAACCTGGTCAGCGATGCGCTGACCGCACTGGTTgg                                                                                                                               <  8:90855/149‑1 (MQ=255)
                     aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCa                                                                                                                                       >  4:126326/1‑141 (MQ=255)
                     aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCa                                                                                                                                       <  3:126326/141‑1 (MQ=255)
                     aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCACATGAACCAGGTCAGCGCATCGCTGAcctgg                                                                                                                                       <  2:198350/141‑4 (MQ=255)
                     aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCACATGAACCAGGTCAGCGCATCGCTGAcctgg                                                                                                                                       >  1:198350/1‑138 (MQ=255)
                       ggAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATGCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGc                                                                                                                              >  1:84450/1‑148 (MQ=255)
                           aTGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCGTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACt                                                                                                                         <  4:21078/149‑1 (MQ=255)
                            tGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTg                                                                                                                        >  3:64717/1‑149 (MQ=255)
                               gcgATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGc                                                                                                                       >  4:151797/1‑147 (MQ=255)
                                      cATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTGCTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgc                                                                                                              >  5:243933/1‑149 (MQ=255)
                                        ttCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTGCTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgccc                                                                                                            <  6:243933/149‑1 (MQ=255)
                                        ttCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgccc                                                                                                            >  6:219474/1‑149 (MQ=255)
                                           gCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgcccctg                                                                                                         <  8:184944/149‑1 (MQ=255)
                                                    acaTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGat                                                                                                <  6:253263/149‑1 (MQ=255)
                                                      aTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGTCTGCCCCTGCCCAGCGata                                                                                               >  6:297331/1‑148 (MQ=255)
                                                      aTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGata                                                                                               >  3:220715/1‑148 (MQ=255)
                                                      aTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGata                                                                                               >  4:66225/1‑148 (MQ=255)
                                                        ggACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATAtct                                                                                            >  2:181544/1‑149 (MQ=255)
                                                              cGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGCCCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGc                                                                                      >  6:4853/1‑149 (MQ=255)
                                                              cGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGc                                                                                      >  2:48889/1‑149 (MQ=255)
                                                              cGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGc                                                                                      >  3:283918/1‑149 (MQ=255)
                                                                  cTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCAtt                                                                                  <  4:220715/149‑1 (MQ=255)
                                                                     gCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaa                                                                               <  4:64717/149‑1 (MQ=255)
                                                                      cTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaaa                                                                              <  1:149260/149‑1 (MQ=255)
                                                                       tGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCaa                                                                                                                           <  4:251411/103‑1 (MQ=255)
                                                                       tGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCaa                                                                                                                           >  3:251411/1‑103 (MQ=255)
                                                                           ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGACCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACt                                                                         >  1:236977/1‑149 (MQ=255)
                                                                              aTTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaaa                                                                              <  2:80573/141‑1 (MQ=255)
                                                                              aTTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGAATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaaa                                                                              >  1:80573/1‑141 (MQ=255)
                                                                               ttCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGc                                                                      >  1:3070/1‑148 (MQ=255)
                                                                                 ccGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCg                                                                     >  1:79183/1‑147 (MQ=255)
                                                                                          gagaAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTc                                                            >  1:6792/1‑147 (MQ=255)
                                                                                           agaAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGa                                                          <  6:11723/148‑1 (MQ=255)
                                                                                                 cGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCggg                                                   <  1:181544/149‑1 (MQ=255)
                                                                                                  ggTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCgggg                                                  <  2:3070/149‑1 (MQ=255)
                                                                                                         gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGcc                                             >  3:6320/1‑147 (MQ=255)
                                                                                                         gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGcc                                             <  2:79183/147‑1 (MQ=255)
                                                                                                         gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                           <  5:219474/149‑1 (MQ=255)
                                                                                                         gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                           <  2:84450/149‑1 (MQ=255)
                                                                                                         gATTCCGAACAGCGTCGCGAAACCCGCATGAAACAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCTAAGAAGCTCGAGGGCGGGGTTGcc                                             <  4:241886/147‑1 (MQ=255)
                                                                                                            tCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGt                                         <  1:263334/148‑1 (MQ=255)
                                                                                                              cGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAc                                      >  5:30092/1‑149 (MQ=255)
                                                                                                              cGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAc                                      <  1:48889/149‑1 (MQ=255)
                                                                                                               gAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAcg                                     >  7:97666/1‑149 (MQ=255)
                                                                                                               gAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGATGACCGCACTGGTTGGACAACTGCAAGGCATGCACCTGCCCAGCGATATCTCGCTGCCATTGAAAAAACTGGCTAAGAAGCTCGAGGCCGGGGTTGCCCAGTCAcg                                     >  2:206720/1‑149 (MQ=255)
                                                                                                                   aGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAg                                                         <  5:263210/125‑1 (MQ=255)
                                                                                                                   aGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAg                                                         >  6:263210/1‑125 (MQ=255)
                                                                                                                    gCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATAGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGc                                >  7:167834/1‑149 (MQ=255)
                                                                                                                     cGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCt                               <  6:30092/149‑1 (MQ=255)
                                                                                                                       tCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTgc                             >  4:136801/1‑149 (MQ=255)
                                                                                                                              aCCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCt                      <  4:6320/149‑1 (MQ=255)
                                                                                                                                  gCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGTCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTg                  <  5:297331/149‑1 (MQ=255)
                                                                                                                                       agtgcggTCAGCGCATCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGc               >  3:33020/6‑147 (MQ=255)
                                                                                                                                         tgcggTCAGCGCATCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTTGGCGAGt           <  4:33020/146‑1 (MQ=255)
                                                                                                                                         ccAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGACGCTGtt                   >  3:249847/1‑141 (MQ=255)
                                                                                                                                            ggTCAGCGCATCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGtt          >  1:60509/1‑147 (MQ=255)
                                                                                                                                                  cGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGCt  >  3:84722/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGCT  >  NC_002947/134249‑134543

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: ‑