Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,141,143 A→G L157P (CTG→CCG)  asmA ← suppressor of OmpF assembly mutants, putative outer membrane protein assembly factor, inner membrane‑anchored periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,141,1430AG94.4% 241.7 / 2.4 72L157P (CTG→CCG) asmAsuppressor of OmpF assembly mutants, putative outer membrane protein assembly factor, inner membrane‑anchored periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (31/37);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCAGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTGCGGCGTCAGT  >  NC_000913/2141008‑2141277
                                                                                                                                       |                                                                                                                                      
gATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCa                                                                                                                                     <  2:68045/139‑1 (MQ=255)
  tatCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCAACTGCTCGTCATCTTCATGCTGGAACACCGGCACg                                                                                                                                   >  2:661116/1‑139 (MQ=255)
   atCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGc                                                                                                                                  >  1:1128460/1‑139 (MQ=255)
   atCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGc                                                                                                                                  <  1:874994/139‑1 (MQ=255)
    tCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCt                                                                                                                                 >  1:590813/1‑139 (MQ=255)
       tGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTc                                                                                                                              >  1:1230363/1‑139 (MQ=255)
        gagaTCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTcc                                                                                                                             >  2:1607206/1‑139 (MQ=255)
          gaTCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           <  2:1603066/139‑1 (MQ=255)
          gaTCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           <  2:1230363/139‑1 (MQ=255)
          gaTCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           <  1:677859/139‑1 (MQ=255)
          gaTCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           <  1:1607206/139‑1 (MQ=255)
             cgcgCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCAc                                                                                                                        <  2:847338/139‑1 (MQ=255)
               cgcTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCt                                                                                                                      >  2:996387/1‑139 (MQ=255)
               cgcTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCt                                                                                                                      <  2:590813/139‑1 (MQ=255)
               cgcTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCt                                                                                                                      <  2:1450521/139‑1 (MQ=255)
               cgcTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCt                                                                                                                      <  1:996387/139‑1 (MQ=255)
                gcTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCtt                                                                                                                     <  1:113649/139‑1 (MQ=255)
                gcTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCtt                                                                                                                     <  2:1128460/139‑1 (MQ=255)
                         gATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGa                                                                                                            <  1:1120742/139‑1 (MQ=255)
                             aaCCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGatata                                                                                                        <  2:973654/139‑1 (MQ=255)
                              aCCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGatatat                                                                                                       >  2:1628865/1‑139 (MQ=255)
                                   cccgGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAAc                                                                                                  >  2:212745/1‑139 (MQ=255)
                                     cgGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGa                                                                                                <  2:823838/139‑1 (MQ=255)
                                              cAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCcgcg                                                                                       <  2:842518/139‑1 (MQ=255)
                                              cAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCcgcg                                                                                       <  1:412929/139‑1 (MQ=255)
                                              cAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCAGCACg                                                                                                                                   <  1:1177558/95‑1 (MQ=255)
                                              cAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCAGCACg                                                                                                                                   >  2:1177558/1‑95 (MQ=255)
                                                aaTGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGAt                                                                                     <  2:1490852/139‑1 (MQ=255)
                                                      ccACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTg                                                                               >  2:1151185/1‑139 (MQ=255)
                                                        aCGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGAc                                                                             <  2:1453203/139‑1 (MQ=255)
                                                         cGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACa                                                                            >  2:97876/1‑139 (MQ=255)
                                                           aTGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGa                                                                          <  2:1165079/139‑1 (MQ=255)
                                                                  ggATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCAGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCa                                                                   <  1:1339366/139‑1 (MQ=255)
                                                                    aTCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           >  1:296317/1‑81 (MQ=255)
                                                                    aTCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           <  2:296317/81‑1 (MQ=255)
                                                                      cTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCcgcg                                                                                       >  1:748043/1‑115 (MQ=255)
                                                                      cTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCcgcg                                                                                       <  2:748043/115‑1 (MQ=255)
                                                                        tGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           >  1:131682/1‑77 (MQ=255)
                                                                        tGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           <  2:131682/77‑1 (MQ=255)
                                                                        tGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTa                                                                                                                    <  1:214229/84‑1 (MQ=255)
                                                                        tGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTa                                                                                                                    >  2:214229/1‑84 (MQ=255)
                                                                        tGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGa                                                                                                                 >  2:27672/1‑87 (MQ=255)
                                                                        tGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGa                                                                                                                 <  1:27672/87‑1 (MQ=255)
                                                                              aTTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGa                                                                                                >  2:1259400/1‑98 (MQ=255)
                                                                              aTTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGa                                                                                                <  1:1259400/98‑1 (MQ=255)
                                                                               tttGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           >  1:274379/1‑70 (MQ=255)
                                                                               tttGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGc                                                                                                                           <  2:274379/70‑1 (MQ=255)
                                                                                 tGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCgg                                                                     >  2:520748/1‑122 (MQ=255)
                                                                                 tGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCgg                                                                     <  1:520748/122‑1 (MQ=255)
                                                                                           tatTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCgg                                          >  1:1236306/1‑139 (MQ=255)
                                                                                              tGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCtt                                                                                                                     <  1:86733/61‑1 (MQ=255)
                                                                                              tGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCtt                                                                                                                     >  2:86733/1‑61 (MQ=255)
                                                                                                    ttGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTc                                 >  2:558918/1‑139 (MQ=255)
                                                                                                     tGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCcgcg                                                                                       <  2:273874/84‑1 (MQ=255)
                                                                                                     tGTCACCTGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCcgcg                                                                                       >  1:273874/1‑84 (MQ=255)
                                                                                                            tGCTCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCa                         <  2:494094/139‑1 (MQ=255)
                                                                                                               tCGTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCg                      >  2:341863/1‑139 (MQ=255)
                                                                                                                 gTCATCTTCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCt                    <  2:624100/139‑1 (MQ=255)
                                                                                                                        tCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGt                                                                                  >  2:1607663/1‑70 (MQ=255)
                                                                                                                        tCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGt                                                                                  <  1:1607663/70‑1 (MQ=255)
                                                                                                                        tCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATc                                                                        <  1:155769/80‑1 (MQ=255)
                                                                                                                        tCATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATc                                                                        >  2:155769/1‑80 (MQ=255)
                                                                                                                         cATGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTg            >  1:1081811/1‑139 (MQ=255)
                                                                                                                          aTGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTgg                                                                                                             <  1:1308438/41‑1 (MQ=255)
                                                                                                                          aTGCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTgg                                                                                                             >  2:1308438/1‑41 (MQ=255)
                                                                                                                            gCTGGAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTgcgg         >  2:553281/1‑139 (MQ=255)
                                                                                                                                gAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTg                                                                               <  2:1584360/65‑1 (MQ=255)
                                                                                                                                gAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTg                                                                               >  1:1584360/1‑65 (MQ=255)
                                                                                                                                gAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTGCGGCGTc     >  1:1417479/1‑139 (MQ=255)
                                                                                                                                gAACACCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTGCGGCGTc     <  1:341863/139‑1 (MQ=255)
                                                                                                                                   cacCGGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTGCGGCGTCAGt  >  1:164330/1‑139 (MQ=255)
                                                                                                                                   cacCAGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTGCGGCGTCAGt  >  2:1226444/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAATGAACCACGATGTTGGGGATCTTGTTCCATTTGCAGGCGAATATTGCGGATTGTCACCTGCTCGTCATCTTCATGCTGGAACACCAGCACGCTGTCCGCCACCTTAAGACTGGATATATCAAACGACCATCCGCGATCGTCTGACAGATCCGGCAAGGTATTGTCGCGCGGTGCAACCGGAGCGTCTTCACTGCGCACCGCTTCCGTCTGCGGCGTCAGT  >  NC_000913/2141008‑2141277

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