Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,826,825 G→C 12.5% A80P (GCA→CCA)  PP_2479 → cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,826,8250GC12.5% 154.0 / 13.2 56A80P (GCA→CCA) PP_2479cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base C (5/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GGGCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGC  >  NC_002947/2826677‑2826968
                                                                                                                                                    |                                                                                                                                               
gggCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTg                                                                                                                                                 >  8:72189/1‑149 (MQ=255)
gggCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTg                                                                                                                                                 >  1:10366/1‑149 (MQ=255)
gggCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTg                                                                                                                                                 <  6:268305/149‑1 (MQ=255)
 ggCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGc                                                                                                                                                <  4:82179/149‑1 (MQ=255)
   cGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCAt                                                                                                                                              <  2:10366/149‑1 (MQ=255)
       tACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                          <  3:260550/149‑1 (MQ=255)
       tACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                          <  6:261092/149‑1 (MQ=255)
           tGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGAt                                                                                                                                      >  1:68403/1‑149 (MQ=255)
               gcgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGTCATCGGTCCATGCACCGATGCCg                                                                                                                                  <  6:249049/149‑1 (MQ=255)
                cgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGa                                                                                                                                 <  2:99962/149‑1 (MQ=255)
                 gcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAg                                                                                                                                >  6:247408/1‑149 (MQ=255)
                       ccAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTc                                                                                                                          <  2:68403/149‑1 (MQ=255)
                              tgtgGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTgcgc                                                                                                                   >  2:125/1‑149 (MQ=255)
                                        cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGAc                                                                                                                                   >  7:165257/1‑123 (MQ=255)
                                        cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGAc                                                                                                                                   <  8:165257/123‑1 (MQ=255)
                                        cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCgg                                                                                                         >  5:61109/1‑149 (MQ=255)
                                          tACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGAt                                                                                                                                      <  4:261770/118‑1 (MQ=255)
                                              gTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTgg                                                                                                   <  6:61109/149‑1 (MQ=255)
                                                  cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGcc                                                                                                 >  7:20637/1‑147 (MQ=255)
                                                  cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGcc                                                                                                 <  8:20637/147‑1 (MQ=255)
                                                    ggCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGa                                                                                             >  4:209048/1‑149 (MQ=255)
                                                       ggTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACgg                                                                                          >  6:30877/1‑149 (MQ=255)
                                                             cccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCa                                                                                    >  3:155096/1‑149 (MQ=255)
                                                              ccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCac                                                                                   <  4:155096/149‑1 (MQ=255)
                                                                  cGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCacta                                                                                 <  1:125/147‑1 (MQ=255)
                                                                       gtggtgTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCg                                                                          <  5:247408/149‑1 (MQ=255)
                                                                        tggtgTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCgg                                                                         >  6:170637/1‑149 (MQ=255)
                                                                           tgtTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGt                                                                      >  3:10346/1‑149 (MQ=255)
                                                                             ttGAGTTCAAGCTGCCGTTCGGCTAGCTGTTTTCGCCCAAGATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCTATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGttt                                                                    <  3:216872/149‑1 (MQ=255)
                                                                               gAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                  >  2:38452/1‑149 (MQ=255)
                                                                                 gTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTccc                                                                 <  2:86839/148‑1 (MQ=255)
                                                                                  ttCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                  >  1:186710/1‑146 (MQ=255)
                                                                                  ttCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                  <  2:186710/146‑1 (MQ=255)
                                                                                  ttCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCta                                                               <  1:38452/149‑1 (MQ=255)
                                                                                            cGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTAc                                                     <  5:30877/149‑1 (MQ=255)
                                                                                                  gCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCac                                                                                   >  4:15550/1‑113 (MQ=255)
                                                                                                          ttttCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgcg                                       >  1:57985/1‑149 (MQ=255)
                                                                                                              cGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGa                                    <  1:134183/148‑1 (MQ=255)
                                                                                                                cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                  <  1:268062/116‑1 (MQ=255)
                                                                                                                cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                  >  2:268062/1‑116 (MQ=255)
                                                                                                                 ccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGa                                 >  4:110312/1‑148 (MQ=255)
                                                                                                                    aCATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATc                             >  5:71962/1‑149 (MQ=255)
                                                                                                                       ttACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTc                          >  5:85915/1‑149 (MQ=255)
                                                                                                                       ttACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGCCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTc                          <  5:170637/149‑1 (MQ=255)
                                                                                                                           cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgc                                        >  3:71571/1‑131 (MQ=255)
                                                                                                                           cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgc                                        <  4:71571/131‑1 (MQ=255)
                                                                                                                           cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGa                      >  8:35659/1‑149 (MQ=255)
                                                                                                                               gACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCaa                  <  3:209048/149‑1 (MQ=255)
                                                                                                                                 cccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAgg                <  6:85915/149‑1 (MQ=255)
                                                                                                                                  ccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTGCTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAggg               <  2:56525/149‑1 (MQ=255)
                                                                                                                                       aaCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGGCATCCTCGCGATCAAGGGTTAtc          <  8:59088/149‑1 (MQ=255)
                                                                                                                                       aaCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtc          <  6:71962/149‑1 (MQ=255)
                                                                                                                                        ctcgtcATCGGTCCATGCACCGATGCCGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtct         >  2:95823/6‑149 (MQ=255)
                                                                                                                                        ctcgtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtct         >  5:196818/6‑149 (MQ=255)
                                                                                                                                         tcgtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGAt                                   >  1:19607/5‑122 (MQ=255)
                                                                                                                                         tcgtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGAt                                   <  2:19607/118‑1 (MQ=255)
                                                                                                                                           gtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGa      >  5:12718/3‑149 (MQ=255)
                                                                                                                                           gtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGa      >  5:51916/3‑149 (MQ=255)
                                                                                                                                               acGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATAGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTCTACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGc  >  6:129555/2‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGGCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGC  >  NC_002947/2826677‑2826968

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: