Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,192,7690GT41.0% 148.6 / 190.6 190intergenic (+228/+192)b3051/b3052putative membrane protein/putative kinase
Reads supporting (aligned to +/- strand):  ref base G (46/66);  new base T (39/39);  total (85/105)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 193 total aligned reads displayed.

TGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAA  >  NC_000913/3192694‑3192841
                                                                           |                                                                        
tGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAg                                                                          <  2:4410162/76‑1 (MQ=255)
tGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAg                                                                          >  1:988667/1‑76 (MQ=255)
gggCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAg                                                                          <  2:1298540/75‑1 (MQ=255)
 ggCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGa                                                                         <  1:5425060/76‑1 (MQ=255)
    gcTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAgg                                                                      <  2:849748/76‑1 (MQ=255)
    gcTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAgg                                                                      <  1:1474670/76‑1 (MQ=255)
      tcatcaACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGaa                                                                    >  2:599994/1‑76 (MQ=255)
      tcatcaACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGaa                                                                    <  1:3927559/76‑1 (MQ=255)
      tcatcaACAATCGGGGAGGCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGaa                                                                    <  2:5403955/76‑1 (MQ=255)
      tcatcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAg                                                                  >  2:2155785/1‑76 (MQ=17)
       catcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAgg                                                                 <  1:4457299/76‑1 (MQ=12)
       catcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCCAAGAAGGGGCCTTGTCTAAGGAAAgg                                                                 >  2:169012/1‑76 (MQ=12)
        atcaCCAATCGGGGTGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAg                                                                  <  2:1822710/76‑1 (MQ=255)
         tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                                               >  1:3806909/1‑76 (MQ=9)
           aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGGTa                                                             >  1:1041603/1‑76 (MQ=9)
            acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGCTat                                                            >  2:2276584/1‑76 (MQ=9)
            acaaGCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGtt                                                              <  2:155713/76‑1 (MQ=255)
             caaTGGGGGAGTCGGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTa                                                             <  2:4203319/76‑1 (MQ=255)
              aaTCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTat                                                            <  1:4140852/76‑1 (MQ=255)
              aaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatga                                                          <  1:5607377/76‑1 (MQ=9)
               aTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatgat                                                         <  1:2575299/76‑1 (MQ=9)
                  ggggAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTatgatg                                                        >  2:1768825/1‑76 (MQ=255)
                  ggggAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTatgatg                                                        >  1:4063453/1‑76 (MQ=255)
                      aGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGc                                                    >  2:2534255/1‑76 (MQ=255)
                       ggcagcaaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                                                  <  2:2030187/74‑1 (MQ=14)
                       ggcagcaaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                                                  <  2:880759/74‑1 (MQ=14)
                       ggcagcaaGGGGCTGAAACGGGAAAGCCCCGCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                                                  <  1:305657/74‑1 (MQ=14)
                       gTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCt                                                   <  1:5833305/76‑1 (MQ=255)
                       gTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCt                                                   >  1:1794060/1‑76 (MQ=255)
                        tCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTc                                                  <  1:660057/76‑1 (MQ=255)
                        tCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTc                                                  <  1:4976629/76‑1 (MQ=255)
                        gcagcaaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                                                 <  2:1407759/75‑1 (MQ=14)
                         gagcaaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGt                                                <  2:4372160/75‑1 (MQ=14)
                          agcagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGt                                                >  1:2756823/1‑76 (MQ=255)
                          agcaaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtc                                               <  2:665396/76‑1 (MQ=17)
                           gcagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGtc                                               <  2:3179802/76‑1 (MQ=255)
                           gcagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGATCGtc                                               >  2:147023/1‑76 (MQ=255)
                           gcaaGGGGCTGAACCGGGAATGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGCAGCTCGtca                                              <  1:3482142/76‑1 (MQ=12)
                            caaggggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                                             <  1:2279783/73‑1 (MQ=17)
                             agagggCTGAAACGGGAGAGCGCCTCCCGAAGCAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGtcat                                             <  2:4150531/73‑1 (MQ=25)
                             aaggggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc                                            <  2:4654580/74‑1 (MQ=17)
                             aaggggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc                                            >  1:5635318/3‑76 (MQ=17)
                             aaggggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc                                            >  1:2327305/3‑76 (MQ=17)
                              gggggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGtcatc                                            <  2:268439/76‑1 (MQ=255)
                              aggggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatca                                           >  2:4882152/2‑76 (MQ=17)
                               ggggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGtcatca                                           <  1:5086807/76‑1 (MQ=255)
                                ggggTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                         <  2:5685444/76‑1 (MQ=14)
                                gggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                         <  1:744285/76‑1 (MQ=17)
                                gggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                         >  1:4063541/1‑76 (MQ=17)
                                gggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                         >  1:1317570/1‑76 (MQ=17)
                                 ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                         <  2:2369382/76‑1 (MQ=255)
                                 ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc                                        >  2:1230714/1‑76 (MQ=17)
                                  gCTGAAACTGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc                                        <  2:4006152/76‑1 (MQ=255)
                                  gCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc                                        <  2:359170/76‑1 (MQ=255)
                                  gCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc                                        <  1:1140675/76‑1 (MQ=255)
                                  gCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt                                       >  1:3266276/1‑76 (MQ=14)
                                   cTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                      <  2:5041497/76‑1 (MQ=17)
                                   cTGAAACCGGAAAGCCCCTCACGAAGAAGGGCGCCTTACAGAAGGAAAGGGTCATGATGAAGCTCGTCATCATACt                                       >  2:3927083/1‑76 (MQ=21)
                                    tGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                      >  2:2931899/1‑76 (MQ=255)
                                    tGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                      >  2:2958681/1‑76 (MQ=255)
                                    tGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGGTATGATGACGCTCGTCATCATACTgg                                     >  2:2583254/1‑75 (MQ=14)
                                      aaaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTggtt                                   <  1:3446698/76‑4 (MQ=14)
                                       caCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGACCTCATACTGGTtg                                  >  1:2844886/2‑76 (MQ=14)
                                       aaggggAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGAtt                                   <  2:5710624/73‑1 (MQ=255)
                                       aaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTtg                                  <  1:4814261/76‑1 (MQ=14)
                                        aCGGGAACGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTtgt                                 <  2:3226906/76‑1 (MQ=14)
                                        aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTtgt                                 <  2:733164/76‑1 (MQ=14)
                                        aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTtgt                                 >  2:4612035/1‑76 (MQ=14)
                                        aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTtgt                                 >  2:3042256/1‑76 (MQ=14)
                                        aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTtgt                                 <  1:2688706/76‑1 (MQ=14)
                                        aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGGTATGATGAAGCTCGTCATCATACTGGTtgt                                 >  1:1735611/1‑76 (MQ=14)
                                           gtaaaGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTgtgt                               <  1:2801182/74‑1 (MQ=255)
                                           ggAAAGCCCCGCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGtgtt                              >  1:5350255/1‑76 (MQ=17)
                                            gAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGtgtt                              <  1:3737645/76‑1 (MQ=255)
                                            gAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGtgtt                              >  2:486716/1‑76 (MQ=255)
                                              aaGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGtgttgt                            <  2:5510514/76‑1 (MQ=255)
                                                 ccccTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGCTATGATGAAGCTCGTCATCATACTGATTGTGTTCTTAc                         >  1:2719737/1‑76 (MQ=255)
                                                  cccTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACtg                       >  1:2655101/1‑75 (MQ=18)
                                                   ccTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACtgt                      >  1:3307005/1‑74 (MQ=18)
                                                     tCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGtgttg                             <  2:558758/68‑1 (MQ=255)
                                                      cccGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTa                    >  2:4765522/1‑76 (MQ=255)
                                                       ccGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTaa                   <  1:4205295/76‑1 (MQ=255)
                                                       ccGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTaa                   <  1:2115325/76‑1 (MQ=255)
                                                         gaagaaGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGt                 >  1:2619494/1‑76 (MQ=255)
                                                         gaagaaGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGt                 >  2:4156343/1‑76 (MQ=255)
                                                          aagaagGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGtt                >  2:4596144/1‑76 (MQ=255)
                                                            gaagGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTc              <  1:361309/76‑1 (MQ=255)
                                                            gaagGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTc              >  1:4298627/1‑76 (MQ=255)
                                                            gaagGGGGCCGTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTc              <  2:5835716/76‑1 (MQ=255)
                                                             aagGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCg             <  2:4217637/76‑1 (MQ=255)
                                                              agGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGc            >  1:5190207/1‑76 (MQ=255)
                                                              agGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGc            >  2:1889576/1‑76 (MQ=255)
                                                              agGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGc            <  1:4313145/76‑1 (MQ=255)
                                                                 gggCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGc         >  2:3600138/1‑76 (MQ=255)
                                                                 gggCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGc         <  1:2566356/76‑1 (MQ=255)
                                                                 gggCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGc         <  1:4358147/76‑1 (MQ=255)
                                                                   gCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCtt       >  2:1992673/1‑76 (MQ=255)
                                                                     cTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCttat     >  1:2014441/1‑76 (MQ=255)
                                                                      ttACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCttatt    >  2:3328698/1‑76 (MQ=255)
                                                                       tACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCttatta   >  2:5193591/1‑76 (MQ=255)
                                                                           |                                                                        
TGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAA  >  NC_000913/3192694‑3192841

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