Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,977,580 A→T 6.8% A1017A (GCT→GCA) ‡ cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,977,5800AT 6.8% 170.3 / 3.1 59A1017A (GCT→GCA) ‡cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base A (32/23);  new base T (2/2);  total (34/25)
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

TCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAGG  >  NC_002947/3977437‑3977727
                                                                                                                                                |                                                                                                                                                   
tcgtggCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTg                                                                                                                                                <  2:133931/149‑1 (MQ=255)
  gtcgCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGc                                                                                                                                              <  6:52824/149‑1 (MQ=255)
         ccgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACCCCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                       >  8:148391/1‑149 (MQ=255)
         ccgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                       <  8:183233/149‑1 (MQ=255)
          cgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                       >  1:348334/1‑148 (MQ=255)
          cgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                       <  2:348334/148‑1 (MQ=255)
                  gCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCgg                                                                                                                               >  4:27543/1‑148 (MQ=255)
                      aGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGGACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATc                                                                                                                           >  5:326557/1‑148 (MQ=255)
                      aGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCg                                                                                                                                         >  5:230745/1‑134 (MQ=255)
                      aGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCg                                                                                                                                         <  6:230745/134‑1 (MQ=255)
                         cgccgTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACg                                                                                                                       >  1:290817/1‑149 (MQ=255)
                                    aGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGGCTGGCTGGCgg                                                                                                                                            <  8:290738/117‑2 (MQ=255)
                                    aGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGGCTGGCTGGCgg                                                                                                                                            >  7:290738/1‑116 (MQ=255)
                                      gTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACACGCCAGGCTGGCTGGCgggt                                                                                                                                          >  3:189642/1‑114 (MQ=255)
                                      gTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACACGCCAGGCTGGCTGGCgggt                                                                                                                                          <  4:189642/117‑4 (MQ=255)
                                        aGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAgc                                                                                                        >  6:21820/1‑149 (MQ=255)
                                         gACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAgcg                                                                                                       >  7:121771/1‑149 (MQ=255)
                                         gACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAgc                                                                                                        <  3:340909/148‑1 (MQ=255)
                                           ccTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGc                                                                                                                                   <  7:32782/119‑1 (MQ=255)
                                           ccTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGc                                                                                                                                   >  8:32782/1‑119 (MQ=255)
                                            ctctGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAg                                                                                                         <  4:151761/144‑1 (MQ=255)
                                            ctctGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAg                                                                                                         >  3:151761/1‑144 (MQ=255)
                                                   aGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCcgctc                                                                                             >  2:284419/1‑149 (MQ=255)
                                                    ggTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCcgctcg                                                                                            <  3:339906/149‑1 (MQ=255)
                                                          cgcgGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTgcgc                                                                                      <  7:148391/149‑1 (MQ=255)
                                                           gcgGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGca                                                                                     <  3:27543/149‑1 (MQ=255)
                                                                 gtgCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACggg                                                                               >  1:187613/1‑149 (MQ=255)
                                                                  tgCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggc                                                                              >  7:103245/1‑149 (MQ=255)
                                                                   gCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcg                                                                             >  3:259975/1‑149 (MQ=255)
                                                                    ccACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcgg                                                                            >  5:154200/1‑149 (MQ=255)
                                                                      acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTACGCGCCCGCTTTGGCACCGAACACCCGCCACCCAGCCTGCCGTGCGGCCTGCTCGGCATCCACGCCTTGACCCTGCACCGCGGCCCGCTCGCTGCCGACACGcgcggc                                                                           >  8:233884/1‑148 (MQ=255)
                                                                      acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcggc                                                                           >  6:310916/1‑148 (MQ=255)
                                                                      acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcggc                                                                           >  2:141069/1‑148 (MQ=255)
                                                                      acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGcc                                                                          <  2:226324/149‑1 (MQ=255)
                                                                        aaaCGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAg                                                                        >  8:234605/1‑149 (MQ=255)
                                                                          aCGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAggg                                                                      >  8:70082/1‑149 (MQ=255)
                                                                            gCCCGTCGATGGCGTTCTGCACCCCGGCACCAGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCCTGCGGCCTGCTCGCCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAgggg                                                                     >  1:232786/1‑149 (MQ=255)
                                                                              ccGTCGATGGCGTTCTGCACCCCGGCACCGTTAGCGCCCGGTTTGGCACCGAACTCCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGtt                                                                  <  3:302458/149‑1 (MQ=255)
                                                                                 tCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGATCACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCGACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAg                                                               <  2:232786/149‑1 (MQ=255)
                                                                                 tCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAg                                                               <  7:70082/149‑1 (MQ=255)
                                                                                 tCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAg                                                               <  6:154200/149‑1 (MQ=255)
                                                                                   gATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAgg                                                                       >  1:111326/1‑139 (MQ=255)
                                                                                   gATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAgg                                                                       <  2:111326/139‑1 (MQ=255)
                                                                                     tGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGTCCCGCTCGCTGCGCACACGggcggc                                                                           <  3:155080/133‑1 (MQ=255)
                                                                                     tGGCGTTCTGCACACCGGCACAGTAAGCGCCCGGTTTGGCAAAGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCAACCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGcgcggc                                                                           >  4:155080/1‑133 (MQ=255)
                                                                                       gCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTc                                                         <  5:21820/149‑1 (MQ=255)
                                                                                              gCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCg                                                  >  3:176807/1‑149 (MQ=255)
                                                                                                      gCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAgcgc                                          <  5:310916/149‑1 (MQ=255)
                                                                                                         ccGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGc                                       >  6:218766/1‑149 (MQ=255)
                                                                                                                      ggTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCAc                           >  1:258257/1‑149 (MQ=255)
                                                                                                                       gTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCAcg                          >  2:93462/1‑149 (MQ=255)
                                                                                                                        tttGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCAcgg                         >  6:161967/1‑149 (MQ=255)
                                                                                                                              aCCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGTGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGt                   >  1:229461/1‑149 (MQ=255)
                                                                                                                              aCCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGt                   >  7:82201/1‑149 (MQ=255)
                                                                                                                               ccGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTc                  <  4:176807/149‑1 (MQ=255)
                                                                                                                                   acacCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAAc              <  6:326557/149‑1 (MQ=255)
                                                                                                                                     acCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAg            >  7:302322/1‑149 (MQ=255)
                                                                                                                                      cccGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAgg           <  1:141069/149‑1 (MQ=255)
                                                                                                                                       ccGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGc          >  5:339674/1‑149 (MQ=255)
                                                                                                                                         gccagccaGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTgc                                                                                        >  4:116827/1‑69 (MQ=255)
                                                                                                                                         gccagccaGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTgc                                                                                        <  3:116827/69‑1 (MQ=255)
                                                                                                                                         gcacgccaggcTGGCTGGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAggg                                                                      <  4:315886/83‑1 (MQ=255)
                                                                                                                                         gcacgccaggcTGGCTGGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAggg                                                                      >  3:315886/5‑87 (MQ=255)
                                                                                                                                              ccagccTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAg   <  2:258257/149‑1 (MQ=255)
                                                                                                                                               cagccTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAgg  >  5:87138/1‑149 (MQ=255)
                                                                                                                                               cagccTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAgg  <  6:87138/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACCGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAGG  >  NC_002947/3977437‑3977727

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: