Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,408,414 A→C 6.2% V65G (GTG→GGG)  PP_2111 ← aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,408,4140AC 6.2% 168.6 / 2.4 58V65G (GTG→GGG) PP_2111aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (31/23);  new base C (2/2);  total (33/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 = 9.69e-01

GACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGATCAACCGCTGC  >  NC_002947/2408281‑2408557
                                                                                                                                       |                                                                                                                                               
gACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCCCATATGGGCAAACACGGCAGGCGGTACGTc                                                                                                                                 >  2:1444563/1‑150 (MQ=255)
gACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCg                                                                                                                                        >  1:552928/1‑143 (MQ=255)
gACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCg                                                                                                                                        <  2:552928/143‑1 (MQ=255)
      cTGGTTGTTATCGAAAATGGGGGTGACATAGGTGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTg                                                                                                                           <  1:1189293/150‑1 (MQ=255)
         gttgttATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGgcg                                                                                                                        >  1:1211742/1‑150 (MQ=255)
            gttATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCAc                                                                                                                     >  2:475242/1‑150 (MQ=255)
             ttATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCAcc                                                                                                                    >  2:1447382/1‑150 (MQ=255)
                tCGAACATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAgg                                                                                                                 >  2:605874/1‑150 (MQ=255)
                  gAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGtt                                                                                                               <  2:1211742/150‑1 (MQ=255)
                    aaaTGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTtgt                                                                                                             <  1:1444563/150‑1 (MQ=255)
                         ggggTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTgtgcggt                                                                                                        <  1:1447382/150‑1 (MQ=255)
                               aCATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGccc                                                                                                  <  2:367614/150‑1 (MQ=255)
                                 aTAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGt                                                                                                >  2:935354/1‑150 (MQ=255)
                                             ccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTcgcg                                                                                    >  2:209682/1‑150 (MQ=255)
                                             ccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCAAATCTGGGGAAAGACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTcgcg                                                                                    >  1:1335566/1‑150 (MQ=255)
                                              cAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                   <  1:475242/150‑1 (MQ=255)
                                               aGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCa                                                                                                                                            >  2:1139987/1‑92 (MQ=255)
                                               aGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCa                                                                                                                                            <  1:1139987/92‑1 (MQ=255)
                                                 tAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTg                                                                                >  2:426016/1‑150 (MQ=255)
                                                        cGCCGGACTTGCCGCTGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGGGGGTTTCAACGTTCCATGCAAACCTGGCAGGCGGTACGTCGGGGTGGCNCACCAGNTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCAc                                                                         <  1:935354/150‑1 (MQ=255)
                                                               cTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAacac                                                                                                                                                <  2:836148/72‑1 (MQ=255)
                                                               cTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAacac                                                                                                                                                >  1:836148/1‑72 (MQ=255)
                                                                  gCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGAt                                                               >  2:314956/1‑150 (MQ=255)
                                                                      cGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGtt                                                                                                               >  2:1311818/1‑98 (MQ=255)
                                                                      cGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGtt                                                                                                               <  1:1311818/98‑1 (MQ=255)
                                                                           cTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATc                                                      <  1:314956/150‑1 (MQ=255)
                                                                            ttGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCg                                                     >  1:491984/1‑150 (MQ=255)
                                                                               aCGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTg                                                  >  2:855232/1‑150 (MQ=255)
                                                                               aCGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTg                                                  <  1:605874/150‑1 (MQ=255)
                                                                                  aTACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAg                                               <  1:209682/150‑1 (MQ=255)
                                                                                   tACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGTGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGt                                              >  2:1266861/1‑150 (MQ=255)
                                                                                     cccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATATCGATAAAGGCATCGTGGCAGTAg                                            <  2:1335566/150‑1 (MQ=255)
                                                                                      ccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGt                                                                                                >  1:1097473/1‑97 (MQ=255)
                                                                                      ccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGt                                                                                                <  2:1097473/97‑1 (MQ=255)
                                                                                      ccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAgg                                           >  2:500469/1‑150 (MQ=255)
                                                                                      ccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAgg                                           >  2:1473134/1‑150 (MQ=255)
                                                                                        aTCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCa                                                <  1:12139/143‑1 (MQ=255)
                                                                                        aTCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCa                                                >  2:12139/1‑143 (MQ=255)
                                                                                            aCGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACgg                                                                                                                                              >  2:1310564/1‑45 (MQ=255)
                                                                                            aCGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACgg                                                                                                                                              <  1:1310564/45‑1 (MQ=255)
                                                                                               gCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCgggg                                                                                                                             <  1:1603004/59‑1 (MQ=255)
                                                                                               gCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCgggg                                                                                                                             >  2:1603004/1‑59 (MQ=255)
                                                                                                    cccTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTg                             >  1:260528/1‑150 (MQ=255)
                                                                                                       tGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCg                          >  1:399737/1‑150 (MQ=255)
                                                                                                          ttGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGTGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTg                       <  1:1266861/150‑1 (MQ=255)
                                                                                                          gtGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTg                       <  2:139048/149‑1 (MQ=255)
                                                                                                           gccgtgTTTGCCCATATGTGGCAAACCCGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTgg                      >  2:306753/6‑150 (MQ=255)
                                                                                                               ggTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGt                  >  2:989803/1‑150 (MQ=255)
                                                                                                                     ccACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCAc                                                                         <  1:1092402/89‑1 (MQ=255)
                                                                                                                     ccACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCAc                                                                         >  2:1092402/1‑89 (MQ=255)
                                                                                                                     ccACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGAt            >  2:579745/1‑150 (MQ=255)
                                                                                                                      cccaTATGTGGCAAACCCGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCa                                                <  2:847561/112‑1 (MQ=255)
                                                                                                                      cccaTATGTGGCAAACCCGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCa                                                >  1:847561/3‑114 (MQ=255)
                                                                                                                                 gCAAACACGGCAGGCGGTCCGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGATCAACCgctgc  >  2:1381504/1‑150 (MQ=255)
                                                                                                                                  cAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAg                                               <  1:472423/104‑1 (MQ=255)
                                                                                                                                  cAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAg                                               >  2:472423/1‑104 (MQ=255)
                                                                                                                                   aaaCACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTg                       <  1:1346628/127‑1 (MQ=255)
                                                                                                                                   aaaCACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTg                       >  2:1346628/1‑127 (MQ=255)
                                                                                                                                   aaaCACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTg                 >  2:529146/1‑133 (MQ=255)
                                                                                                                                   aaaCACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCATCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTg                 <  1:529146/133‑1 (MQ=255)
                                                                                                                                     acacGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATc                                      <  2:790997/110‑1 (MQ=255)
                                                                                                                                     acacGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATc                                      >  1:790997/1‑110 (MQ=255)
                                                                                                                                      cacGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGga               <  2:219143/132‑1 (MQ=255)
                                                                                                                                      cacGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACGGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGga               >  1:219143/1‑132 (MQ=255)
                                                                                                                                       |                                                                                                                                               
GACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGATCAACCGCTGC  >  NC_002947/2408281‑2408557

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: