Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 964,468 C→A 6.8% P92T (CCG→ACG)  phnE → phosphonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947964,4680CA 6.8% 188.4 / 3.9 59P92T (CCG→ACG) phnEphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base A (2/2);  total (27/32)
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.37e-01

GTCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGC  >  NC_002947/964328‑964610
                                                                                                                                                |                                                                                                                                              
gTCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACggctggc                                                                                                                                        >  3:59307/1‑149 (MQ=255)
   agcagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTg                                                                                                                                     >  7:266429/1‑149 (MQ=255)
    gcagATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGc                                                                                                                                    <  7:211838/149‑1 (MQ=255)
      agATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTg                                                                                                                                  <  8:249644/149‑1 (MQ=255)
          tGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGc                                                                                                                                >  6:57984/1‑147 (MQ=255)
             cAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGcc                                                                                                                           >  1:85326/1‑149 (MQ=255)
              aGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCCAGCAAGCCAGCCg                                                                                                                          <  3:14160/149‑1 (MQ=255)
                 ttCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGc                                                                                                                       >  1:334800/1‑149 (MQ=255)
                   cACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGcc                                                                                                                      <  2:337179/148‑1 (MQ=255)
                    aCCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCgggatc                                                                                                                     <  4:59307/148‑6 (MQ=255)
                       aGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                 <  2:133327/149‑1 (MQ=255)
                       aGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                 >  7:3361/1‑149 (MQ=255)
                       aGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCggga                                                                                                                       >  5:100296/1‑140 (MQ=255)
                       aGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCACGGCTGGCCAGCAAGCCAGCCggga                                                                                                                       <  6:100296/143‑4 (MQ=255)
                          gcCTTCTGGCCGCCGGCACATGATCCCGATTTTCTTGCACTGCTGTATGAATCTACCTTGCAGACCCTGGCCGTTGCGACGATTGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTc                                                                                                                 <  1:115758/149‑1 (MQ=255)
                                tGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGtgc                                                                                                        >  4:352958/1‑149 (MQ=255)
                                      ccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCtt                                                                                                   >  1:3504/1‑148 (MQ=255)
                                      ccgGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTc                                                                                                  <  8:266429/149‑1 (MQ=255)
                                             atgatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                           >  2:151076/1‑149 (MQ=255)
                                             atgatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTACTTCACgcggc                                                                                           <  8:324127/149‑1 (MQ=255)
                                               gatgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGcggcgg                                                                                         >  6:358963/1‑149 (MQ=255)
                                                atgCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcggg                                                                                        <  5:228403/149‑1 (MQ=255)
                                                   ccGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTAGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCggcgggcg                                                                                      <  5:352942/148‑1 (MQ=255)
                                                          tCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGcgggccgggc                                                                              >  8:286521/1‑149 (MQ=255)
                                                            ttGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGctc                                                                            >  1:19975/1‑149 (MQ=255)
                                                              gCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGctc                                                                            <  2:3504/147‑1 (MQ=255)
                                                                aCTGTTGTATGAAGCTACCTTGCAGAGCCTGGCGGTGGCGAC‑‑‑TGCCGGCTTGGCGCGGGCGTTGTTGTTGGCGATCCCGGCTGGTTTGGTGGCCAGCCGTGGGTTGTCGTTGCGTGCTGCTTCACGCGGGGGGCGGGGCGGGCTCTGGt                                                                        <  2:85326/149‑1 (MQ=255)
                                                                aCTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGATGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCCGGCTCTGGt                                                                        >  1:90381/1‑149 (MQ=255)
                                                                    ctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTcgcg                                                                    >  2:139749/1‑149 (MQ=255)
                                                                    ctgTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTcgc                                                                     <  8:3361/148‑1 (MQ=255)
                                                                          gAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGACATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                              >  6:234085/1‑149 (MQ=255)
                                                                           aaGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                              <  7:286521/148‑1 (MQ=255)
                                                                           aaGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                              <  1:151076/148‑1 (MQ=255)
                                                                           aaGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGtt                                                              <  3:352958/148‑1 (MQ=255)
                                                                               tACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                           <  4:198041/115‑1 (MQ=255)
                                                                               tACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACgcggc                                                                                           >  3:198041/1‑115 (MQ=255)
                                                                                  cTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                      >  4:358138/1‑149 (MQ=255)
                                                                                  cTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGc                                                      <  4:189609/149‑1 (MQ=255)
                                                                                     gCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCg                                                    <  3:330743/148‑1 (MQ=255)
                                                                                      cAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGt                                                  >  3:113324/1‑149 (MQ=255)
                                                                                          cccTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGt                                              <  4:113324/149‑1 (MQ=255)
                                                                                          cccTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGt                                              <  5:358963/149‑1 (MQ=255)
                                                                                           ccTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                             <  2:19975/149‑1 (MQ=255)
                                                                                           ccTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTg                                             >  4:115998/1‑149 (MQ=255)
                                                                                                             tGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATTCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTCCGCCTGCCGGTGCGTGGGCTACTGATTTTc                               >  1:43704/1‑149 (MQ=255)
                                                                                                               ccGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTg                            >  4:6192/1‑149 (MQ=255)
                                                                                                                   cATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTa                        <  3:358138/149‑1 (MQ=255)
                                                                                                                    aTGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTa                                        <  3:348565/132‑1 (MQ=255)
                                                                                                                    aTGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTa                                        >  4:348565/1‑132 (MQ=255)
                                                                                                                     tGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGc                      <  2:101565/149‑1 (MQ=255)
                                                                                                                                  gttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATc          <  2:90381/148‑1 (MQ=255)
                                                                                                                                  gttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCg         >  5:373289/1‑149 (MQ=255)
                                                                                                                                  gttgttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCg         >  1:62162/1‑149 (MQ=255)
                                                                                                                                     gttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTg                                                                          >  7:198533/1‑81 (MQ=255)
                                                                                                                                     gttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTg                                                                          <  8:198533/81‑1 (MQ=255)
                                                                                                                                      ttgGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTg     <  5:57984/149‑1 (MQ=255)
                                                                                                                                         gCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTggg   <  2:43704/148‑1 (MQ=255)
                                                                                                                                         gCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGc  >  3:65810/1‑149 (MQ=255)
                                                                                                                                         gCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGc  <  2:334800/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GTCAGCAGATTGCCAGCTTCACCAGCGCCTTCTGGCCGCCGGCACATGATGCCGATTTTCTTGCACTGCTGTATGAAGCTACCTTGCAGACCCTGGCCGTGGCGAC‑‑‑TGCCGGCATGGCGCTGGCGTTGTTGTTGGCGATCCCGGCTGGCTTGCTGGCCAGCCGTGCCTTGTCCTTGCGTGCTGCTTCACGCGGCGGGCGGGCCGGGCTCTGGTCGCGACTGTTGCGCCTGCCGGTGCGTGGGCTACTGATTTTCCTGCGTAGCGTGCCGGAAATCGTCTGGGC  >  NC_002947/964328‑964610

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: