Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,246,756 T→C N5S (AAC→AGC)  cirA ← colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,246,7560TC98.3% 211.9 / NA 60N5S (AAC→AGC) cirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (35/24);  minor base G (1/0);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAAC  >  NC_000913/2246608‑2246900
                                                                                                                                                    |                                                                                                                                                
aaTGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACCAACACAGCCCGACCCGTACGAAAGGGCTc                                                                                                                                                >  2:564654/1‑151 (MQ=255)
 aTGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCa                                                                                                                                               <  1:637170/151‑1 (MQ=255)
 aTGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCa                                                                                                                                               <  2:405097/151‑1 (MQ=255)
  tGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCa                                                                                                                                               >  1:619813/1‑150 (MQ=255)
   gACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCaa                                                                                                                                              <  1:657708/150‑1 (MQ=255)
        tGATACTGGCAGGTGCATCTTTAAGGTTTTGTTCCACGGAAAATTAAGTGACAACCATCGTTTCGCCATCATCATTGTCCGCTAACACAAGGCATGCACAAGAAATAGCGGACAAACACAGCCCCACCCGCACGAAAAGGGTCaaac                                                                                                                                            >  1:500938/1‑145 (MQ=255)
        tGATACTGGCAGGTGCATCTTTAAGAATTTTGTCCACTGAAGATGCCAGGCCAACCAACCTTTCGCCATCATCATCGACCCCGAACACAAGCCATGCACAAGAAATAACGCAGAAACACAGCCCGATCCGTCCGCAACGGCTCGACCtca                                                                                                                                         >  1:26352/1‑148 (MQ=255)
         gATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTaa                                                                                                                                         <  1:1164681/149‑1 (MQ=255)
            aCTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACAtt                                                                                                                                    <  2:118877/151‑1 (MQ=255)
            aCTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACAt                                                                                                                                     <  1:715005/150‑1 (MQ=255)
            aCTGGCAGGTGCAACTTTAAACTTTTGTTTCACGGGAGATGCAGTGACAACCATCGTTTTGCCAACATCATCGACCGCTAACACAGGCCATTCACAAGAAAAAGCGGACAAACACAACCCGACCCCTACGAAAAGGCCCAACCTAAACAtt                                                                                                                                    >  1:824639/1‑151 (MQ=255)
               ggCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTc                                                                                                                                   <  2:802300/149‑1 (MQ=255)
                 cAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCat                                                                                                                                >  2:310838/1‑150 (MQ=255)
                 cAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCat                                                                                                                                >  1:1254299/1‑150 (MQ=255)
                  aGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCatat                                                                                                                              <  1:262201/151‑1 (MQ=255)
                   ggTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATAtc                                                                                                                             <  2:876247/151‑1 (MQ=255)
                      gCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATAtctc                                                                                                                           <  2:595408/150‑1 (MQ=255)
                          cTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCAGCACGAAAGGGCTCAACCTAAACATTCCATATCTCAATg                                                                                                                       >  1:564899/1‑150 (MQ=255)
                                    ttGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACg                                                                                                            >  2:1208314/1‑151 (MQ=255)
                                      gTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGa                                                                                                           >  1:1110959/1‑150 (MQ=255)
                                       ttCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaa                                                                                                          >  1:267220/1‑150 (MQ=255)
                                        tCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaa                                                                                                         <  2:1228163/150‑1 (MQ=255)
                                           aCGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGGCCGCTAACACAGGCCATGCACAAGCAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACAGTCCATATCTCCATGAGGTAACTACg                                                                                                            <  2:564929/144‑1 (MQ=255)
                                           aCGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataa                                                                                                     <  2:330016/151‑1 (MQ=255)
                                            cGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataa                                                                                                     >  1:1060544/1‑150 (MQ=255)
                                             ggAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaa                                                                                                    >  1:503301/1‑150 (MQ=255)
                                              gAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaaat                                                                                                  >  1:297087/1‑151 (MQ=255)
                                                 gATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATggg                                                                                               <  1:1259459/151‑1 (MQ=255)
                                                  aTGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGt                                                                                              >  2:25236/1‑151 (MQ=255)
                                                      aGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATc                                                                                          <  1:25241/151‑1 (MQ=255)
                                                        tGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGc                                                                                        >  2:639490/1‑151 (MQ=255)
                                                         gACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCt                                                                                       >  2:635888/1‑151 (MQ=255)
                                                             aCCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTcaca                                                                                   <  1:919859/151‑1 (MQ=255)
                                                                 tCGTTTCGTCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtctt                                                                               >  2:332707/1‑151 (MQ=255)
                                                                 acGTTTAGCCATCATCATCGACCGCTAACACAAGTCATGCACAAGAAATAGCGGACAAACCCAGCCCGAGCCGTAAGAAAGAGCACAACCTAAGCATTCAATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtc                                                                                 >  1:245225/2‑149 (MQ=255)
                                                                           atcatcatcGACCGCTAACACAGGCCATGCACAAGAACTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTcc                                                                     >  1:636093/1‑151 (MQ=255)
                                                                             catcatcGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                                   >  1:78211/1‑151 (MQ=255)
                                                                                catcGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCtt                                                                 <  1:287253/150‑1 (MQ=255)
                                                                                   cGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCg                                                              <  2:888922/150‑1 (MQ=255)
                                                                                      ccGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACg                                                          >  2:194133/1‑151 (MQ=255)
                                                                                       cGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACgg                                                         >  1:458693/1‑151 (MQ=255)
                                                                                             cacaGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCAt                                                   >  1:333126/1‑151 (MQ=255)
                                                                                              acaGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATc                                                  >  1:1022594/1‑151 (MQ=255)
                                                                                                aGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgc                                                <  2:722167/151‑1 (MQ=255)
                                                                                                aGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgc                                                <  2:1242887/151‑1 (MQ=255)
                                                                                                 ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgcg                                               >  1:406595/1‑151 (MQ=255)
                                                                                                      tGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                           >  2:1240543/1‑150 (MQ=255)
                                                                                                       gCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                          <  1:882549/150‑1 (MQ=255)
                                                                                                               aaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttt                                 >  1:1265457/1‑151 (MQ=255)
                                                                                                               aaTAGCGCACAAACACAGCCCGACCCGTACGAAAGGGCCCCACCTAAACATTCCATATATCCATGAGGTAACTACGACACTAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCACCATCGCGGTGGACTTAttttt                                 >  1:1173914/1‑151 (MQ=255)
                                                                                                                             acaGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatc                    >  2:1273052/1‑150 (MQ=255)
                                                                                                                               aGCCCGACCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcatc                 >  1:94915/1‑151 (MQ=255)
                                                                                                                                     aCCCGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGa            >  1:553744/1‑150 (MQ=255)
                                                                                                                                      cccGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa          <  2:298686/151‑1 (MQ=255)
                                                                                                                                      cccGTACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaa           >  2:326130/1‑150 (MQ=255)
                                                                                                                                          tACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                           <  2:220540/130‑1 (MQ=255)
                                                                                                                                          tACGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                           >  1:220536/1‑130 (MQ=255)
                                                                                                                                           aCGAAAGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCgg     <  2:739053/151‑1 (MQ=255)
                                                                                                                                              aaaGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAAc  >  2:563247/1‑151 (MQ=255)
                                                                                                                                                aGGGCTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAAc  >  2:1125268/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAAC  >  NC_000913/2246608‑2246900

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