Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,081,564 C→T 6.6% E7E (GAG→GAA PP_5333 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,081,5640CT 6.6% 178.9 / 2.8 60E7E (GAG→GAAPP_5333hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base T (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTCGATTTCAGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTGGCCGCCATAGG  >  NC_002947/6081417‑6081712
                                                                                                                                                   |                                                                                                                                                    
gTCGATTTCAGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACt                                                                                                                                                     <  5:230926/149‑1 (MQ=255)
   gATTTCAGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAATTcat                                                                                                                                                  >  3:184468/1‑147 (MQ=255)
    aTTTCAGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCgaa                                                                                                                                                 >  5:53498/1‑149 (MQ=255)
     tttCAGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAATTcattc                                                                                                                                                >  7:166430/1‑145 (MQ=255)
     tttCAGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCgaa                                                                                                                                                 <  2:125473/148‑1 (MQ=255)
         aGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCgaatgaa                                                                                                                                             <  8:276740/148‑1 (MQ=255)
                tGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTGGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGa                                                                                                                                      <  6:303930/148‑1 (MQ=255)
                tGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGa                                                                                                                                      <  3:317915/148‑1 (MQ=255)
                 gTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAATTCATTCGAGTTGCCGAg                                                                                                                                     <  8:166430/148‑1 (MQ=255)
                    cGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAATTCATTCGAGTTGCCGAGCAt                                                                                                                                  <  4:184468/148‑1 (MQ=255)
                     gAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCAtt                                                                                                                                 <  7:119574/148‑1 (MQ=255)
                     gAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTc                                                                                                                                <  5:323195/149‑1 (MQ=255)
                      aGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTc                                                                                                                                >  7:75562/1‑148 (MQ=255)
                         gcCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAAc                                                                                                                            <  6:53498/149‑1 (MQ=255)
                         gacTTCTGGTCCGTCCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGTTACTTGCCAACATCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCTAATGAATTGCAGAGCATTCTAAc                                                                                                                            <  6:296035/147‑1 (MQ=255)
                          ccTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAAc                                                                                                                            >  1:337113/1‑148 (MQ=255)
                           cTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGCGCATTCTAAc                                                                                                                            >  7:269188/1‑147 (MQ=255)
                           cTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAAc                                                                                                                            >  1:317903/1‑147 (MQ=255)
                           cTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAAc                                                                                                                            >  5:184133/1‑147 (MQ=255)
                           cTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCa                                                                                                                          <  5:113627/149‑1 (MQ=255)
                           cTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCa                                                                                                                          <  7:225818/149‑1 (MQ=255)
                                 gTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCATGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATc                                                                                                                    <  4:90855/149‑1 (MQ=255)
                                   ccGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGAAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAg                                                                                                                  >  3:315003/1‑149 (MQ=255)
                                     gACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCg                                                                                                                >  5:44861/1‑149 (MQ=255)
                                     gACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGAAACTCGAAGGAATTGACGCGCATTCTAACCAGGAATCAGCg                                                                                                                >  8:137987/1‑149 (MQ=255)
                                      aCCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGt                                                                                                               >  8:7495/1‑149 (MQ=255)
                                         tGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTc                                                                                                            >  8:258401/1‑149 (MQ=255)
                                           ggCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTc                                                                                                            <  2:66333/147‑1 (MQ=255)
                                           ggCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCt                                                                                                           <  6:243503/148‑1 (MQ=255)
                                           ggCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCGCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCtt                                                                                                          <  5:156290/149‑1 (MQ=255)
                                             cATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTtga                                                                                                        <  7:7495/149‑1 (MQ=255)
                                                  gCTCCAGCCAGTCGGTTTCGTGCAGTTCGGCATGGGCCAGTGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGTCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAAt                                                                                                   <  7:235536/149‑1 (MQ=255)
                                                  gCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAAt                                                                                                   >  8:3784/1‑149 (MQ=255)
                                                       agccagTCGGGTTCGGGAAGTTCGGCATGGGCCAGGGTACTGGCCAACAGGAATGGGACGAGTAAAAAGTCACTCATGAAAATGCGCGTCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTa                                                                                              <  7:137987/149‑1 (MQ=255)
                                                         ccagTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTat                                                                                                                      >  8:204353/1‑122 (MQ=255)
                                                         ccagTCGGTTTCGGGCAGTTCGGCATGGGCACGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGAAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCg                                                                                                >  6:128085/1‑145 (MQ=255)
                                                          cagTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCg                                                                                                >  3:172799/1‑144 (MQ=255)
                                                               ggTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGAAGATCACca                                                                                      <  6:44861/149‑1 (MQ=255)
                                                                gTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGGAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACcag                                                                                     >  7:102233/1‑149 (MQ=255)
                                                                 tttCGGGCAGTTCGGCATGGGCCAGGGTACTTGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATTAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACca                                                                                      <  8:256936/147‑1 (MQ=255)
                                                                 tttCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACcagc                                                                                    <  8:75562/149‑1 (MQ=255)
                                                                   tCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACcagcag                                                                                  >  4:137514/1‑149 (MQ=255)
                                                                              ggCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCAc                                                                       <  2:337113/149‑1 (MQ=255)
                                                                                          ggTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTtg                                                                                                         <  3:269616/103‑1 (MQ=255)
                                                                                          ggTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTtg                                                                                                         >  4:269616/1‑103 (MQ=255)
                                                                                            tACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCAcc                                                          >  7:47272/1‑148 (MQ=255)
                                                                                                 gCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTa                                                                                              <  6:241964/107‑1 (MQ=255)
                                                                                                 gCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCGAACCAGTAATCAGCGTTTCTTGATGAATCCGTa                                                                                              >  5:241964/1‑107 (MQ=255)
                                                                                                          aaGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCccgccg                                            >  5:60209/1‑148 (MQ=255)
                                                                                                            ggggACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCt                                         <  8:47272/149‑1 (MQ=255)
                                                                                                              ggACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGAc                                                                   <  4:162372/121‑1 (MQ=255)
                                                                                                                  gAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAg                                   >  6:141709/1‑149 (MQ=255)
                                                                                                                          aGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGtat                           <  8:102233/149‑1 (MQ=255)
                                                                                                                                gCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATAccc                     >  3:60902/1‑149 (MQ=255)
                                                                                                                                       aaTGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGta                            >  6:139380/1‑135 (MQ=255)
                                                                                                                                        aTGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTg             <  7:3784/149‑1 (MQ=255)
                                                                                                                                            tCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTGgccg         >  4:110404/1‑149 (MQ=255)
                                                                                                                                            tCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTGgccg         >  3:137901/1‑149 (MQ=255)
                                                                                                                                              ggCAACTCGAATGAATTGCCGAGCATTCTATCCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTGgccgc        >  1:27522/1‑148 (MQ=255)
                                                                                                                                                cAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTGGCCGCCAt     <  4:172799/149‑1 (MQ=255)
                                                                                                                                                   cTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTGGCCGCCATAgg  <  4:60902/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCGATTTCAGGCATTTGTTCGAGCGCCTTCTGGTCCGACCTGGGCATCAGCTCCAGCCAGTCGGTTTCGGGCAGTTCGGCATGGGCCAGGGTACTGGCCAACAGCAAGGGGACGAGGAAAAAGGCACGCATGAAAATGCTCGGCAACTCGAATGAATTGCCGAGCATTCTAACCAGTAATCAGCGTTTCTTGATGAATCCGTAGATCACCAGCAGGACCACGGCACCGACCAGCGCACCAAAGAAGCCCGCCGCCTGCCCAGCCTGGTATATACCCAGCGCCTGGCCGCCATAGG  >  NC_002947/6081417‑6081712

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